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Abstract 


The nonlinear coupled boundary value problems (BVPs) have very important and interesting as¬ 
pects in the kingdom of Nonlinear Analysis due to not only the theoretical aspects but also the 
applications which they have in almost every field of science. Problems with coupled boundary 
conditions (BCs) appear while studying mathematical biology, Sturm-Liouville problems, reaction 
diffusion phenomena, chemical systems, and Lokta-Volterra models. 

This thesis has two parts. In the first part, the existence results are established for the first-order 
and the second-order nonlinear coupled BVPs subject to nonlinear coupled BCs. Also in the same 
part, the existence results are established for the second-order nonlinear coupled BVPs when the 
nonlinear functions have dependence on the first-order derivative. 

Multiple approaches are available in the literature to investigate the existence of solutions of 
nonlinear BVPs, but lower and upper solutions (LUSs) approach is one of the strongest. In this 
approach the original problem is modified logically to a new problem, known as modified problem, 
then the theory of differential inequalities with the combination of well-known existence results are 
applied to establish the existence of solution of the modified problem. Finally the solution of the 
modified problem leads to the solution of the original problem. Moreover in the first part of the thesis 
the treatment of the many different first-order and the second-order nonlinear BVPs are unified by 
developing the idea of coupled LUSs. Under this idea, some monotonicity assumptions are imposed 
on the arguments of the nonlinear BCs in the presence of the existence of a lower solution and an 
upper solution to unify the classical existence results for very important types of BVPs, like periodic, 
anti-periodic, Dirichlet, and Neumann. Several examples are discussed to support the theoretical 
results. 

The subject fractional calculus being a generalization of integer-order calculus has numerous 
applications in almost every field of science. Due to the intensive use of fractional order differential 
problems (FODPs) in almost every field of science based including, but not limited to, fluid dynamics, 
physics, aerodynamics, chemistry, mathematical biology, image processing, and psychology, there is 
a strong motivation for the researchers to develop the reliable and efficient numerical methods to 
find the approximate solutions of FODPs. 
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In the second part of the thesis, we consider a generalized class of multi-terms fractional order 
partial differential equations (FOPDEs) and their coupled systems. We develop a new numerical 
method and generalize the corresponding Jacobi operational matrices of integrals and derivatives 
considered on a rectangular plane. By means of the operational matrices, the considered problem of 
fractional order is reduced to an algebraic one. Being easily solvable, the associated algebraic system 
leads to finding the solution of the considered problem of fractional order. Validity of the method 
is established by comparing our Matlab software simulations based obtained results with the exact 
solutions in the literature yielding negligible errors. 
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Preface 


This thesis comprises of two parts and addresses the different types of problems having some very 
interesting and productive impact in the world of nonlinear coupled boundary value problems (CB- 
VPs). 

The first part deals with the existence criterion for the solutions of nonlinear coupled systems of 
ordinary differential equations (ODEs) together with nonlinear coupled boundary conditions (CBCs). 
In the kingdom of differential equations, the study of nonlinear CBVPs have very important and 
interesting aspects for the researchers due to not only the theoretical aspects but also the rich 
applications in the existing literature. The applications of CBCs are observed while studying Sturm- 
Liouvillc problems, mathematical biology, chemical systems, and reaction-diffusion equations [1, 2, 
5, 37, 178, 188]. Also the differential systems with CBCs are studied while dealing with interaction 
problems, reaction diffusion phenomena, and Lokta-Volterra models [6, 117, 128]. 

The study of nonlinear CBVPs inherits complications owing to its completely different behaviours 
of underlying physical phenomena as compared to the uncoupled BVPs. For the sake of establishing 
the existence theory for nonlinear CBVPs a number of techniques are developed. The strongest 
of among is to establish the existence of solutions is the technique of LUSs. In this technique, the 
existence of solutions of the considered problems is established by modifying them logically and then 
the theory of differential inequalities with combination of well known existence results are applied 
to the modified problems. 

In 1893, Picard [145] gave birth to a lower solution while looking for the solution of ODEs of the 
type 


u 2 (t) = gi(t,u 2 (t)), t £ [c, d ], 
u 2 (c) = 0, u 2 (d) = 0. 


( 0 . 0 . 1 ) 


1 
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He built a monotone sequence of successive approximations of the form ai 0 < a\ x < aq 2 < with 
a\ 0 > 0 on (c,d) that converges to the solution of the problem (0.0.1). Further the existence of the 
first approximation aq 0 was proved that satisfies the following inequality 

ot" 0 (t) +gi(t,a lo (t)) > 0, t G (c, d), 
ai 0 ( c ) = °’ ai o (d)=0. 

Later on a± 0 is known by the name of a lower solution. But the major development of the method has 
been observed in 1931 owing to the work of Scorza Dragoni [45]. He discussed the existence results 
for the second-order Diriclilet BVPs of the type u 2 (t) = gi(t,u 2 (t),u 2 (t)), 112 (c) = B\, 112 (d) = 
B 2 , i £ [c, d], and assumed the existence of a lower solution aq and an upper solution /3q, such 
that aq,/?i € C 2 ([c,d\) with aq(f) < fii(t) for all t € [c, d\. The developed method ensured the 
existence of at least one solution u\, such that oq(f) < u±(t) < /3i(t), for all t G [c,d\. Recently, 
many researchers investigated the existence of solutions of BVPs with certain boundary conditions 
(BCs) applying LUSs approach [7, 38, 59, 118, 146, 185]. 

The exchange of letters between Gottfried Leibniz and Marquis de L’ Hopital in 1695 gave birth 
to the idea of fractional calculus. L’ Hopital inquired what would be the consequences if the order 
of differentiation assumed to be non-integer \ instead of n € N. During September 1695, Leibniz 
entertained the question of L’ Hopital with a historical sentence “this is an apparent paradox from 
which, one day, useful consequences will be drawn...” [129]. Many well known mathematicians had 
made considerable development in the field of fractional calculus, Laplace (1812), J.P.J. Fourier 
(1822), N.H. Abel (1823-1826), J. Liouville (1832), B. Riemann (1847), A.K. Grlinwald (1867-1872), 
A.V. Letinkov (1868-1872), Heaviside (1892-1912), J. Hadamard (1892), H. Weyl (1917), Erdelyi 
(1939-1965), H. Kober (1940), and M. Riesz (1949) are the names of few. The contributions of 
Gorenflo and Mainardi can be seen in [71]. 

Integer-order derivatives have clear physical and geometrical interpretations. For instance, the 
first-order derivative gives slope of a tangent line, the second-order derivative elaborates the concavity 
of the differentiable objects. In classical physics the first-order derivative exhibits the velocity of the 
moving objects and the second-order derivative yields the acceleration. It is expected that being 
a generalization of the conventional differential calculus, the fractional calculus should have wider 
meanings, but still it is very miserable in this context. Owing to this deficiency, for a long time, the 
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subject fractional calculus could not get the attention of the researchers working on applied problems, 
they gave it a little respect by considering it nothing but curiosity of the mathematicians. But later 
on they observed the huge advantages of the fractional calculus in almost every field of applied 
sciences. For example in fluid dynamics a considerable advantage was examined when conventional 
Maxwell model [130] was replaced with a fractional order Maxwell model. Fractional order models 
of non-Newtonian fluids played a valuable role in the explanation of viscoelastic properties of non- 
Newtonian fluids [20, 46, 64]. 

Some material is available in the literature that explains the geometrical and physical inter¬ 
pretation of the fractional operators. Trustworthy geometrical and physical interpretation of the 
fractional integration and fractional differentiation is proposed by I. Podlubny [147]. In 1998, Tor- 
bati and Hammond [131] observed the geometrical and physical interpretation of the fractional 
operators in Fractal geometry by considering the fractional operators as linear filters. To explain 
the geometrical interpretation, they used the domino ladder network (a series of capacitors and 
resistors which are connected in various relationship), and Fractal Cantor’s set. 

In classical calculus, derivatives and integrals are uniquely computed. The similar situation 
exists in the case of fractional integrals. For example, many researchers like Marichev, Samko, and 
Kilbas [160], Podlubny [148], Oldham and Spanier [142], Miller and Ross [129] used the similar 
definition to compute the fractional integrals. But the situation is much different and complicated 
in the case of fractional order derivatives because several different competing definitions exist in 
the literature. For instance, a few of those include, the Riemann-Liouville approach, the Caputo 
approach, the Hadamard approach, the Marchaud approach, the Granwald-Letinkov approach, the 
Erdelyi-Kober approach, and the Riesz-Fcller approach. These existing definitions coincide only for 
some particular cases but in general these are not alike. Among of all these approaches to define 
fractional differentiation and fractional integration, the approach of Riemann-Liouville gets more 
attention. But the approach of Riemann-Liouville does not properly address the physics of some 
fractional order problems in the presence of certain boundary and/or initial conditions. Also this 
approach can exhibit the derivative of a constant term other than zero. To overcome this problem 
Caputo proposed an alternate definition of fractional derivative in 1967, and it is used in fluid 
dynamics in 1969 to explain the theory of viscoelasticity. 
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Perhaps, the first application of fractional calculus in physical problems was observed due to 

the work of Niels Henrik Abel (1823) while finding the solution of integer-order equation, known 

as tautochrone problem [8]. In this problem, the curve of an object (frictionless wire), lying in a 

1 

vertical plane, was determined by using the operator D§ , and assuming the dependence of the time 
position not on the starting point. 

According to the lecture notes of Ross [156], the Abel’s solution of the tautochrone problem was 
proved so valuable because it prompted Liouville to do some logical work on the subject fractional 
calculus, then in 1832, the first valid and logical definition proposed by Liouville [119] on fractional 
differentiation was added in the literature. 

During 1832 to 1855, Liouville proposed two different approaches to define the derivatives of 
fractional order. In the first approach the concept of nth order derivative of exponential function 
D n e ax = a n e ax was extended to establish the idea of fractional order differentiation, As a result, 
the following relation was obtained 

D a e ax = a a e ax . (0.0.2) 


Then a function h(x) was expressed in the form of exponential series by using (0.0.2) as 

OO 

D a h{x) = Y M“e a ” x . 

n—0 

This was the first proposed definition of Liouville [157]. But this definition has some limitations 
regarding the values of a, because the series exhibits divergent behaviour for some values of a. In 
[47], Liouville suggested another definition to define fractional order derivative as 


D"/ = ( i rr(« + «) a .-a-a, a > Q 
T(a) 


Later on the use of this definition was seen in potential theory but it also had some limitations, 
because it did not deal with a wide class of functions. In 1847, Riemann proposed the definition 
for integral of fractional order of a function g(x ) by using the complementary solution of nth order 
derivative of integer order as 

D- a g(x) =[ (x-t) a ^g(t)dt + <t>(x). 

r(a) Jb 

In 1880, Cayley commented that the complementary function of Riemann was of indeterminate form. 
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From 1835 to 1850, Mathematicians recorded a rich discussion on the definition of derivative of 
fractional order. George Peacock preferred the generalization of Lacroix instead of Liouville, while 
others gave importance to Liouville’s definition. Augustus De Morgan asserted that both approaches 
dealt with a class of more general system. The debate was settled in the mid of nineteenth century. 
According to Harold Thayer Davis [48], “ The mathematicians at that time were looking for a plausible 
definition of generalized differential but, in fairness to them, one should note they lacked the tools to 
examine the consequences of their definition in the complex plane." 

For many years fractional calculus had been discussed as an abstract concept of mathematics. 
But currently, the applications of fractional calculus can be observed in almost every field of science, 
including physics, bio-chemistry, engineering, image processing, bio-mathematics, and many others. 
It has been realized round the world that the models based on fractional derivative exhibit better 
results than integer-order models. Fractional derivatives inherit nonlocal nature, so it is an excellent 
tool to get better understanding of hereditary properties of different processes and materials. R. L. 
Bagley [21] presented the first Ph.D thesis on the applications of fractional calculus in Viscoelas¬ 
ticity models. Recently, the applications of fractional calculus have been observed in psychology 
to determine the time variation of emotions of mankind [9, 161]. The applications of FODPs can 
be seen in dynamics and control systems [10], Marine sciences and Wave dynamics [143, 149], Heat 
transfer models [72, 179], Diffusion processes [132], Solid mechanics [158, 171], and many more. 

In the literature many monographs and books are available to study the applications and theory 
of fractional calculus. The comprehensive overview and historical aspects of the subject fractional 
calculus have been discussed in the book of K.B. Oldham and J. Spanier [142], The detail study 
on FODPs and derivatives of fractional order is available in the book of K.S. Miller and B. Ross 
[129]. The book of I. Podlubny [147] exhibits the nature and applications of FODPs and fractional 
calculus. In addition many other works on the application of fractional calculus and FODPs have 
been added in the existing literature including the work of A. Carpinteri and F. Mainardi [39], R. 
Hilfer [84], A.A. Kilbas, H.M. Srivastava, J.J. Trujillo [101], J. Sabatier, O.P. Agrawal and J.A. 
Tenreiro Machado [162], and V. Lakshmikantham, S. Leela, and J. Vasundhara Devi [121]. 

The focus of attention of the second part of the thesis is the development of the reliable and 
efficient approximating techniques to compute the approximate solutions of the generalized class of 
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multi-terms FOPDEs and its coupled systems. Currently, a number of numerical techniques are 
available to deal with the approximate solutions of FODPs. The names of the most frequently used 
methods are; Differential transform method [11, 12, 60], Spectral method [22, 24, 49, 50, 51, 52], 
Adomian decomposition method [25, 86], Reproducing Kernal method [73], Homotopy method [100]. 

In this thesis, the Spectral method is employed to compute the approximate solutions of FOPDEs; 
it approximates the solutions in a series form. Further, this method has capability to reduce the 
FOPDEs into a problem of solving a system of algebraic equations which are simple in handling 
by any computational software and indeed makes finding the approximate solutions much easier. 
Different approaches are available to convert the FOPDEs into a system of algebraic equations, 
including collocation approach, Spectral tau approach, Galerkin approach. Galerkin approach is 
similar to Spectral tau approach in the way that differential equation is enforced. However, to 
satisfy the BCs, there is no need of the test functions. Therefore, BCs are applied in the presence 
of a certain supplementary set of equations [133, 150]. While finding the numerical solutions of 
differential equations using collocation approach, firstly, a suitable finite or discrete representation 
of the solutions is obtained using polynomials interpolation based on some proper nodes like well- 
known Gauss-Lobatto nodes. Secondly, the system of algebraic equations is obtained by discretizing 
the considered differential equation. 

Operational matrix method is among the variants of the spectral methods [24, 50, 51]. Under 
this method the fractional order derivatives involved in FODPS are replaced with the operational 
matrices of derivatives of fractional order to obtain the system of algebraic equations. This method 
deals accurately with FODPs when constraints are in the form of some initial conditions. On the 
other hand when fractional order problems have constraints in the form of non-local and local BCs, 
it is not successful. This method is very much effective when the FOPDEs have one spatial variable 
rather than two or three spatial variables. 

The arrangement of the thesis is as follows; In Chapter 1, we discuss a comprehensive introduction 
of coupled lower and upper solutions (CLUSs) that establishes the basis for our further works. We 
also include some results from Fixed Point Theory (FPT), fractional calculus, approximation theory, 
and matrix theory. Section 1.1 deals with the definitions of lower and upper solutions for the first- 
order and the second-order ordinary nonlinear differential systems along with some classical results 
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on the existence of solutions for periodic, anti-periodic, Dirichlet, and Neumann BVPs. Section 
1.2 describes the detail introduction on the development of CLUSs. Moreover some examples are 
discussed in the same section to elaborate the idea of the construction of CLUSs. In Section 1.3, 
a brief note on the development of Nagumo conditions is given. In Section 1.4, we mention the 
development of the fixed point results with a historical note and also mention some basic definitions 
and some results from FPT that play a very important role to establish the existence of solutions 
for our considered nonlinear coupled BVPs. In the remaining part of the Chapter, we recall some 
results from fractional calculus and approximation theory. Also we discuss some properties of the 
orthogonal JPs that are very useful for the development of the numerical results for our later works. 

Chapter 2 deals with the existence results of the first-order nonlinear ordinary coupled BVPs of 
the type 


u[{t) = fi{t,u 2 (t)), te[o, 1], 

u' 2 (t) = f 2 (t,u 1 (t)), t € [0,1], 


(0.0.3) 


subject to nonlinear coupled boundary condition (CBC) 


^(«i(0),U2(0),tt 1 (l),« 2 (l)) = (0,0), (0.0.4) 

where the nonlinear functions f\,f 2 '■ [0,1] xl->l and </> : ]R 4 —> ]R 2 are continuous. In Section 
2.1, we construct the method of lower and upper solutions to establish the existence results. In 
Section 2.3, we ensure the existence of at least one solution (ui,w 2 ) of (0.0.3) such that (ai,a 2 ) ^ 
(wi,« 2 ) ^ (/3i,/? 2 ) by assuming the existence of a lower solution (a l7 a 2 ) and an upper solution 
(/3 i 5 /?2), such that (a 1 , 0 : 2 ) ^ (/3i,/3 2 ) on [0,1] x [0,1]. In Section 2.5, the existence results are 
also discussed when lower and upper solutions are in reverse order, such that (ai,a 2 ) ^ (/3 i,/? 2) 
on [0,1] x [0,1]. Our chosen boundary condition (0.0.4) is very much generalized because it unifies 
the treatment of many different nonlinear and linear first order problems [65, 98, 126]. For instance 
if (/) (x, y, z,w) = {x — C, y — D), then (0.0.4) implies the initial conditions, wi(0) = C, ^ 2 ( 0 ) = 
D , C,D € R. On the same fashion, if <f>{x,y,z,w) = (x~y,z — w), then (0.0.4) implies the 
periodic conditions, «i(0) = «2(0), Mi(l) = ^ 2 ( 1 ). Similarly, if 4> (x,y, z,w) = (x + y,z + w ), then 
(0.0.4) implies the anti-periodic conditions, «i(0) = — «2(0),wi(l) = — ^ 2 ( 1 )- On the same lines, we 
can establish the nonlinear-coupling among the arguments of the boundary condition (0.0.4). To 
satisfy the classical existence criterion for periodic, anti-periodic, and among others nonlinear BVPs, 



we extend the concept of coupled lower and upper solutions explained in [66] for the two dimensional 
case in Sections 2.2 and 2.5.1. Our main tools to establish the existence of solutions of (0.0.3) are 
Arzela-Ascolli theorem and Schauder’s fixed point theorem. In Sections 2.4 and 2.5.3, we check the 
validity of our theoretical results by considering some examples. 

In Chapter 3, we extend the findings of Chapter 2 for the nonlinear ordinary second-order coupled 
BVPs the type 


u'l (t) = fi{t,u 2 (t)), te[o,i], 

u 2 {t) = ie[0,l], 

with nonlinear CBCs 

<t> (ui(0),U2(0),Ul(l),U 2 (l),W , l(0),W2( 0 )) = (0,0), 

^ ^Ui(O), ti 2 (0), ui(l), w 2 (l), Wi(l), = (0,0)) 

and 

-«i (t) = h(t,u 2 (t)), t G [0,1], 

-«2 (t) = f 2 {t,ui(t)), t G [0,1], 

with nonlinear CBCs 

/* (mi(0),m 2 (0),w 1 (0),u 2 (0),u 1 (1),u 2 (1)) = (0,0), 

v (m(o), u 2 (o)) + (m(i), u 2 (i)) = (o, o), 


(0.0.5) 


( 0 . 0 . 6 ) 


(0.0.7) 


( 0 . 0 . 8 ) 


where fi, f 2 '■ [0,1] x R — > R, <f>, ip, fi : ]R 6 —> R 2 , and v : K 2 — > K 2 are continuous functions. The 
existence results for the second-order systems have been discussed by many researchers considering 
coupled and non coupled BCs [14, 15, 40, 74, 86]. In [14], a sequence of modified non-singular 
problem on a compact subset of (0,1) has been constructed and the existence of positive solutions 
of (0.0.7) has been examined by constructing a cone of positive functions, and applying the Guo- 
Krasnosel’skii fixed point theorem. The BCs in [14] are linear and uncoupled. But we examine the 
existence of solutions of (0.0.7) in the presence of nonlinear CBCs. Moreover our BCs (0.0.8) are 
much generalized. For example, if n(j,k,l,m,n,o) = (l — n,m — o) and v(j,k) = (— j,—k ), then 
(0.0.8) implies the periodic BCs and if n(j,k,l,m,n,o) = (l + n,m + 0) and v{j,k) = ( j,k ), then 
(0.0.8) implies the anti-periodic BCs. Definitely, in order to obtain solutions satisfying these BCs 
and lying between a lower solution and an upper solution, some additional conditions are required. 



9 


For example, in the periodic case it suffices that 

MO) > «i(l), a 2 (°) > = ai(0),a 2 (l) = a 2 (0), 

0[{O) < /?i(l), 0>) < 0a(l),0i(l) = 0r(O),0 2 (l) = /?2(0)), (0.0.9) 

and in the anti-periodic case it suffices that 

a'i(O) > -/3i(l), « 2 (0) > -0^(1), -ai(0) = 0i(l), -o 2 (0) = 0 2 (1), 

p'M < -<*[( 1), 0^0) < —a 2 (l),ai(l) = —/?i(0),a 2 (l) = -0 2 (O)). (0.0.10) 

So to ensure the existence of (0.0.9) and (0.0.10), we extend the idea of CLUSs presented in [66] for a 
nonlinear second-order coupled problems applying some monotonicity assumptions on the arguments 
of nonlinear functions /z, and v in Section 3.2.2. The nronotonicity assumptions on the arguments 
of BCs (0.0.8) do not verify the existence criterion for the nonlinear BVPs having Diriclilet and 
Neumann BCs. So, to achieve the existence results for the second-order nonlinear Dirichlet and 
Neumann BVPs, we use the BCs (0.0.6). For instance, if <j>(j,k,l,m,n,o ) = (S — j, W — k) and 
ip(j,k,l,m,n,o ) = (X — l,Y — m) with S,W,X,Y e R, then (0.0.6) implies the Dirichlet BCs. 
Similarly, if <j)(j, k, l, to, n,o) = (n — S,o — W) and ip(j, k, l, to, n, 6) = (X — n,Y — 6), then (0.0.6) 
implies the Neumann BCs. Similarly, we develop the nonlinear-coupling among the arguments of 
BCs (0.0.6) and (0.0.8) and develop the existence of at least one solution (ni,w 2 ) for the nonlinear 
coupled ordinary differential systems of the type (0.0.5) and (0.0.7) in Sections 3.1.3 and 3.2.3. 
In Section 3.1.2, we extend the idea of coupled lower and upper solutions as discussed in [68] for 
the second-order nonlinear coupled differential systems of the type (0.0.5) to verify the existence 
criterion ((0.0.11)-(0.0.12)) for Dirichlet and Neumann nonlinear BVPs under some nronotonicity 
assumption on the arguments of the nonlinear functions <j> and if) 

«i(0) < S < A(0), a 2 (0) < W < fc(0), 

a i(l) < X < /Ml), a 2 (l) < X < /8 2 (1), (0.0.11) 

and 


0i(O) < S < M0), 0^(0) < W < a'(0), 
«i(l)<X<0;(l) } a 2 (l) <Y</3' 2 (1). 


( 0 . 0 . 12 ) 
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In the light of the above discussion, we claim that our results are much generalized because they 
unify the treatment of many nonlinear second-order coupled differential systems, like periodic, anti- 
periodic, Diriclilet and Neumann. In Sections 3.1.4 and 3.2.4, we check the validity of our theoretical 
results by taking some examples. 

In Chapter 4, we extend the results of Chapter 3 by introducing the derivative terms in the 
nonlinear functions f\ and f 2 and establish the existence results for nonlinear coupled systems of 
the type 


-w'i'(f) = fi{t,u 1 (t),u 2 (t),u' 1 (t),u' 2 (t)), tG [0,1], 
~u 2 (t) = h{t,u 1 {t),u 2 (t),u' 1 (t),u 2 (t)), tG [0,1], 


(0.0.13) 


subject to nonlinear CBCs of the type (0.0.6). 

Where /j, / 2 '■ [0,1] xRxRxRxR-fK are continuous functions. Problems with nonlinearity 
dependent on derivative with certain types of BCs, like periodic, anti-periodic, Direchlet, Neumann, 
and multi-points has been discussed by many authors [26, 41, 68, 102, 186, 189]. In [102], the 
existence results for the problems having nonlinearity in first derivative with periodic and multi-point 
BCs were discussed using lower and upper solutions approach, the BCs used in [102] were uncoupled 
and linear. We extend the results for the nonlinear coupled BCs using the same approach. We discuss 
the existence results in the presence of the existence of a lower solution (au, 0 : 2 ) and an upper solution 
(Pi, /3 2 ), which are not necessarily constant, that is, there exist (ou, a 2 ),(/3i, /3 2 ) €E C 2 [0,1] x C 2 [ 0,1] 
such that 


-a'i'(f) < fi (t,a 1 (t),a 2 (.t),a' 1 (t),a' 2 (t) S J , ie [0,1], 
~oi 2 {t) < h (t,a 1 (t),a 2 (t),a' 1 (t),a 2 (t)) , tG [0,1], 
and 

-Pi (t) > fi (t,Pi(t),P 2 (t), P'i(t),P' 2 (t)^ , t G [0,1], 
- P 2 {t ) > f-2 (t,Pi(t),p 2 (t),p' 1 (t),p' 2 (t)j , t G [0,1], 


(0.0.14) 


(0.0.15) 


with (ai,a 2 ) ^ (Pi,P 2 ) on [0,1] x [0,1], and for a more generalized class of BCs (0.0.6), that unify 
the treatment of many second-order nonlinear coupled systems of the type (0.0.13) applying some 
monotonicity assumptions on the arguments of nonlinear functions </>, and ip, as discussed in Chapter 
3. We develop the method of CLUSs for the nonlinear CBVPs by extending the idea presented in 
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[68] and ensure the existence of at least one solution (iq, u 2 ), such that (aq, a 2 ) ^ («i, u 2 ) A (/?i, P 2 ) 
on [0,1] x [0,1] assuming the nonlinear functions /1 and f 2 satisfies the Nagurno conditions relative 
to (ai,« 2 ), (/?i, /S 2 ) in Section 4.1.4. In section 4.1.5, we check the validity of the theoretical results 
by considering some examples. 

In Chapter 5, we develop a reliable and an efficient numerical scheme for getting the approximate 
solutions of multi-terms FOPDEs and their coupled systems of the type 


d ai U(x,y) _ d 11 U(x,y) d 3l U(x,y) 

dx^ ~ ai dy^ + ° 2 dx<*/*dy<n/* +Gl ^ y) ’ 

corresponding to the following initial conditions with j = 0, 1 ,..., n 

UW(0,y) = H j (y), , 

and a coupled system of a generalized class of FOPDEs 

d ai U(x,y) d 11 U{x,y) d l2 V(x,y) d l3 V{x ,y) 


dx ai 


=a 1 


a 2 - 


dx T' 2 


a 3 - 


dy^ 


dy^ 

d 3i U{x,y) d es V(x,y) 

4 dx ei / 2 dy ei / 2 5 dx B 3 / 2 dy e3 / 2 1 ’ ’ 

d a2 V(x,y) d pl V{x,y) d p2 U(x,y) d p3 U(x,y) 
=01 7; -r 0 2 —-f 0 3 


dx a2 


dy pl 

d ei U(x,y) 


dx P2 
d B2 V(x,y) 


4 dx ei / 2 dy ei / 2 5 dx e2 / 2 dy S2 / 2 
corresponding to the following initial conditions with j = 0,1, ...,n, 


dy p i 

+ G 2 (x, y), 


(0.0.16) 


(0.0.17) 


(0.0.18) 


UV\0,y) = H j (y),vV\0,y) = G j (y), 


(0.0.19) 


where n < a\, a 2 < n+1, aj,bj, j = 1,..., 5 are all real constants and G\{x, y), G 2 (x, y), U(x, y), V ( x , y) € 
cm a] x [0, A]). Our developed numerical scheme is based on the operational matrices of integrals 
and derivatives of fractional order for a particular kind of polynomials, namely, the two-parametric 
orthogonal JPs. In Section 5.1.1, we exhibit the notion of two dimensional shifted orthogonal JPs 
and develop its orthogonality relationship with the aid of a weight function. In Section 5.1.3, we 
develop the new operational matrices of integrals and derivatives of fractional order and then in 
Section 5.2 these operational matrices are applied to (0.0.16) and (0.0.18) to develop a reliable and 
stable numerical technique for a generalized class of multi-terms FOPDEs and its coupled systems. 

We check the validity of our developed numerical technique in Section 5.2.3 by taking several test 
examples. 
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In Section 5.3, we extend the idea to three-dimensional shifted orthogonal JPs and developed its 
orthogonality relationship in Section 5.3.1. Based on the operational matrices as discussed in Chapter 
5, we develop an efficient numerical method for the multi-terms high-dimensional FOPDEs. In Sec¬ 
tion 5.3.4, we discuss the application of the operational matrices for multi-terms three-dimensional 
heat conduction equation having mixed derivative terms of the type 


d^U(t,x,y) 

dt ai 


n d 01 U(t,x,y) ( n d 02 U(t,x,y) 
1 dx * + 2 dyfo 


+ C 3 


d 03 U(t,x,y) 

n ^3 o ^3 

ox 2 oy 2 


+ Gi(t, x, y), 


corresponding to the following initial condition 


( 0 . 0 . 20 ) 


U{0,x,y)=f(x,y). (0.0.21) 

where C\,C- 2 and C 3 are constants, 0 < ay < 1, and t € [0, A], x G [0, A], y G [0, A]. 

A significant number of test problems are chosen to check the validity of the developed algorithms 
included in this thesis. Exact and approximate solutions of the considered fractional order problems 
are compared to check the accuracy of the proposed methods. We also develop the analytical 
relations to attain the error bounds of approximate solutions. The convergence of the proposed 
methods is studied in the form of tables and plots for each test problem. 



Chapter 1 

Preliminaries 

Throughout this thesis, we are interested in the development of the existence results of nonlinear 
CBVPs corresponding to nonlinear CBCs applying CLUSs approach. In the same thesis, we develop 
the numerical technique based on the operational matrices of integrals and derivatives of fractional 
order to treat with the generalized class of multi-terms FOPDEs and its coupled systems. The 
analytical results related to fractional calculus are taken from [129, 142, 148]. Let us begins with 
the idea of lower and upper solutions method. 

1.1 Lower and upper solutions method 

While studying the nonlinear CBVPs for the first-order and the second-order ordinary differential 
systems (ODSs) of the type 

u[{t) = fi(t,u 2 (t)), t G [0,1], 

( 1 . 1 . 1 ) 

u 2 (t) = f2{t,ui(t)), te[o,i], 

-u'l (t) = fi{t,u 2 (t)), t G [0,1], 

( 1 - 1 - 2 ) 

-u 2 (t) = f 2 (t,u 1 (t)), t G [0,1], 

= /l(t,Ul(t),U2(t),%(0> u 2(i))> te [0,1], 

a- 1 - 3 ) 

-u 2 (t) = f 2 (t,U l(t),«2(*),Wl(t),W2(*)). t G [0, 1], 
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subject to certain linear or nonlinear CBCs on the compact interval [0,1] C R, the properties of 
LUSs for (1.1.1),(1.1.2), and (1.1.3) are often used to investigate the existence of solutions. In LUSs 
approach, the existence of solutions of the considered problems is established by modifying them 
logically and then the theory of differential inequalities with combination of well known existence 
results are applied to the modified problems. The initial development of LUSs method for ODEs 
was observed in 1893 in the works of Picard [145]. But the major development of the method had 
been observed in 1931 owing to the research article [45] of Dragoni in which the existence results of 
the second-order Dirichlet BVPs were discussed. 

LUSs method ensures the existence of at least one solution (ui,u 2 ), such that (/3i, /3 2 ) R («i, u 2 ) R 
(oq,a 2 ), satisfying (aq,a 2 ) A (/3i,f3 2 ) on [0,1] x [0,1]. The case, when LUSs are in reverse order 
has also been discussed by many authors [66, 103, 104]. Recently, the author [17] investigated the 
existence of solutions of the first-order nonlinear coupled system subject to generalized nonlinear 
CBC, when LUSs are in reverse order, such that (/3i,/? 2 ) ^ (aq,a 2 ) on [0,1] x [0,1]. Moreover in 
Chapter 2, the existence results for (1.1.1) are also discussed when LUSs are in reverse order, that 
is, (oq, a 2 ) >: (/3i, /? 2 on [0,1] x [0,1]. Let us note the definition of LUSs 

Definition 1.1.1. [Lower and upper solutions for the first-order nonlinear coupled sys¬ 
tems] Assume Si and / 2 are continuous functions on [0,1] x R. The order pair of functions 
(ctq, a 2 ), (/3i, f3 2 ) € C 1 [0,1] x C' 1 [0,1] are known to be LUSs of (1.1.1) respectively, if they sat¬ 
isfy the following set of inequalities 

< Si ( t,a 2 (t )), 

a' 2 (t) < f 2 (i, o:i(t)), 

(1.1.4) 

> Si , 

P 2 (t) > Si , 

for all t € [0,1]. If the inequalities in (1.1.4) are strict then the order pair of functions (aq, a 2 ), 
and (/3 i,/ 3 2 ) are called strict LUSs respectively. 

Definition 1.1.2. [Lower and upper solutions for the second-order nonlinear coupled 

systems] Assume Si an d S 2 are continuous functions on [0,1] x R. The order pair of functions 
(oq, a 2 ), (/3i, /3 2 ) € C 2 [0,1] x C 2 [0,1] are known to be LUSs of (1-1.2) respectively, if they satisfy 
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the following set of inequalities 

-<*i (t) < fi (t,a 2 {t)), 

-«2 (t) < h , 

(1.1.5) 

-Pi (t) > fl (t,p2{t)) , 

-P 2 (t) > h (t,Pi(t )), 

for all t G [0,1]. If the inequalities in (1.1.5) are strict then the order pair of functions (aq, a 2 ), 
and (01,02) are called strict LUSs respectively. 


Definition 1.1.3. [Lower and upper solutions for the second-order nonlinear coupled 
systems with nonlinearity depending on the first derivative] Assume fi and f 2 are contin¬ 
uous functions on [0,1] x R 4 . The order pair of functions (oti,a 2 ), (0i,02) G C' 2 [0,1] x C 2 { 0,1] are 
known to be LUSs of (1.1.3) respectively, if they satisfy the following set of inequalities 


-ai (t) < fi (t,ai(t),a 2 (t),a 1 (t),a 2 (t)^ , 
~a 2 (t) < /2 (t, oi(t), a 2 (t), aq(f), a 2 (^)) > 
-0i (t) > fi {t,0i(t),0 2 (t),0[(t),0 2 (t)^j , 
-0 2 (t) > fi (t,0i(t),0 2 (t), 0[(t), 0 2 (t)^J , 


( 1 . 1 . 6 ) 


for all t G [0,1]. If the inequalities in (1.1.6) are strict then the order pair of functions (ai,a 2 ), 
and (0i,02) are called strict LUSs respectively. 

We require certain relationships between BCs and pair of LUSs (( 01 , 02 ), (0i,02)) to build up 
the specific connection at the end points 0 and 1 of [0,1]. In the light of well known classical results 
these relationships has the following form 


• In case of the first-order nonlinear coupled system for the periodic BCs, («i( 0 ) = rti(l), 1(2(0) = 1(2(1)), 
a pair of functions ( 01 , 02 ) and (0i,0 2 ) satisfies 

ai(0) < oi(l), o 2 (0) < o 2 (l) 

0i(O)>0i(l), 02(0) > 02(1). 

• Similarly, for the anti-periodic BCs («i( 0 ) = — Ki(l), 112(6) = —u 2 (l)), a pair of functions 



16 

( ai,a 2 ) and (/?i,/ 3 2 ) satisfies 

«i(0) < —/3i(l), 0 : 2 ( 0 ) < — /?2(1) 

/?1 (0) > — Oi(l), 02 ( 0 ) > — 02 ( 1 ). 

• In case of the second-order nonlinear coupled system for the periodic BCs, Ui (0) = u\ (1), u[ (0) = 

and w 2 (0) = u 2 (l),u 2 (0) = w 2 (l), a P aa ' °f functions (oi,a 2 ) and (/?i,/ 3 2 ) satisfies 

Oi(0) > Oi(l), o 2 (0) > o 2 (l), oi(l) = oi(0), o 2 (l) = o 2 (0), 

P'i{0) < /?i(l), /3 2 (0) < P' 2 {1),Pi( 1) = Pi (0),/? 2 (l) = /3 2 (0)), 

• on the same fashion, for the anti-periodic BCs, «i(0) = —ui(l),u , 1 (0) = —u[( 1), and u 2 (0) = 
—« 2 (1 ),u 2 (0) = —u 2 (l), a pair of functions (oi,a 2 ) and (/3 i,/ 3 2 ) satisfies 

o'i(0) > -A(l), o 2 (0) > -or(0) = ft(l),-a 2 (0) = &(1), 

ft(0) < -o'^l), ^(0) < -0^(1),0,(1) = —/M0 ),o 2 (1) = -/3 2 (0)). 

• For the Dirichlet BCs, (mi(0) = S, Ui(l) = X, u 2 (0) = IF, u 2 ( 1) = F), X, Y, S, IF G K, we have 

oi(0) < S < A(0), o 2 (0) < IF < /3 2 (0), 
oi( 1) < X < ft(l), o 2 (l) < Y < /3 2 (1), 

• similarly, for the Neumann BCs, (^4(0) = S, Wi(l) = X, u 2 (0) = IF, « 2 (1) = F), X, F, S', IF G 
K, we have 

Pi(0) <S< «;(0), /&0) < IF < O 2 (0),. 
ai(l)<X<^i(l), o 2 (l) < F < /? 2 (1). 

1.2 Coupled lower and upper solutions 

In this section, we elaborate the idea of coupled LUSs for the first-order and the second-order BVPs 
subject to some specific nonlinear BCs. We discuss how the ida of couple LUSs makes possible for 
us to unify the treatment of many different first-order and the second-order BVPs with periodic, 
antiperiodic, Dirichlet, and Neumann BCs. We also discuss the fundamental results on which the 
main results of our thesis depend. Let us note the definitions of coupled LUSs. 
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Definition 1.2.1. [Coupled LUSs for the first-order BVPs][66] Assume /i be a nonlinear 
continuous function on [0,1] x R, and ^ be a nonlinear continuous function on R 2 . The functions 
ai,/?i € C' 1 [0,1] are known to be coupled LUSs for 


wiW = /i(<,«i(t))> te[o,i], (i- 2 - 1 ) 

with nonlinear boundary condition (BC) of the type 

4>(ui( 0),ur(l)) = 0, (l- 2 - 2 ) 

if a.\ is a lower solution and /3\ is an upper solution for (1.2.1), such that a.\ < f3\ on [0,1], satisfying 

max{0(ai(O),ai(l)) ,0(ai(O),/3i(l))} < 0 < nrin{^ (^i(O),/?i(l)), 0 (/3i(0), ai(l))}. (1.2.3) 

Remark 1.2.1. It is important to note that the BC (1.2.2) is very much generalized. For example, if 
</> (j, k) = j — k, then (1.2.2) yields periodic BCs, (ui(0) = wi(l)), and if (j) (j, k) = j + k, then (1.2.2) 
implies anti-periodic BCs, (ui(0) = —rti(l)). 

Remark 1.2.2. We note that (1.2.3) generalizes the classical results for periodic and anti-periodic 
BCs by applying the monotonicity assumptions on the arguments of (1.2.2). For instance, if (f) is 
monotone decreasing in the second argument then (1.2.3) satisfies the classical existence results, 
ai(0) < ai(l), and /?i(0) > /3i(l) for the periodic BVPs. On the same fashion, if (f) is monotone 
increasing in the second arguments then (1.2.3) verifies the classical existence results, ai(0) < 
—/3i(l), and /?i(0) > —«i(l) for the anti-periodic BVPs. 

Remark 1.2.3. Coupled LUSs for the problem (1.2.1)-(1.2.2) can be also defined when LUSs are in 
reverse order, such that, ot\ > /?i on [0,1]. In this case, (1.2.3) becomes 

max{<^(ai(0),ai(l)) ,0(/3i(O),ai(l))} < 0 < min{^ (/3i(0), /?i(l)), 4> (ai(0), /?i(l))}. (l- 2 -4) 

(1.2.4) also generalizes the classical existence results for the periodic and anti-periodic BVPs by 
applying the monotonicity assumptions on the arguments of (1.2.2), when LUSs are in reverse 
order, such that, a,\ > /?i on [0,1]. 

Example 1.2.1. Consider the following boundary value problem (BVP) 


u^ft) = —2ui(t) + Asin(lOf), t € [0,1], 


(1.2.5) 



18 


with two-point nonlinear boundary condition (NLBC), defined as 

0(u 1 (O),u 1 (l))=u?(O)-u?(l), (1.2.6) 

where X is a positive integer. Let a\(t) = —2A, and fii(t) = 2A, with cti(l) < fiift), for all t £ [0,1]. 
Obviously, — 2A, and 2A are LUSs respectively of (1.2.5), because 

0 ^( 1 ) = 0 < fi(t,ai(t)) = —2(c*i (t)) + A sin (101), t £ [0,1] 

= 4A + A sin (101), 

(1.2.7) 

/3[(t) = 0 > f\{t, fii (1)) = -2(/3i(1)) + A sin (101), t £ [0,1] 

= —4A + A sin (101). 

Moreover, in the light of (1.2.3) and (1.2.7), —2A, and 2A are coupled LUSs respectively of (1.2.5)- 
(1.2.6), because, if f is monotone increasing in the second argument, then (1.2.3) yields 

<A(a 1 (0),A(l))=a?(0)-^(l) 

= —8A 3 - 8A 3 

< 0 , 

<A(/3 1 (0),a 1 (l))=^ 3 (0)-a 3 (l) 

= 8A 3 + 8A 3 
> 0 . 

Similarly, 

iff) is monotone decreasing in the second argument, then (1.2.3) implies 

4> ( a i(0), «i(l)) = a 3 (0) - a 3 (l) 

= -8A 3 + 8A 3 

= 0 , 

0(/3 1 (O),/3 1 (l))=/3 3 (O)-/3 3 (l) 

= 8A 3 - 8A 3 
= 0. 

Example 1.2.2. Consider the following BVP 


u[(t ) = 4 u\{t) + A cos (51), 1 £ [0,1], 


(1.2.8) 
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with two-point NLBC defined as 

0 (ui(O),ui(l)) = wi(l)iti(0) - u?(0)ui(l), 


(1.2.9) 


where X is a positive integer. Let a\(t) = 3A, and fii(t) = —3A, with a\(t) > fiiit), for all t £ [0,1]. 
Obviously, 3A, and — 3A are the LUSs respectively of (1.2.8), because 


ai'ifit) = 0 < fi(t, «i(i)) = 4 (ai(t)) 5 + A cos (5 1), t £ [0,1] 
= 972A 5 + A cos (5 1), 

P[{t) = 0 > fi(t,/3i(t)) = 4(/3i(l)) 5 + A cos (5 1), t £ [0,1] 
= —972A 5 + A cos (51). 


( 1 . 2 . 10 ) 


Moreover, in the light of (1-2.4) and (1.2.10), 3A, and — 3A are coupled LUSs respectively of (1.2.8)- 
(1.2.9), because, if <j> is monotone increasing in the first argument, then (1.2.4) yields 


4 > (c>!i(0), «i(l)) = ai(l)ai(0) - ai(0)ai(l) 


= 9A 2 - 243A 5 


< 0 , 

<t> (/3i(0), /?i(l)) = ft(l)ft(0) - #( 0)h 
= 9A 2 + 243A 5 
> 0 . 

Similarly, 

if <t> is monotone decreasing in the first argument, then (1.2.4) implies 

(MO), a,(l)) = ai (l)/3i(0) - (0)ai(l) 

= —9A 2 - 243A 5 

< 0 , 

<^( Ql (0),ft(l)) = /3i(l)ai(0) - a?(0)ft(l) 

= -9A 2 + 243A 5 

> 0 . 

Definition 1.2.2. [Coupled LUSs for the second-order BVPs][68] Assume /i be a nonlinear 
continuous function on [0,1] x R, and </>,'!/’ are the nonlinear continuous functions on M 3 . The func¬ 
tions ai,fii £ C 2 [0,1] are known to be coupled LUSs for 
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u"{t) = fi (t,U!{t )), ie[0,1], (1.2.11) 

with nonlinear BCs 


^(ui(0),ui(l),w / 1 (0)) = 0, 
^(tii(0),ui(l), «i(l)) = 0, 


( 1 . 2 . 12 ) 


if a\ is a lower solution and /3\ is an upper solution for (1.2.11), with a\ < (3\ on [0,1], satisfying 


max{</> (/?i(0), /3r(l), /3((0)), <p (/3i(0), ou(l), /3((0))} < 0 
min{^(ai(0),ar(l),ai(0)) ,<()(ai(0),/3i(l),a , 1 (0))} > 0, 
and (1.2.13) 

nraxO (/?r(0), /?r(l), P[(l )), ip (oi(0), /3i(l), p[{l))} < 0 
min{0(o; 1 (O),ai(l),a , 1 (l)) ,V’(A(0),ai(l),a , 1 (l))} > 0. 

Remark 1.2.4. It is worth mentioning that (1.2.12) generalizes most of the usual linear and nonlinear 
BCs. For example, if cp ( j , k, l) = C — j, and ^ ( j , k,l) = D — k, then (1.2.12) gives Dirichlet BCs, 
(«i(0) = C, Ui(l) = D). Similarly, if (j) (j , k, l) = l — C, and if) ( j , k, l) = D — l, then (1.2.12) gives 
Neumann BCs, (wi(0) = C,u' 1 ( 1) = D). 

Remark 1.2.5. We observe that (1.2.13) generalizes the classical existence results for Dirichlet and 
Neumann BVPs by applying the monotonicity assumptions on the arguments of (1.2.12). For 
instance, if </> is monotone increasing in the second argument and ip is monotone increasing in 
the first argument then (1.2.13) satisfies the classical existence results, «i(0) < C < /3r(0), and 
ai(l) < D < /3i(l) for the Dirichlet BVPs. On the same fashion, if (p is monotone increasing in 
the second argument and ip is monotone increasing in the first argument then (1.2.13) satisfies the 
classical existence results, /3((0) < C < 04 ( 0 ), and 0^(1) < D < /3((1) for the Neumann BVPs. 

Remark 1.2.6. (1.2.12) and (1.2.13) do not allow us to generalize the classical existence results for 
periodic and anti-periodic BVPs. But if we replace (1.2.12) with the following nonlinear BCs 


^(ui(0),'u , 1 (0),m' 1 (1)) = 0, 
ip (ui(0)) +«i(l) = 0, 


(1.2.14) 
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then (1.2.14) allows us to verify the existence criterion for periodic and anti-periodic BVPs under 
the monotonicity assumptions on the arguments of (1.2.14). For this case, (1.2.13) becomes 


max{<()(/3 1 (O),^(O),/3((l)),0(A(O),^(O),a / i(l))}<O 
min{</> («i(0), 0 ^( 0 ), (a)(1)) , (f> («i(0), 0^(0), /3((1))} > 0, 
and (1-2.15) 

max{V> (/3i(0)),V’(ai(0))} + «i(l) = 0 
min{'i/’ (ai(0)),^(/3i(0))} + /3i(l) = 0. 

Example 1.2.3. Consider the following BVP 

u" (t) = uf (t) + cos (t), t e [0,1], (1.2.16) 


with two-point nonlinear BCs defined as 

<t> (iri(O), «i(l), Wi(0)) = ui(0)ui(l) - (wi(l)u , 1 (0)) 3 , 
t/>(ui(0),ui(l),ui(l)) = ui(0)ui(l) - (ui (l)ui(0)) 5 . 


(1.2.17) 


Let «i (t) = sin(i) — 6, and = sin(f) + 6, with a\(t) < fi\{t), for all t G [0,1]. Obviously, 
sin(f) — 6, and sin(f) + 6 are LUSs respectively of (1.2.16), because 


= - sin (t) > fi{t, ai(t)) = (aii(f)) 3 + cos (f), t G [0,1] 

= (sin(f) — 6) 3 + cos(t) 

< — sin(f), 

(1.2.18) 

Pi(t) = - sin (t) < /i(t,/3i(t)) = (/3i(t)) 3 + cos (t) , t G [0,1] 

= (sin(i) + 6) 3 + cos (t) 

> — sin(i). 


Moreover, in the light of (1.2.13) and (1.2.18), sin(t) —6, and sin(f) + 6 are coupled LUSs respectively 
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of (1.2.16)-(1.2.17), because, iff is monotone increasing in the second argument, then (1.2.13) yields 

<H/M0),/M1),#(0)) = /?i(0)/3i(l) ^ (/3i(l)/3i(0)) 3 
= 6sin(l) + 36 — (sin(l) + 6) 3 

- 279.17 
<0, 

0(ai(O),ai(l),o; , 1 (O)) = ai(0)ai(l) - (ai(l)ai(O)) 3 
= 6sin(l) — 36 — (sin(l) — 6) 3 
= 106.32 
> 0. 

Similarly, 

if f is monotone increasing in the first argument, then (1.2.13) implies 

f (/3i(0), /3i(l), /?i(l)) = /3i(0)#(l) - (/3i(l)/3((0)) 5 

= 6cos(l) — (sin(l) + 6)° 

- 14984.89 
<0, 

^(ai(0),ai(l),ai(l)) = ai(0)a' 1 (l) - (ai(l)ai(O)) 5 
= —6 — (sin(l) — 6) 5 
= 3646.83 
> 0. 

Example 1.2.4. Consider the following BVP 

u'((t) = 5 u'i(t) + ( u\(t ) + l) 2 — sin 2 (nt ), t £ [0,1], (1.2.19) 

with two-point nonlinear BCs defined as 

<t> ( u i(0), wi(1),m / 1 (0)) = (ui(0)sin(u!(0)) +cos(w' 1 (0))ui(l)), 

( 1 . 2 . 20 ) 

f [ui (0), u\ (1), u[ (1)) = (ui (0) tan (m (0)) + cos («i(1)) ui (1)). 

Let «i (t) = —t 2 — t, and fd i (t) = t 2 + t, with a\(t) < /?i (t), for alltG [0,1]. Obviously, —t 2 — t, and 
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t 2 +t are the LUSs respectively of (1.2.19), because 

a"(t) = — 2 > fi{t, ai{i),a'i{t)) = 5 a[(t) + {a\(t) + l) 2 — sin 2 (nt ), t € [0,1] 

= —5(21 + 1) + (t 2 + t — l) 2 — sin 2 ( 7 rf) 

= t 4 + 2 1 3 — t 2 — 12 1 — (4 + sin 2 (7 rf)) 

<- 2 , 

/?i (*) = 2 < fi{t,Pi(t),= 50(0) + (/?i(i) + l) 2 - sin 2 ( 7 rf), t € [0,1] 

= 5 (21 + 1) + (t 2 + i + l) 2 — sin 2 ( 7 it) 


( 1 . 2 . 21 ) 


Moreover, in the light of (1.2.13) and (1.2.21), — t 2 — t, and t 2 + t are coupled LUSs respectively of 
(1.2.19)-(1.2.20), because, if <f is monotone increasing in the second argument, then (1.2.13) yields 

<t> ( 0 ,( 0), /Ml), #(0)) = 03,(0) sin 03,(0)) + cos (#(0)) /3,(1)) 

= —2cos(l) 

- 1.0806 


^(ai(0),ai(l),c/,(0)) = (a,(0)sin(a,(0)) + cos (c/,(0)) a,(l)) 
= 2cos(l) 


= 1.0806 


Similarly, 

if 4> is monotone increasing in the first argument, then (1.2.13) implies 

if 08,(0), 0,(1),0((1)) = (0i(0) tan (0,(0)) + cos (0((1)) 0,(1)) 
= 2cos(3) 


- 1.98 



and 
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4> (ai( 0 ), ai(l), 07 ( 1 )) = (07 ( 0 ) tan (07 ( 0 )) + cos ( 07 ( 1 )) «i(l)) 
= —2 cos(3) 

= 1.98 

> 0 . 

We need the following lemma [ 66 ] for the extension of our later work. 


Lemma 1 . 2 . 1 . Let L : C 1 [0,1] —» Cg [0,1] x K. be defined by 

[Lui\(t) = — Mi( 0 ) + A J ui(s)ds,a«i( 0 ) + 6 wi(l)^ , 

where A, a, and b are real constants with A > 0, such that 


(a + be 0 , 


and 


C<}[0, 1] = {w e C^O, 1] : «7(0) = 0} . 
Then L _1 exists and is continuous and defined by 

[L~\y,T )](*) = e~ M A + y(t ) - A f e x ^y(s)ds, 

Jo 

with 

A _ 7 + Ab/g 1 e A(s ~ 1} y(s) ds - by( 1) 
a + be~ x 


Proof. Choose 


and 


y(t) = ui(t) - mi(0) + A [ iti(s)ds, 
Jo 


7 = arti(O) + bui(l). 

In the light of (1.2.24)-(1.2.25), (1.2.22) can also be written as 

[Lui\{t) = (y(t), 7 ). 


Differentiating (1.2.24) w.r.t. t, we have 


( 1 . 2 . 22 ) 

(1.2.23) 

(1.2.24) 

(1.2.25) 

(1.2.26) 


y'(t) = u\(t) + A«i(f). 


(1.2.27) 
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Multiplying (1.2.27) with integrating factor e xt , we have 

e xt y'(t) = (ui (t)e Xt )', 

then after integrating and taking the limits of integration from 0 to t, (1.2.28) becomes 

Ui(t) = Ui(0)e~ xt + y(t) - A f e A(s-t) y(s)ds, 

Jo 

Mi(0) can easily be determined with the help of (1.2.25) as 

7 = (a + be~ x ) u\(0) + by(l) — bX f e x( ' s ~ 1) y(s)ds , 

Jo 

which further implies 

/ \ 7 + bX e x ( s ~^y(s)ds — by(l) , 

ui(0) = --^, a + be~ x ? 0, 


a + be 


for simplicity of notation, let 


7 +tA; 0 V<-'y«)<i.-Mi) o+fc -A ^ 0 

a + be~ x 

Using (1.2.31) in (1.2.29), we have 

u\(t) = Ae~ xt + y(t) — A f e x( - s ~^y(s)ds. 

Jo 

Equation 1.2.26 implies 

[L~ x (: ) = ui(t). 

Hence, (1.2.31) proves the result. 

The extension of Lemma 1.2.1, is the following 
Lemma 1.2.2. Let L : U 2 [0,1] —x Cq[ 0, 1] xlxR be defined by 

\Lufi[(t) = ^u'-l( t) — u\ (0) — A J Ui(s)ds,o«i(0) + 6ui(l),cui(0) + dui(l) 

where A > 0 is a real constant. Moreover, d , c, b and a are also treated as real constants, 

(ad — be) (e~'^ x — e v ^ 0, 

and here 


(1.2.28) 

(1.2.29) 


(1.2.30) 

(1.2.31) 

(1.2.32) 

(1.2.33) 

□ 


(1.2.34) 
such that 


C 0 2 [0,l] = { Wl eC 2 [0,l]:i(0)=0}. 
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Then L 1 exists and is continuous and defined by 

[L~ 1 (y,'y,d)](t) = ^Ae^ t + Be-' rxt + l -j\^ t -^y{s)As 

with 


A 


and 


B 


Proof. Choose 


(ad — be) ( e^ 


^2(5 ^a + be — d(a + be ^ 




) y(s)ds 


J r e Vx{l ~ s) y(s)ds + d(a + be-' /x ) 

bfc + de-^ J e Vx{s - 1 ) y(s)ds S j, 


(ad — bc)[e '/A _ e \/A 


26 (a + be ' /I ) - d (a + fce' 7 * 
J e Vx( - 1 - s) y(s)ds + d{a + be' /x ^ J e^' - 1 )|/(s)ds 

b^c + de^J e Vx{s ~ 1) y(s)d. 


y(t) = u^t) - Ui( 0) - A / ui(s) ds, 


7 = aui(0) + 6ui(l), 


<5 = cui(O) + dui(l), 

In the light of (1.2.35)-(1.2.37), (1.2.34) can also be written as 

[L (wi)](t) = (y(t), 7 ,S). 

After differentiating the Equation 1.2.35, we have 


(1.2.35) 

(1.2.36) 

(1.2.37) 

(1.2.38) 


y'(t) = u'((t) — Xui(t), A > 0. 


(1.2.39) 
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The general solution of (1.2.39) can be easily determined using variation of parameters technique 
along with integration by parts and taking limits of integration from 0 to t, we have 

Ul (t) = Ae^ + Be\ f e' /X(t - s) y(s)ds + \ f e VX( ~ s ~ t) y(s)ds. (1.2.40) 

2 Jo 2 J 0 

A and B can be easily determined with the help of (1.2.36) and (1.2.37) as 


( c 


7 = ( a + be ^ A + (a 


be 


— y/\ 




[ e' /X( 1 - s) 2 /(s)ds + e' /X(s - 1 ) 2 /(s)ds 

Jo 


5 = (c + de^) A+ (c + de v7 *) B + ^ e ^ 1 S ^y(s)ds + e' 7 ^ 7s 1 ) y(s)ds^ 

Solving (1.2.41), we have 


A = 


(ad — be) (t 


,y/X _ p— \[\ 


and 


B = 


^ 2 c> (a + be ' 7 ^) — d (a + be 
J e^ 1 ~ s) y(s)ds + d(a + be^ V ^j J e VX ^- 1) y(s)ds 

bfc + de-^ J e^ s - 1) y(s)ds^, 

■ (^2<5 (a + l 


e 


)- d (° 


(ad — be) ( 

J e^ 1 ~ s) y(s)ds + d(a + be' / ^ e' /X< - s ~ 1) y(s)d. 

27 (c + de^ + b (c + de' 71 ) J e^ 1 ~ s) y(s)ds 
bfc + de^j e' /X(s - 1) y(s)ds^. 


(1.2.41) 


(1.2.42) 


(1.2.38) can also be written as 


ui(t) = [L-l (j/, 7 ,d)](t). (1.2.43) 

Hence (1.2.42) proves the required result. □ 

A fundamental result to ensure the existence of at least one solution U\ of (1.2.1) corresponding 
to generalized nonlinear BCs (1-2.2) under the approach of coupled LUSs is the following (see [ 66 ]) 
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Theorem 1.2.3. Assume that a± and pi are the CLUSs for the BVP (1.2.1)-(1.2.2). In addition, 
suppose that the functions 

4> ai (x) := 0(ai(O),a:), 

: = </H/3i(0),z), 

are monotone on [ai(l),/3i(l)], then (1.2.1)-(1.2.2) has at least one solution, u\ £ [ai,/3i]. 
Example 1.2.5. Consider the following BVP 

u' x {t) = — 2m(t) + Asin(lOi), t 6 [0,1], (1.2.44) 

with two-point NLBC, defined as 

<t> M0), Ul (l)) = u?(0) - u?(l), (1.2.45) 

where X is a positive integer. Let a\(t) = —2A, and fti{t) = 2A, with a\(t) < fiiit), for all t £ [0,1]. 
Obviously, —2A, and 2A are the LUSs respectively of (1.2.44) as already been varified in Example 
1.2.5 

Moreover, it has already been varified in Example 1.2.5, that — 2A, and 2A are the CLUSs respectively 
of (1.2.44)-(1.2.45). 

Further the functions, (x), <j>p 1 (x) : R 2 —> R) are monotone non-decreasing on [aii(l), /?i (1)] = 
[—2A, 2A] because for every x\ < X 2 in [—2A, 2A], we have 

</> ai (xi) ^ 0a 1 (^ 2 ) 

0(ai(O),a;i) ft </>(a i(0),x 2 ) 
xf — af(0) <x\ — Oi(0) 

Xi < x 2 , 

and 

0 ft Oft) ^ 0ft Oft) 

0(A(O),a:i) ft 0(/3i(0),x 2 ) 

x\ - Pf(0) <x\- 

X[ < X 2 . 

Hence according to the Theorem 1.2.3, (1.2.44)-(1.2.45) has at least one solution u\ £ [ai,/?i]. 
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1.3 Nagumo condition 

In this section, we talk about the existence results of a nonlinear BVPs to the case when the 
nonlinear term is of the form /) (t. u-\ ,u\), that is also depends on u\ . Existence of well-ordered 
LUSs are enough to ensure the existence of solutions of a nonlinear BVPs when nonlinear term 
/i(f, ui,u[) is independent of u\ . But the existence of well-ordered LUSs is not sufficient to ensure 
the existence of solutions for the case when nonlinearity f\ also depends on u' 1 . For example, 
u\(t) := 4 — a/4 — t is a solution of a nonlinear differential equation u”(t) = 2{u' 1 (t)) 3 on [0,4). 
Clearly, Ui(t) := 4 — \/4 — t is bounded on [0,4) but its derivative is unbounded on [0,4). This is 
due to the fact that the nonlinear term f 1 (t,Ui,u' 1 ) = 2(u' 1 (t)) 3 grows very rapidly with respect to 
the gradient of the solution. Moreover, in 1954, Nagumo [137] took an example of Dirichlet problem 
and proved that the considered problem had no solution although the ordered lower and upper 
solutions exist. Later on Habets and Pouse [87] extended the results of Nagumo to the periodic 
and separated BCs and proved that the method of LUSs is not valid without considering Nagumo 
conditions. Therefore to ensure the existence of a solution, we need to find a prior bound on the 
gradient of a solution. For this a very fruitful assumption is used, that is, for nonlinear BVPs, 
the nonlinearity fi grows not faster than quadratically with respect to the gradient. Consider the 
following differential equation on the interval (0,1) 

u"{t) = /i(£,ui,«i), £€(0,1). (1.3.1) 

The nonlinear term fi(t, iq, u \) in the presence of Nagumo condition exhibits the maximum growth 
rate with respect to the gradient of the solution that is quadratic. We need to estimate u\ on (0,1) 
to ensure the existence of solution of (1.3.1) in the presence of some BCs. For this we will impose 
very useful growth condition which is the result of Nagumo [138] and Bernstein[28] works, such that, 
for all £ € (0,1) and u[ el, we have 

|/i(i,ui,ui)| < c(«i)(l + |Mi] 2 ). (1.3.2) 

To observe that (1.3.2) gives boundness of lit'/, assuming boundness of |wi|. For this assume |ui(£)| < 
N, for all t € (0,1), sq be the extremum point of u\(t) and Si € (0,1) with sq < si- If u'i(t) > 0 on 
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[s 0 ,«i], then in the light of (1.3.2), we have 

” < c (^(t))K(t)), t€[8 o,ai], (1.3.3) 

which further implies that 

^ ln (l + uf(t)) - ^C(ui) < 0, t, G [s 0 ,si], (1.3.4) 

where C(u\) = f^ 1 c(s)ds. So, \ ln(l + — C(ui) is decreasing on [s 0) si]. Hence 

^ ln(l + uf(t)) < C{m) - C(ui(s 0 ). (1.3.5) 

(1.3.5) implies that u[ is bounded on [so,Si]. If u\ (t) < 0 on [s2,S3], where 0 < S2 < S3 < 1. Then 
(1.3.2) implies 

<(i)>-c( Ul (t))( 1 + K(t)| 2 ). (1.3.6) 

Hence 

u"(t)wi(f) < -c(tii(t))(l + |ui(i)| 2 )«i(t). t G [s 2 ,s 3 ]. (1.3.7) 

Hence u\ is bounded on (0,1). The important point is that this argument does not depend on the 
BCs. If Nagumo condition is violated then u\ need not be bounded. For example, the following 
BVP [105], 

+ (1 + ( u i(t)) 2 ) 3 ^ 2 = t G (0,2), 
ui(0) = 0, u 1 (2) = 0, 
has a solution the upper half of the circle 

(i-l) 2 +«? = l, 

having infinite derivatives at t = 0 and t = 2. 

Now we give the general definition of Nagumo condition relative to upper and lower solutions. 

Definition 1.3.1. (see, [27]) Suppose /1 : [0,1] xR 2 ->K and ai,/?i : [0,1] —> R are continuous 
functions provided (t) < on [0,1], and suppose that there is a function g : [0, 00) — > (0,oo), 
such that 


\fl(t,Ui,Ui)\ < ffdu'il) 


(1.3.8) 
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for all t G [0,1], ai(t) < ui < fii(t), u[ G R. with 

r°° s ds 

/ —— > max fdAt) — min a\(t), (1.3.9) 

Jm g( s ) o <t<i o<t<i 

where 

M := max{|/3i(l) - «i(0)|, |ai(l) - /3i(0)|}, (1.3.10) 

where function g is known as Nagumo function . 

To observe that Nagumo condition yields boundness to the derivative u\ of the solution u\ of 
(1.3.1), the following theorem is required [27, 102]. 

Theorem 1.3.1. Suppose that fi : [0,1] xl 2 -) 1 satisfies a Nagumo condition on [0,1] with 
respect to aq,/3i with aq(t) < fii(t). Assume Ui G C 2 ([0,1],K.) is a solution of (1.3.1) satisfying 
a\(t) < Ui(t) < for all t G [0,1]. Then for any u\ G C 2 ([0,1],K), there exists a positive 

constant Mi that depends on aq,/3i and Nagumo function g only, such that 

\u'i(t)\<Mi, t G [0,1]. (1.3.11) 


Proof. Choose Mi > M due to (1.3.9) so that 


r Ml sds 
m d(s) 


> max BAt) 
0<t<l 


min oq (t). 
0<t<l 


(1.3.12) 


Let ui be the solution of u'{ = fi{t,u i,wi) on [0,1] satisfying oq(t) < ui(t) < fii(t). Then according 
to mean value theorem there exists s 0 € [0,1], such that 


|wi(s 0 )| = |mi( 1) - wi(0)| < M. 


Suppose on contrary (t)| > Mi. Then there exists [si,S2] C [0,1] such that the following cases 

hold: 


1. «i(si) = M,u' 1 (s 2 ) = M l5 and M < u^ft) < Mi, t G (si,s 2 ), 

2. u'^si) = Mi,u / 1 (s 2 ) = M, and M < u[(t) < Mi, t G (si,s 2 ), 

3. «i( s i) = -M, m' 1 (s 2 ) = -Mi,and — Mi < u'i(t ) < — M, t G (si,s 2 ), 

4. Wi(si) = -Mi,w' 1 (s 2 ) = — M, and - Mi < u'i(t) < —M, t G (si,s 2 ). 
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Firstly, we consider case(l) on [si,s 2 ] and in the light of (1.3.8), we can write 


5(KW I) 
<(*KW 


< Wi(t), 


sWW) 
<( f K(0 
sWW) 


dt < 


Mi(t)dt, 


dt < mi(s 2 ) — ui(si), 


which further implies that 


u i(® 2 ) sc [ s 

~rx fs ui(s 2 ) — Ui(si) < max /?i(i) — min ai(t), 

K-0 51 s ) 


0<t<l 


0<t<l 


it further follows that 


r Ml sds_ 
m g{s) 


< Wi(s 2 ) 


Ui(si) < max 6i(t) — min «i(i), 

0<t<l 0<t<l 


which contradicts to (1.3.12). 


Now consider case (2) and again with the help of (1.3.8), it can be written as 


u '( OKW 


sK(*)) 
rs2 /<(iM(i) 


<«i(i), i € [si,s 2 ] 


rS2 


sKW) 

sWW) 


dt < 


j,i(t)dt, 


dt < ui(s 2 ) — Mi(si), 


which further implies that 




min ai(t). 

0<t<l 


it further follows that 


Jm 1 d( s ) 


min ai(t), 

0<t<l 


which again contradicts to (1.3.12). 
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Now, consider case (4), and again using (1.3.8), we have 

\u"(t )| = |/i(t,ui,ui)| < g{\u[(t)\), t€ [si,s 2 ], 
\u"(t)\ui(t) = |/i(t,ui,wi)| > ff(|ui(i)|)wi(i). 


it further implies that 


^(KWI) 
<( f KW 

s(- u 'i(*)) 


> (t), 




< Wi(si) - Wr(s 2 ), 


which further leads to 


—sds 


_ sds _ _ 

9{s) Jmx g(s) 


< Ml(si) - Wl(s 2 ), 


Finally, we have 


f Ml sds 

/ —^ < rti(si) — Ui(s 2 ) < max 3\(t) — nrin a\(t). 

Jm g(s) o<t<i o<i<r 


Again there is a contradiction. Case (3) can be dealt is a similar way. Hence, it can be concluded 
that (1.3.11) is valid. □ 


1.4 Fixed Point Theory 

Fixed point theory (FPT) plays an important role in establishing the classical results of modern 
mathematics. It is a gorgeous combination of geometry, analysis, and topology. Theorems relating 
to the properties and existence of fixed points are identified with the name of fixed point theorems. 
Fixed point theorems has numerous applications in almost every field of science, the names of few 
are, chemistry, physics, engineering, economics, game theory, biology, medical science, and much 
more, see for example [30, 32, 79, 134, 172, 190]. In this Section, the short note on the development 
of FPT along with some fixed point results is presented. 
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1.4.1 Development of Fixed Point Theory 

Assume /i ba a mapping on a nonempty set X such that /i : X —> X, if fi(x) = x, then the 
point x £ X is known to be a fixed point (FP) of the map f±. On the other hand, if /i ba a 
multi-valued map, i.e., f\ : X —> 2 X , and x € fi{x) then the point x G X is known to be a FP of 
the map /i- The main advantage of the FPT is the ability to transform the wide class of equations 
occurring in a variety of physical formulations into a fixed point inclusions or equations. In 1837, 
Liouville [124] gave birth to the idea of FPT by developing the theory of successive approximations 
used to investigate the existence of solutions of DEs. The other notable contribution regarding 
the existence of solutions of DEs applying FPT was due to the work of Picard [151] in 1890. But 
the major breakthrough in the held of nonlinear analysis was become possible in twentieth century 
due to the pioneering work of well renowned mathematician Stefan Banch [33]. He gave the very 
important result, known as Banach fixed point theorem, which states that “a contraction mapping 
on a complete metric space has a unique fixed point’’. In 1912, Brouwer [34] presented a very strong 
result that led the foundations of classical FPT, that states that “a continuous map on a closed unit 
ball in K" has a fixed point”. In 1930, this idea was extended in Schauders fixed point theorem, that 
states that “a continuous map on a convex compact subspace of a Banach space has a fixed point” 
[165]. Afterward a number of researchers, for example (Tychonoff [173], Kakutani [108], Lefsclretz 
[125], Tarski [174], Edelstein [63], Kannan [109], Clratterjea [43]) enriched the FPT by extending 
these distinguished results. Banach emphasized the the map /i should be continuous on space X but 
in 1968, Kannan [109] presented the fixed point theorem for operators that need not be continuous. 

1.4.2 Some basic definitions and known results 

Some known results form functional analysis and basic definitions is given in this section. For 
detailed study the reader can consult [54, 89]. 

Definition 1.4.1. [Compact] A subset X of a Banach space © is said to be compact if every 
sequence in X has a subsequence that converges to an element contained in X. 

Definition 1.4.2. [Relatively Compact] A subset X of a Banach space © is said to be relatively 
compact if and only if its closure X is compact. 
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Definition 1.4.3. [Compact map] Assume X be a subset of a Banach space © and fi : X —> ©. 
If /i maps every bounded subset of X into a relatively compact subset of © then /i : X — > © is 
compact. 

If /i : X —> © is continuous and compact then it is completely continuous map. 

Theorem 1.4.1. [Arzela-Ascoli Theorem] If 13 be the subset o/M", C(15) represents the Banach 
space of continuous functions and X be a subset of C(15) then X is said to be relatively compact if 
and only if it is equibounded and equicontinuous. 


1.5 Fractional Calculus 

In this section, we discuss some preliminary concepts of fractional calculus, some results from matrix 
theory, and some results related to approximation theory. We give the general idea of fractional 
integrals and derivatives with the aid of Euler Gamma function. 


1.5.1 Fractional order differentiation and integration 

The generalization of factorial function to Gamma function leads to the concept of derivatives of 
fractional order. 


Definition 1.5.1. [148] In 1729, Euler developed the idea of Gamma function while working on the 
interpolation problem for factorial function. It is defined as 


cOO 

T(x)= t x ~ 1 e~ t dt, ieR\{0,-l,-2, •••}. 
Jo 


The Cauchy’s integral [159] formula for n-folcl integral assists us to define fractional order differen¬ 
tiation and integration. That is 


i m g(s) = 


S f‘Sm-1 


g(so^dsQ • • • ds rn —2^m—i 


b J b 


(m — 


lb ^ ~ l9 ^ dtl n G N ’ 


(1.5.1) 


where g £ Li [a, 6], a, 6 £ R. If integral on right hand side in (1.5.1) converges then Euler’s Gamma 
function permits us to replace m in (1.5.1) with an arbitrary a £ M. The Riemann-Liouville 
fractional-integral of order a can be defined with the help of this natural extension as 
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Definition 1.5.2. [101, 148] The integral of fractional order a £ R+ of a function (p £ (L 1 [[c, d], K]) 
with [c, d] C K. in terms of Riemann-Liouville approach is defined as, 

X c l f(s) = -^-f (s - t) a ~ l g>(t)dt, (1.5.2) 

r (a) J o 

provided that the integral on right hand side exists. 

(1.5.2) implies Riemann version for finite lower limit of integral and it implies Liouville version 
of fractional-integral for infinite lower limit of integral. The relation 1.5.2 is known as Riemann- 
Liouville fractional integral for the case when lower limit of integral is zero. In 1823, Abel solved his 
celebrated integral equation using fractional integral operator Is. At later stage he extended the 
idea for fractional order a lies in 0 and 1, i.e. 0 < a < 1. 

Now it is the time to develop the idea for derivatives of fractional order and their important 
properties that are frequently used in the approximation theory. 

If the notation T> n is used to represent the n-th order derivative then fundamental theorem for 
calculus of integer order can be written as 

Vl a g{t) = g(t), (1.5.3) 

As an outcome of the repeated application of (1.5.3), we can write 

V n l n a g{t) = g(t). (1.5.4) 

Applying T> n on (1.5.4) and replacing n by to — n with n < m, we can write 

V n g(t) = V n V m - n l™- n g(t) = V m l™~ n g{t). (1.5.5) 

Now we are in a position to define the following definition because (1.5.5) is still justified and 
workable if n is replaced by a. 

Definition 1.5.3. [159] If a £ R, to = [a], and g £ C[0,1], then Riemann-Liouville derivative of 
fractional order a can be defined as 

V a g = V™i:~ m g. (1.5.6) 

By taking the left inverse of Riemann-Liouville integral of fractional order we get Riemann- 
Liouville derivative of fractional order but in general the converse does not hold with exception 
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of some classified class of functions. Now we discuss the following lemma and some limitations of 
Riemann-Liouville differential operator of fractional order. 

Lemma 1.5.1. [148] If a £ (n — 1, n), and m £ N, then we have the following residt 

771 fj- _ n \a—j 

iZKait) = g(t) - £ r{a _ J + 1) K~ j g(a)- (1-5.7) 

The development of the differential operator of fractional order by Riemann-Liouville played 
some meaningful role in developing the theory of differentiation and integration of arbitrary order 
but his developed theory have some flaws because it does not address properly to the real world 
physical models. It is to be observed that the Riemann-Liouville derivative of a constant is never 
zero, it is zero only for the case when the lower limit of the integral is infinite. But in several 
physical models the requirement of the limit is finite. It is to be also observed that while adopting 
Riemann-Liouville approach for differential equations of fractional order subject to initial conditions 
directs us to the initial conditions having derivative of fractional order at lower limit. These kind 
of problems can be solved successfully using some valid numerical techniques. But until now, we do 
not have any well-known physical interpretation to initial conditions prepared by using fractional 
order derivatives. So, Riemann-Liouville approach do not properly suit to handle the models based 
on real phenomenon. In 1967, Caputo presented another definition of derivative of fractional order 
to minimize the problems occurred in Riemann-Liouville approach. Later on in 1969, Caputo used 
his developed approach in viscoelasticity theory. Now it is the time to give the formal definition of 
Caputo derivative of fractional order. 


Definition 1.5.4. [159] The Caputo left derivative of fractional order a for a function (p(s ) € 
C m [b , d] can be expressed as 

= r(m^a) l (<-»)■+!-■> • m - 1 < « <m, m e N, (1.5.8) 

provided the right side is point wise defined on (6, oo). Moreover, if a is not an integer then 
m = \a] + 1, and if a is an integer then n = a. 


By the use of Equation 1.5.2 and Equation 1.5.8, it can be easily observed that 


c V a v l = F(1 + ?) y l ~ 

u o y wi i 7 


T(1 + 1 — a) 
and c VqQ = 0, if D is constant. 


ZoV = 


r ( 1 + 0 l+c 

\ y 


r(i +1 + a) 


(1.5.9) 
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Now we define a lemma which shows that the left inverse of Riemann-Liouville integral of fractional 
order gives us a Caputo derivative of fractional order. 

Lemma 1.5.2. [159] Assume g be a continues function and «el, then 

c V a a l«g(t) = g(t). (1.5.10) 

In general the converse of Lemma 1.5.3 does not hold except for a restricted class of functions. 
For this we have the following result. 

Lemma 1.5.3. [159] Assume g € AC m [a,b], a > 0, and m = [a] then 

2 a a ro(i)] = g(t) E ' DM p^ > (t ay. 

i—0 l - 


1.6 Approximation Theory 

In 1853, a well-known Russian mathematician, Chebysliev gave birth to the approximation theory, 
while investigating the problem of linkages , devices that translate the linear motion of a steam engine 
into the circular motion of a wheel, considered the following problem: 

If /i be a continuous function on [a, b] then is it possible to construct a polynomial pi(x) = 
£" =0 c iA °f degree at most n, such that the error (max a < x <f> \fi(x) — pi(x)\) is minimized? 

The construction and choice of pi(x) is very important to get the accuracy in numerical results. 
Currently, many scientists are using the orthogonal polynomials to develop the approximation theory 
and efficient numerical schemes for integer and fractional order differential equations. The names of 
few of among are, Hermite polynomials, Laguerre polynomials, Jacobi polynomials, and Legendre 
polynomials. In this thesis we are using the generalized two parametric shifted Jacobi polynomials 
(JPs), which is defined as 

1.6.1 Shifted Jacobi Polynomials [75, 163] 

A well reputed scientist Carl Gustav Jacob Jacobi (1804-1851) defined the Jacobi polynomials which 
are actually the solutions of homogeneous differential equations of the second order of the form 


(1 — z 2 )u'[(x ) + (a — /3 — (a + /3 + 2 )x)u' 1 (x) + n(n + /3 + a + \)ui{x) = 0. 


( 1 . 6 . 1 ) 
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If we choose (3,a> —1, and nSN then a polynomial of order n is solution of (1.6.1). It is defined 
as 

q(M( x ) = r (/3 + n + jO y' ( n r ( n + m + P + a + 1) ( x- 1 \ m . 

Qn U n\Y(fi + a + n +1) (mj r(/3 + n+l) V 2 J ' 

It is to be noted that JPs are orthogonal on the interval [—1,1] with respect to the weight function 
(1 — x)^(l +x) a . Under the transformation x = A — 1, we can obtained the two-parametric shifted 
JPs, defined as 




(-iy- k T{i + a + l)T(i + k + P + a + l) k . 


k =0 


T(k + a + l)r(z + /3 + a + l)(z — k)\k\rj k 


z\ « = 0,1,2,3, .... 


In the light of (1.6.3), the following relations are obtained [75] 


o (M (o) - ( iv ; + Q + ^ 
i; r(* + l)r(a + l) 

Qf“\ A)= ^ + /3 + 1 ) 


r(z + i)r(/3 + i)’ 


and 


max ze y 0 A \\Q { l]f > {z)\ < 


where D(i, Q = r( ^r(c+i) and C = max{a, 7 ). 

The orthogonality condition of JPs on [0, A] is as under 
r-A 


f Q { lf{z)Q { lf{z)Rf a \z)dz = W%f 8 {i , 

Jo 


where the weight function R^’ a \z) has the following form 


(1.6.3) 


(1.6.4) 


(1.6.5) 


( 1 . 6 . 6 ) 


R ( f’ a \z) = (A-zfz a , 


(1.6.7) 


and 


W. 


(P,< 


A ,j 


A /3+a+l r (j + /3+ + a + l) 

(2J + /3 + a + l)T(j + l)T(j + /3 + a + 1)' 


( 1 . 6 . 8 ) 


Due to orthogonal behaviour of JPs, any function fi(z) € C[0, A] can be expressed in terms of a 


linear combination of shifted JPs, as 


h(0 = E Cb< 2A,b Q) (u)- 

b =0 


(1.6.9) 
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Due to the practicality of the numerical problems we are interested in finding the truncated sum, 
therefore we can write (1.6.9) in terms of truncated sum, as 

M 

h{ z )-^Z c bQ { A,b\ z )- ( 1 . 6 . 10 ) 

6=0 

The obtained numerical results approaches to the exact results as TV —► oo. Therefore with the help 
of (1.6.6), (1.6.7), and (1.6.8), the coefficients c(b) can be easily determined. In vector notation 
(1.6.10) has the following form 

( 1 . 6 . 11 ) 

Where M = m +1 is actually a size of vectors and we use it as a scale level to develop the examples 
of the proposed numerical schemes. The coefficient vector is B^ and finally SIrepresents the 
function vector. The reader can consult [75, 163] for the detail study of JPs. 


1.6.2 Error Analysis 


Practically, It is very important to take the finite terms while dealing to the series representation 
of a sufficiently smooth function /i G C[a , b ]. Assume Y\ M {z) denotes the space of M terms shifted 
JPs and the best approximation of fi(z) is f(M)( z ) in the given space Then it leads us to 

the following relation 

ll/i(2) - /(m) ( 2 ) 11 2 < ||/i( 2 ) - E (m) {z )|| 2 , (1.6.12) 


where E( M \ ( 2 ) is an arbitrary polynomial having degree < M. Then following the same proce¬ 
dure as mentioned in [76], the following relation can easily be obtained representing the error of 
approximation 

A', 

(1.6.13) 


where 


Whiz) - E {M ){z)\\ 2 < jfM+i-. 
qm+i 


K\ = — max 


r/i ( 2 ) | - 


(1.6.14) 


4z7[a?b\'dz M + 1 ‘ 

But in our case K 1 will be calculated on the domain [0, A] instead of [a, b\. Hence (1.6.14) becomes 

1 d M+1 

Kl = i.”o a 5i l a?m / ' w| - 
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1.6.3 Spectral Accuracy 

The decay of the expansion coefficients and spectral accuracy can be investigated with the help of 
the following lemma [31, 88]. 


Lemma 1.6.1. If IljvM taken as the space of orthogonal-shifted JPs of N terms then for any 
fi(z) = Ejlo h iQ { ff\ z ) ^ rijvMj we have the following rusult; 

TS M oo 

I hj I- II /i (M) II and II h(z) - ^2 d iQ { f,i a) ( z ) f = (1.6.15) 

'W' 1=0 l=M 

where pi = 1(1 + (3 + a + 1) and di = fi(z)Q^ a \z)Ri a ’^\z)dz. Moreover, K is dealt as a 

constant and M can be chosen in a way such that /i (2M) e IU*)- Also we have the equality 

/ \ M 


f 1(M) 


K 


(/ 3 .«) 




Lfl ( M-l)( Z ) = 


K 


(h,a) 


( z ) 


h( z ), 


where L is the Sturm-Liouville operator and /i (0) = fi(z). 


It is to be noted that in the presence of some specific conditions JPs are well suited to expand an 
arbitrary function defined on a finite interval. Also the Jacobi spectral expansion and the convergence 
of any function defined on a finite interval depends upon the power decay of Jacobi spectral expansion 
coefficients. Moreover it is evident form Lemma 1.6.1 that for any function fi(z) £ <27°°[0,1], it is 
possible to recover the spectral decay using expansion coefficients. We mean to say that |d;| decays 
faster than any algebraic order of pi. It is worth to mention that this result does not depend on the 
BCs. For detail study on the analytical results of JPs we refer the reader to [31, 88]. 


1.7 Matrix Theory 

Our developed numerical scheme is based on the operational matrices of integrals and derivatives 
of fractional order that are used to transform the considered FOPDEs into a system of easily solv¬ 
able algebraic equations. In this section we are interested in defining these special type of matrix 
equations. 
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1.7.1 Sylvester Equation 

Let A be a square matric of order m with real entries and B be a square matrix of order n with real 
entries then Sylvester equation can be written in the following form as 

AX + XB = C. (1.7.1) 

In the light of Sylvester theorem [123, 164] for every matrix C with real entries of order m x n, (1-7.1) 
has a unique solution X of order m x n if and only if the intersection of a (A) and a(B) generates 
an empty set, where a(A) and a(B) represents the spectrums of the matrices A and B respectively. 
Numerous applications regarding to Sylvester equations are available in the literature, the names 
of few of them are, control theory, model reduction, filtering, signal processing, image restoration, 
decoupling techniques for ordinary and partial differential equations, implementation of implicit 
numerical methods for ordinary differential equations, block-diagonalization of matrices, and much 
more, see for emaple [42, 53, 61, 62, 77, 164]. Some direct methods like Stewart method [29], and 
Hessenberg-Schur method [61, 77] are used to solve (1.7.1). Under these methods the coefficient 
matrices in (1.7.1) are transformed into Schur or Hessenberg form and then the corresponding 
system of linear equations are solved applying back-substitution technique. But these methods are 
still unable to provide an explicit formula for the solutions. 


1.7.2 Vandermonde matrix 


Vandermonde matrix is known by the name of a famous researcher Alexandre-Theoophile Vander¬ 
monde, is a matrix of scalars Ki, K 2 , • • ■ ,k„ and defined as 

^ Ki k\ ■ ■ ■ 

K2 ' ' ' k 2 


H{k 1, AC 2 , - - • ,K n ) = 


(1.7.2) 


\n n ••• n n n J 

where «i < K 2 < • • • K n - These matrices have very important role in approximation problems, like, 
interpolation, moment and least square problems. Its applications can also be observed in digital 
signal processing while computing discrete fourier transform and inverse discrete fourier transform 
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[69, 78, 106]. The following Lemma is used to calculate the inverse of a Vandermonde matrix and 
will be used for our later work to calculate the inverse of operational matrices developed using JPs. 


Lemma 1.7.1. The inverse of a matrix (1-7.2) exists and is defined as 


where 


c (s,r) — ( 1 ) 


H 1 = [c (s , r) ], 


i E 

0 < ni <... < n n — k < m 


Krai ' ' ' Kra„_; 


V 


j (fi'n 

0 < n < m 
n 

For the detail study of the above Lemma, the reader can consult [107]. 


(1.7.3) 



Chapter 2 

First-Order Nonlinear Coupled 
Systems With Nonlinear Coupled 
Boundary Conditions 


The aim of this chapter is to investigate the existence of solutions of nonlinear first-order ordinary 
coupled systems with nonlinear CBCs using LUSs approach in the presence of the existence of a 
pair of coupled lower solution ( 011 , 0 : 2 ) and upper solution (/?i,/ 3 2 ) such that ( 011 , 0 : 2 ) ^ {PuPi) on 
[0,1] x [0,1]. We also discuss the existence of solutions when lower and upper solutions are in reverse 
order, that is, (oi,o 2 ) >r (Pi,P 2 ) on [0,1] x [0,1]. 

The applications of nonlinear BVPs can be observed in almost every held of science, including 
population dynamics, mechanics, biotechnology, optimal control, physics, ecology, harvesting, and 
many more [139, 191, 192], 

While dealing with nonlinear ordinary differential systems(ODSs) mostly authors only focus their 
attention on the differential systems with uncoupled BCs [18, 44, 91]. However, very rare research 
work is available where the differential systems are coupled not only in the differential systems but 
also through the BCs [16, 180]. We deal to the latter case. 

In this chapter, we deal to the existence of solutions of the first-order nonlinear coupled systems 
of the type 


u[(t) = t £ [0,1], 

u' 2 {t) = t £ [0, 1], 


( 2 . 0 . 1 ) 


44 
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subject to nonlinear CBC 

^(« 1 (0),U2(0),« 1 (1),U 2 (1)) = (0,0), (2.0.2) 

where the nonlinear functions fi, f 2 '■ [0,1] xM-}R and <f> : R 4 —> M 2 are continuous. 

The very interesting and productive aspect of this chapter is the generalization of the classical 
existence results that had been discussed in [120, 126, 140]. We mean to say if (j)(x,y,z, w) = 
(x — z,y — w), then (2.0.2) implies the periodic BCs. Also if <f>(x,y,z,w) = (x + z,y + w), then 
(2.0.2) implies the anti-periodic BCs. Definitely, in order to obtain a solution satisfying some initial 
or BCs and lying between a lower solution and an upper solution, we need additional conditions. 
For example, in the periodic case it suffices that 

c*i(0) < «i(l), a 2 (0) < a 2 (l), 

ft(0) > 0i(l), 02(0) > 0 2 (1), (2.0.3) 

and in the anti-periodic case it suffices that 

«i(0)<— 0r(l), a 2 (0) <-02(1), 

0i(0) > -oi(l), 02(0) > — 0 : 2 ( 1 ). (2.0.4) 

We extend the idea of CLUSs presented in [66] for a nonlinear first-order coupled systems in 
Section 2.2 that allows us to get a solution in the sector [oi,0i] x [[ 02 , 02 ] or [0i,ai] x [0 2 ,a 2 ] 
and unifies the treatment of the classical results (2.0.3) and (2.0.4) applying some nronotonicity 
assumptions on the arguments of the boundary function <f>. 

This chapter is organized as: In Section 2.1, the idea of LUSs is recalled. In the same section, the 
result presented in the Preliminary chapter in the form of the Lemma 1.2.1 is extended to calculate 
the inverse of the continuous operator L. In Section 2.2, the idea of CLUSs is discussed that unify 
the treatment of many linear and nonlinear BVPs. In Section 2.3, the existence of at-least one 
solution is ensured applying coupled LUSs approach. In Section 2.4, the validity of the theoretical 
results developed in Section 2.3 is checked with the help of an example. In Section 2.5, the existence 
of at-least one solution is ensured when lower and upper solutions are in reverse order. Moreover, 
in the same section, the validity of the theoretical results is also checked by the help of an example. 
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2.1 Lower and upper solutions method 

In this section, we recall the concept of lower and upper solutions for the nonlinear coupled systems 
of the type (2.0.1) and a very useful result in the form of a lemma that is very necessary to prove 
the main results of the Section 2.3 and the Section 2.5.2. 

Definition 2.1.1. We say that a couple of function ( 0 : 1 , 02 ) G C 1 [0,1] x C ,1 [0,1] is a lower solution 
of (2.0.1) if 

< /1 (t,o 2 (f)), te[0,l], 

< / 2 (t,on(t)), t G [0,1]. (2-1.1) 

In the same way, an upper solution is a function (/?i,/3 2 ) G C ,1 [0,1] x C 1 [0,1] that satisfies 
the reversed inequalities in (2.1.1). In what follows we shall assume that (oi,o 2 ) ^ (/?i,/ 3 2 ), if 
Oi (t) < and o 2 (t) < /3 2 (t), for all t £ [0,1] or (oi,o 2 ) y (/?i,/3 2 ), if Oi (t) > f3\ (t) and 

o 2 (t) > /3 2 (t), for all t G [0,1]. 

For u,v G C 1 [0,1], we define the set 

[u,i>] = {rc G C 1 [0,1] : u(t ) < w(t) < v(t),t G [0,1]} . 
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_ 7 + A6/ Q 1 e A(s 1 '>y(s)ds - by(l) 
a + be~ x 

5 + Xd e x ( s ~^z(s) ds — dz(l) 
c + de~ x 


Proof. Choose 


y(t) = ui(t) - ui(0) + A [ Mi(s)ds, 

Jo 

z(t) = u 2 {t) - U 2 (0) + A f U 2 (s) ds, 

Jo 

7 = oui(0) + bu i (1), 

and 

<5 = cu 2 ( 0) + du 2 (l)- 
In the light of (2.1.4)-(2.1.7), (2.1.2) can also be written as 

[L(u!,u 2 )](t) = ( y(t),z(t ), (7,5)). 

Differentiating (2.1.4) w.r.t. t, we have 

y\t) = u[(t) + Xui (t). 

Multiplying (2.1.9) with integrating factor e xt , we have 

e xt y'(t) = (m(t)e M )', 

then after integrating and taking the limits of integration from 0 to t, (2.1.10) becomes 


u 1 (t)=u 1 (0)e xt + y(t)~ A f e A(s t) y{s)ds , 

Jo 


ui(0) can easily be determined with the help of (2.1.6) as 


7 = (a + be A ) Ui(0) + by(l) — bX f e x ^ s ^y(s)ds, 

Jo 


(2.1.4) 


(2.1.5) 

( 2 . 1 . 6 ) 


(2.1.7) 


( 2 . 1 . 8 ) 


(2.1.9) 


( 2 . 1 . 10 ) 


( 2 . 1 . 11 ) 


, 7 + 6A C 1 e A ^ s ^y(s)ds — by(l) , 

«i 0 = -- ^^ -— > a + be~ x ^ 0, 

a + be A 


for simplicity of notation, let 


7 + 6A L 1 e A ^ s 1 ^y(s)ds — by(l) , 

A=J - ^ a + be~ x ^ 0. 


; + be~ x 


Using (2.1.13) in (2.1.11), we have 


ui(t) = Ae xt + y(t) — X f e x( - s ^ y(s)ds. 

Jo 


( 2 . 1 . 12 ) 


(2.1.13) 


(2.1.14) 
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Similarly along the same lines, it can easily be shown that 


u 2 (t) = Be xt + z(t) — A f e x ( s ^z(s)ds , 

Jo 

(2.1.15) 

with 

S + dA [* e x ( s ~ 1 )z(s)ds ~ dz(l) . 

B= -^, c + de-V 0; 

c + de A 

(2.1.16) 

(2.1.8) can also be written as 


[L - 1 (y(t),z(t),( 7,5))] = (u 1 (t),u 2 (t)). 

(2.1.17) 

Hence, (2.1.13)-(2.1.16) prove the result. 

□ 


2.2 Coupled lower and upper solutions 

In this section we define the concept of CLUSs that unifies the treatment of many linear and non 
linear first order BVPs in the presence of some nronotonicity assumptions on the arguments of the 
boundary function f. Now we look towards the definition of CLUSs. 

Definition 2.2.1. We say that a couple of functions ((cti, 02 ), (/?i, P 2 )) G C ,1 [0,1] x C' 1 [0,1] are 
coupled lower and upper solutions for the problem (2.0.1) and (2.0.2) if ( 01 , 02 ) is a lower solution 
and (/?r, /S 2 ) an upper solution for the system (2.0.1) with ( 01 , 02 ) p, (Pi,^ 2 ), such that 

</»(oi(0),a 2 (0),oi(l),o 2 (l)) ^ (0,0) ^ 0(/3 i(O),/3 2 (O),/5i(1),/5 2 (1)) (2.2.1) 

and 

</»(o 1 (0),O2(0),/3i(1),/3 2 (1))^ (0,0)^^(/3i(0),/3 2 (0),oi(1),o 2 (1)). (2.2.2) 

2.3 Existence of at least one solution 


In this section we develop a result which ensures the existence of at least one solution of the problem 
(2.0.1)-(2.0.2). The statement and proof of the result is as under: 

Theorem 2.3.1. Assume that a couple of functions ( 01 , 02 ) and (pi,fi 2 ) are coupled lower and 
upper solutions respectively for the coupled BVPs of the type (2.0.1)-(2.0.2). In addition, suppose 
that the functions 


^(ai,a 2 )( ; U2/) := 0(ai(O),a 2 (O),£ ,y ), 
^(/h.ftoOu^) ; = (A(0), #2(0), £,j/), 


are monotone on [ai(l),/3i(l)] x [ 02 ( 1 ), fi 2 (I)\, then the coupled BVP (2.0.1)-(2.0.2) has at least one 
solution (ui,u 2 ) G [ai,/3i] x [o2,/3 2 ]. 
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Proof. Let A > 0, and consider the modified coupled BVP 

u^t) + Xui(t) = F*(t,ui(t),u 2 {t)), fe [0,1], 
u 2 (t) + Xu 2 (t) = Ff(t,u 1 {t),u 2 (t)), t£ [0,1], 

4 >* (mi(0),m 2 (0),mi(1),m 2 (1)) = («i(0),u 2 (0)), (2.3.1) 


with 


Ff(t,u 1 (t),u 2 (t)) = < 


fi {t,p2(t)) + A/?i (t) 

fl (t,U 2 (t)) + \p! (t) 

fi (t,a 2 (t)) + A/?i(t) 

fl (t)) + \Ul(t) 

fi (t,u 2 (t)) + Ami (t) 


fi (t,a 2 {t)) + Ami (t) 
fi (t,p2(t)) + Aai (t) 
fi (t,u 2 (t)) + Aai(i) 
fi (t, a 2 (t)) + A«i(t) 


if (mi( t), m 2 (f)) ^ (Pi{t),P 2 {t)), 
if a 2 (t) < u 2 (t ) < fi 2 (t),ui(t) > /3i(t), 
if u 2 (t) < a 2 (t),ui(t) > Pi(t), 

if M 2 (i) > /3 2 (t),ai(t) < Mi(t) < 

if a 2 (t) < u 2 (t) < p 2 (t), 

oil(t) < Mi(t) < /3i(t), 

if u 2 (t) < a 2 (t),ai(t) < Mi(t) < Pi (t), 
if u 2 (t) > p 2 (t),ui(t) < ai(t), 
if o; 2 (t) < M 2 (t) < P 2 (t),ui(t ) < ai(t), 
if m 2 (t) < a 2 (t),ui(t) < «i(t), 


and 


F 2 (t,ui(t),u 2 (t)) 


fi (L Pi(t)) + Xp 2 {t) 
fi (t,ui(t)) + Xp 2 (t) 
/ 2 (t,ai(i)) + Xp 2 (t) 
h + Xu 2 {t) 

fi {t,ui(t.)) + Xu 2 (t) 


h (t, Q-l(t)) + Au 2 (t) 
/ 2 (t,Pi(t)) + Xa 2 {t) 
fi (t,ui(t)) + Xa 2 (t) 
t / 2 (L 04 (f)) + A a 2 (t) 


if (Mi(t),M 2 (t)) ^ (Pi(t),P 2 (t)), 

if on (t) < Mi(t) < Pi(t),u 2 (t ) > /? 2 (f), 

if Mi(t) < a 1 (t),M 2 (t) > /3 2 (t), 

if Mi (t) > pi(t),a 2 (t) < u 2 (t) < P 2 (t), 
if «i(f) < «i(t) < Pi{t), 
a 2 {t) < u 2 (t) < P 2 {t), 
if Mi (f ) < ai(t),a 2 (t) < u 2 (t) < p 2 (t), 
if Mi (t) > Pi(t),u 2 (t) < a 2 (t), 
if oii(t) < Mi(t) < Pi(t),u 2 (t) < a 2 (t), 
if Mi (t) < «i (i),u 2 (f) < a 2 (t), 


4>* (x, y, z, w) = p( 0 , (x, y )) - <f> (x, y, z, w ), 


and 

f (Pi(t),P 2 (t)) if (x,y) js (P!(t),p 2 (t)), 
p(t,(x,y)) = < (x,y) if (ai(t),a 2 (t)) ^ (x,y) ^ (Pi{t), p 2 {t)), 

[(ai(f),a 2 (f)) if (x,y) ^ («i (t),a 2 (t)). 

Note that if (mi,m 2 ) G [ai,/3i] x [a 2 ,/3 2 ] is a solution of the coupled BVP (2.5.3), then (mi,m 2 ) is a 
solution of the coupled BVP (2.0.1)-(2.0.2). 

For the sake of simplicity we divide the proof into three steps: 

We define the mappings 

L,N : C 1 [0,1] x C 1 [0,1] ->• C] [0,1] x C] [0, l]xR 2 , 
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by 


[L(m, u 2 )](t) = ^«i(i) — ui( 0 ) + A J ui(s) ds,u 2 (t) — u 2 (0) + X J u 2 (s)ds, 

(ui(0),tt 2 (0)) 


and 


t ^ 

[N(u l7 u 2 )](t) = (y F*(s,u 1 (s),u 2 {s))ds, J F 2 * (s,ui(s),u 2 (s))ds, 


<£* (ui(0),u 2 (0),tti(l),u 2 (l)) 


Since .F*(s,Ui(s),u 2 (s)) and _F 2 (s, Ui(s), u 2 (s)) are bounded on [0,1] x K 2 and integral is a con¬ 
tinuous function on [0,1]. Further <j>* being constant function is continuous. Therefore [ N(u\,u 2 )] 
is continuous on [0,1]. Further, the class {N(ui,u 2 ) : u\,u 2 G C 1 [0,1]} is uniformly bounded and 
equicontinuous. Therefore in view of Arzela-Ascoli theorem { N{u\,u 2 ) : u\,u 2 G C 1 [0,1]} is rela¬ 
tively compact. Consequently N is a compact map. Moreover from Lemma 2.1.1 with a = 1, b = 0, 
c = 1 , and d = 0 , L -1 exists and is continuous. 

On the other hand, solving (2.5.3) is equivalent to finding a fixed point of 

L-'N : C 1 [0,1] x C 1 [0,1] -► C^O,!] x C^O.l]. 

Now the Schauder fixed point theorem guarantees the existence of at least a fixed point since L~ l N 
is continuous and compact. 

Step 2: It remains to show that (ui,u 2 ) G [ai,/3i] x [a 2 ,/? 2 ]. 

We claim that (ui,u 2 ) A (/?i,/3 2 ). If (wi,it 2 ) ^ (Pi,P 2 ), then iq ^ and/or u 2 ^ fi 2 . If rq ^ /3i, 
then there exist some ?’o G [0,1], such that tq — attains a positive maximum at ?’o G [0,1]. We 
shall consider three cases. 

Case 1. r o G (0,1]. Then there exists £ G (0,r 0 ), such that 0 < ui(t) — /?i (t) < rq(ro) — /?i(ro), for 
all t G [£,ro). This yields a contradiction, since 


Pi{r 0 ) - /?i(0 < wi(?’o) - ui(0 

J ^/i(s,/32(s))-A(ui(s)-/3i(s))^ds 


< 


fi(s,/3 2 {s))ds= / Px(s)ds = /3i(r 0 ) - /?i(0- 


Case 2. ro = 0 and <j)p is monotone nonincreasing. Then rq(0) — /3i(0) > 0 or it 2 (0) — /? 2 (0) > 0, 
and in view of ( 2 . 2 . 1 ), we have 


(u(0),u(0)) = (j)* (ui(0),u 2 (0),u 1 (l),u 2 (l)) 

= (/3i(0),/3 2 (0)) - </>(A(0),/3 2 (0),ui(l),u 2 (l)) 

— (&(()), ft(0)) - </>(/?i(0),/3 2 (0),/? 1 (1),/3 2 (1)) 

— (ft (0), &(<))), 


(2.3.2) 


a contradiction. 
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Case 3. Similarly <pfj is monotone nondecreasing. We shall change the inequality (2.3.2) by (ui(0), u 2 (0)) ■< 
(/3i(0), /S 2 (0 ))—<f (/3i(0), /? 2 (0), ai(l), 012 ( 1 )) and again we get a contradiction. Consequently, (u\,u 2 ) di 
(/3i,/ 3 2 ), for all t G [0,1]. Similarly, we can show that (ui,u 2 ) t (ai,a 2 ), for all t G [0,1]. 

Step 3: Now, it remains to show that (ui,« 2 ) satisfies the boundary condition (2.0.2). 

For this, we claim that 

M0),a 2 (0)) < K(0),« 2 (0)) -<f>( Ul (0),« 2 (0), U 1 (1),W 2 (1))) 1 (/3i(0), /3 2 (0)). (2.3.3) 

If (m(0),it 2 (0)) - <^(ui(0),u 2 (0),ui(l),w 2 (l)) ^ (/3i(0),/3 2 (0)), then 
(ui(0),m 2 (0)) = ( j >* (mi(0),m 2 (0),Ui(1),u 2 (1)) 

= p(0, (ui(0),u 2 (0))) - 0(«i(O),u 2 (O),ui(l),u 2 (l)) 

= (/3 1 (0),/3 2 (0)). 

If <bp{x. y) is monotone nonincreasing, then we have 

(«i( 0 ),u 2 ( 0 )) - <j> (ui( 0 ), u 2 ( 0 ), Ui(l), w 2 (l)) 

= (/?i(0), &(0)) - </> (/3i(0), /? 2 (0), Ul (l), u 2 (1)) 

= (A( 0 ),/? 2 ( 0 )) - (j>p (ui(l),u 2 (l)) 

d (/3i(0),/3 2 (0))-^(/3 1 (l),/3 2 (l)) 

= (/3i(0),/3 2 (0)) - ^>(/3 1 (0),/3 2 (0),/3 1 (l),/3 2 (l)) 
d (/?i(0),/3 2 (0)), (2.3.4) 

a contradiction. Similarly if <j)p(x, y) is monotone nondecreasing, then we get the same contradiction. 
Consequently, (2.3.3) holds. Hence the system of BVPs (2.0.1)-(2.0.2) has a solution (u% r u 2 ) G 
[<* 1 , Pi} x [a 2 )^ 2 ]- n 

2.4 Application of the theoretical results 


In this section, we check the applicability of the developed result 2.3.1. 

Example 2.4.1. Let 


fi{t,u 2 (t)) = -2u 2 (t) + 7 sin(wt), t G [0,1], 
f 2 {t,ui(t)) = -2u\(t) + 7 COs(wt), t G [0,1], 
y, z, w) = (a; 3 - y 3 , z 3 - w 3 ), 

where 7 and ui are positive integers. Let a\(t) = —2"/, a 2 (t ) = —7 and fd\(t) = 27 , f3 2 (t) = 7 . 

In the light of the methodology adopted in Example 1 . 2.5 of the Preliminary chapter, it is evident 
that a couple of functions (ai,a 2 ), (/3i,f3 2 ) are lower and upper solutions respectively of the system 
(2.0.1). On the same fashion as mentioned in Example 1 . 2 . 5 , it can be shown that a couple of 
functions («i,a 2 ) and (/3i,/3 2 ) satisfy the system of inequalities (2.2.1)-(2.2.2). Hence by Theorem 
2 . 3 . 1 , the coupled BVPs of the type (2.0.1)-(2.0.2) has at least one solution (u\,u 2 ) G [ai,/3i] x 

[a 2 ,/3 2 ]- 
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2.5 Lower and upper solutions in reverse order 


In this section, we discuss the existence of solutions of the coupled BVPs of the type (2.0.1)-(2.0.2) 
when a pair of lower and upper solutions are in reverse order, i.e., ( 011 , 0 : 2 ) cl (Pi,P 2 )- Let us see 
the definition of coupled lower and upper solutions for the case when (oi(t), 02 (f)) >r (Pi(t), P 2 (t)), 
for all t G [0,1]. 


2.5.1 Coupled lower and upper solutions 

Definition 2.5.1. We say that a couple of functions ((oi, 02 ), (Pi, P 2 )) G C , 1 [0,1] x (C 1 [0,1] are 
coupled lower and upper solutions for the coupled BVP (2.0.1)-(2.0.2), if ( 01 , 0 : 2 ) is a lower solution 
and (Pi,p 2 ) an upper solution for the system ( 2 . 0 . 1 ) with ( 01 , 02 ) h (Pi,Pi), such that 

</»(o 1 (0),o 2 (0),o 1 (l),o 2 (l)) ^ (0,0) ^ 0(/3 1 (O),/3 2 (O),/5 i(1),/3 2 (1)) (2.5.1) 

and 

^(/3 1 (0),/3 2 (0),o 1 (l),a 2 (l))^ (0,0)^^(o 1 (0),a 2 (0),/3 1 (l),^ 2 (l)). (2.5.2) 


2.5.2 Existence of at least one solution 


Theorem 2.5.1. Assume that a couple of functions (oi, o 2 ) and (pi,p 2 ) are coupled lower and upper 
solutions such that ( 01 , 02 ) cl (Pi, @ 2 ) f or the coupled BVP (2.0.1)-(2.0.2). In addition, suppose that 
the functions 

<t>(a u ai)(x,y) : = 0( x ,y,oc 1 (l),a 2 (l)), 

<l>(Pi,fo)( x >y) '■= <t>(x,y,Pi(l),P 2 (l)), 

are monotone in [/?i(0), oi(0)] x [/^(O), o 2 (0)], then the system (2.0.1)-(2.0.2) has at least one solution 
(ui,U 2 ) G [/3l,Oi] x \p 2 , o 2 ]. 

Proof. Let A > 0, and consider the modihed coupled BVP 

■u^t) - Xui(t) = Ff(t,ui(t),u 2 (t)), t G [0,1], 
u 2 (t) - Xu 2 (t) = Ff(t,ui(t),u 2 (t)), t G [0,1], 

(f* (tii(0),u 2 (0),u 1 (l),u 2 (l)) = (u 1 (l),u 2 (l)), (2.5.3) 


with 


Ff(t,ui(t),u 2 (t)) = < 


fi(t,p 2 (t))~Xpi(t) 
fi (t,u 2 (t)) - XPi(t) 
fi (t,a 2 (t)) - XPi(t) 
fi (t,P 2 (t)) - Xui(t) 
fi (t,u 2 (t)) - Xui(t) 


fi (t,a 2 (t)) - Xui(t) 
fi (t,P 2 (t)) - Xai(t) 
fi (t,u 2 (t)) - Aoi (t) 
fi (t,a 2 (t)) - Aoi(£) 


if (m(t),u 2 (t)) ^ (Pi(t), P 2 (t)), 
if o 2 (£) < u 2 (t) < p 2 (t),ui(t) > Pi(t), 
if u 2 (t) < a 2 (t),ui(t) > Pi(t), 
if u 2 (t) > P 2 (t),a i(£) < u\(t) < Pi(t), 
if a 2 (t) < u 2 (t) < p 2 (t), 

01 (t) < ui(t) < Pi(t), 
if u 2 (t) < a 2 (t),ai(t) < ui(t) < pi(t), 
if u 2 (t) > p 2 (t),ui(t) < oi(£), 
if a 2 (t) < u 2 (t) < P 2 (t),ui(t) < Oi(£), 
if u 2 (t) < a 2 (t),ui(t) < oi(£), 
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and 


F^(t,u 1 (t),u 2 {t)) = < 


h (t,Pi{t)) - A fc(t) 
h (t,U!(t)) - A p 2 (t) 
h (t,ai(t)) - A fi 2 (t) 
h (t,Pi(t)) - A u 2 {t) 
h - A u 2 (t) 


fi (t,a i(t)) - A u 2 (t) 
h - A a 2 (t) 

fi - A a 2 (t) 

fi (*, o;i(t)) - A a 2 (t) 


if (wi(i),u 2 (0) ^ {Pi(t),Pi(t)), 
if qi( t) < m(t) < fh(t),u 2 (i) > /3 2 {t), 
if m(t) < ai(t),u 2 (t) > p 2 (t), 
if m(t) > pi(t),a 2 (t) < u 2 (t) < /3 2 (t), 
if aiW < «i(t) < Pi(t), 

a 2 (t) < u 2 (t ) < fi 2 (t), 
if ui(t) < a\(t),a 2 (t) < u 2 {t) < P 2 {t), 
if m(t) > pi(t), u 2 {t) < a 2 (t), 
if oii(t) < ui(t) < f3i(t),u 2 (t ) < a 2 (t), 
if ui(t) < cei(t),u 2 (t) < a 2 (t), 


<f>* (x, y, z, w) = p{ 1, {x, y)) - <$> (x, y, z, w ), 


and 

( (Pi(t),p 2 (t)) if {x,y) ^ (Pi(t),P 2 (t)), 
p(t,(x,y)) = S (x,y) if («i (t),a 2 (t)) < (x,y) < (Pi(t), P 2 (t)), 

[(ai(f),a 2 (i)) if {x,y) ^ («i(t),a 2 (t)). 

Now the rest of the proof is analogous to the proof of the Theorem 2.3.1 using Lemma 2.1.1 with 
a = 0, 6 = 1, c = 0, and d= 1. We left it as an exercise. □ 

2.5.3 Application of the theoretical results 

Example 2.5.1. Let 


fi{t,u 2 (t)) = 4u 2 (t) + 7 sin(wf), t & [0,1], 

/ 2 (t,ui(t)) = 4uf(t)+7C0s(wt), t€ [0,1], 

<t>{x, y, z, w) = (yw - xz,y + xz), 

where 7 and u> are positive integers. Choose 07 (f) = 37 , a 2 (t) = 27 , P\{t) = — 37 , p 2 (t) = — 27 . 

In the light of the methodology adopted in Example 1.2.8 of the Preliminary chapter, it is evident 
that a couple of functions ( 07 , a 2 ), {Pi, p 2 ) are lower and upper solutions respectively in the reverse 
order of the system (2.0.1). On the same fashion as mentioned in Example 1.2.8, it can be shown 
that a couple of functions ( 07 , a 2 ) and {Pi,P 2 ) satisfy the system of inequalities (2.5.1)-(2.5.2). 
Hence by Theorem 2.5.1, the system of BVPs (2.0.1)-(2.0.2) has at least one solution (u\,u 2 ) € 
\Pi, a-i] x [ f} 2 , a 2 \. 



Chapter 3 

Second-Order Nonlinear Coupled 
Systems With Nonlinear Coupled 
Boundary Conditions 


In this chapter, we extend the ideas of Chapter 2 for a second order nonlinear coupled BVPs 
corresponding to nonlinear CBCs in the presence of the existence of well ordered LUSs, such that 
Oi,a 2 ) ^ on [ 0 , 1 ] x [ 0 , 1 ]. 

The study of differential systems with CBCs is a very important area of analysis due to not only 
the theoretical aspects but also the rich applications in the existing literature; including reaction- 
diffusion equations, Sturm-Liouville problems, mathematical biology, chemical systems and engi¬ 
neering, see for example [1, 2, 5, 37, 117, 178, 188]. In 1982, Leung discussed the applications 
of coupled BVPs in mathematical biology by studying the following reaction-diffusion system for 
prey-predator interaction [117]: 

dv,i(t,x) _ + m (a + fi(ui,u 2 )) , t > 0, x G 13 C R", 

vt ( 3 . 0 . 1 ) 

^dt X ^ = a ' 2 ^ u ’ 1 + (r + f'2 (ui, U2 )), t> 0 , lellCl", 

subject to CBCs 


Ifrf = “ p( m i) ~ q ( U2 ^ = ° ° n (3.0.2) 

where fi , / 2 : R 2 — > ffi. have Holder continuous partial derivatives up to second order in compact 
sets, p is a unit outward normal at 815 and p and q have Holder continuous first derivatives in 
compact subsets of [0, oo). The functions U\{t,x) and u 2 (t,:r) respectively represent the density 
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of prey and predator at t > 0 and at position x = (xi, X 2 , ..., x n ). In 1978, Aronson studied the 
biochemical system using same type of CBCs [2]. In 2009, Cardanobile and Mugnolo studied the 
following parabolic system with CBCs [37]: 


x) = Cu\{t, x), t > 0, x €13 & R", 
fi\dl3 € y, jr— € 

OU2 


(3.0.3) 


where y 1 - is closed subspace of L 2 (dfl;W), the unknown function u-\ takes values in a separable 
Hilbert space 11’ and C is an elliptic operator with operator-valued symbol. 

This chapter is organized as: In Section 3.1, the method of LUSs for the second order nonlinear 
coupled systems is briefly discussed. In the same section, the inverse of the continuous operator L is 
determined by extending the result of the Chapter 2 discussed in the Section 2.1. In Section 3.1.2, 
the concept of CLUSs is discussed that unify the treatment of the classical existence results of many 
linear and nonlinear BVPs, like Dirichlet and Neumann by imposing some monotone assumptions 
on the arguments of the boundary functions <p and ip. In Section 3.1.3, the existence of at-least one 
solution is ensured for the nonlinear coupled BVPs applying coupled LUSs approach. In Section 
3.1.4, the validity of the theoretical results is checked by taking some examples. In Section 3.2.1, the 
method of LUSs is recalled for the nonlinear second-order coupled BVPs corresponding to nonlinear 
coupled BCs. In Section 3.2.2, the idea of CLUSs is discussed that generalizes the treatment of the 
classical existence results for very important BVPs, like periodic and anti-periodic by imposing some 
monotonicity assumptions on the arguments of the boundary functions /x and v. In Section 3.2.3, 
the existence of at-least one solution of the second-order coupled BVPs is established using coupled 
LUSs approach. In Section 3.2.4, some illustrative examples are considered to ensure the validity of 
the theoretical results. 

3.1 Second-order coupled system with nonlinear coupled bound¬ 
ary conditions (I) 


In this section, we investigate the existence of solutions of the following second order ordinary 


nonlinear coupled BVPs of the type 


«i (t) = fi(t,u 2 (t)), te[o,i], 

U 2 W = Mt,ui{t)), t € [0,1], 


(3.1.1) 
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with nonlinear CBCs 

( j > (mi(0),m 2 (0),mi(1),w 2 (1),w , i(0),U2(0)) = (0,0), 

^ ^wi(0),u 2 (0), ui(l),u 2 (l),u 1 (l), w 2 (l)) = (0,0), 


(3.1.2) 


where /i,/ 2 : [0,1] x K. — > R. and (j), ip : R 6 —► ]R 2 are continuous functions. The existence results 
for the second order BVPs has been established by many authors in the presence of certain types 
of BCs, see for example [16, 86, 95, 96, 99, 181, 186, 193]. In contrast to previous studies our BCs 
are very much generalized because they unify the treatment of many important BVPs like Dirichlet, 
Neumann, periodic and anti-periodic. The detail note on the general and broad behaviour of our 
chosen BCs 3.1.2 and 3.2.2 has been comprehensively discussed in the introductory chapter of this 
thesis. 

The study of differential systems becomes more interesting and some times complicated when cou¬ 
pling is defined among the arguments of the BCs. Systems with CBCs can be applied to interaction 
problems, reaction diffusion phenomena and Lotka-Volterra models, see for example [6, 128, 117]. 

Zhou and Xu [193] established the existence and multiplicity of positive solutions of the following 
nonlinear coupled system of BVPs 


-u"(t) = f(t,v(t)), f 6(0,1), 

~v"{t) = g{t,u(t)), fG( 0,1), 

u(0) =u(0) =0, 

cra(?y) = u(l), av(rj) = u(l), ry € (0,1), and 0 < arj < 1, 


(3.1.3) 


with uncoupled linear-multi point BCs applying the fixed point index theory in cones. Wei and Sun 
[184] extended this study by considering the fully nonhomogeneous uncoupled multi-point BCs of 
the type 


u(0) = au(rf), tt(l) = Pu(rj), 

v(0) = av(r]), u(l) = fiv(rj), r] € (0,1), and 0 < /3 < a < 1. 


(3.1.4) 


Asif and Khan [16] extended the results for a multi point linear coupled BCs and investigated the 



57 


existence of a positive solution to the following four-point coupled BVPs by constructing a cone 

-x"(t) = f(t,x{t),y(t)), ie( 0,1), 

~y"(t) = g(t,x(t),y{t)), t G (0,1), 

(3.1.5) 

a;(0) = 0, x(l) = ay (£), 
y(0) = 0, 2 /( 1 ) = Px(r)), 

where £, r/ G (0,1), 0 < a/ 3£?7 < 1, /, </ : [0,1] x [0, oo) x [0, oo) —> [0, oo) are continuous and allowed 
to be singular at t = 0 and t = 1. On the other hand, we discuss the more generalized approach 
regarding the existence results for a second order nonlinear coupled systems by introducing the 
nonlinear coupling among the arguments of the BCs. We also claim that our considered BCs 3.1.2 
are very much generalized as compared to the BCs studied in [16, 193]. 

3.1.1 Lower and upper solutions method 

In this section, we recall the concept of LUSs and a very useful result in the form of a lemma that 
addresses the inverse of a compact map. 

Definition 3.1.1. We say that a couple of functions (on, 02 ) G C 2 [0,1] x C 2 [0,1] is a lower solution 
of (3.1.1) if 

> fi {t,a 2 (t)) , t G [0,1], 

a 2 (t) > h , t € [0,1]. (3.1.6) 


In the same way, an upper solution is a couple of functions (^ 1 ,^ 2 ) G C 2 [0,1] x C 2 [ 0,1] that 
satisfies the reversed inequalities in (3.1.6). In what follows we will write (ai,a 2 ) If? (Pufc), if 
ai(t) < j3\(t) and a 2 (t) < /3 2 (t), for all t € [0,1]. 

For U\,U 2 G C 2 \ 0,1], we define the set 

luiM = {wi G C 2 [0,1] : m(t) < Wi(t) < u 2 (t),t G [0,1]} . 


Let us construct the statement of the following lemma along with its proof . 

Lemma 3.1.1. Let L : C 2 [0,1] x C 2 [0,1] — > Cq [0,1] x Cq[0, 1]xR 2 x JR 2 be defined by 

[L{ui,u 2 )](t) = fui(t) - «i(0) - A [ wi(s) ds, u 2 (t) - u 2 (0) - A [ u 2 (s)ds, 


(aui(0) + 6ui(l), cm 2 (0) + du 2 ( 1)) , (Frn(0) + F Mi(l), Gm 2 (0) + Hu 2 ( 1)) ) ,(3.1.7) 
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where A, a , b, c, d, E, F, G and H are real constants with A > 0, such that 

{ad - be){EH - FG ) (e”^ + 0, 

here 

Co [0,1] = {w\ € C 2 [0,1] : wi(0) = 0} . 
Then L _1 exists and continuous, defined by 


[L 1 {y,z,y,S,nX)] = (Cie^* + C 2 e ^+ \ J S) y(«) ds +^ e VA(s *'?/(«) ds, 


with 


C 3 e^ + Cie-^ + l f e VX(t “ 3) 

2 Jo 


1 f y/\ (s—t) 

2 Jo 1 

z(s) ds + ^ f e v ^' ( - s ~ t> z{s) ds), 
2 Jo 


Ci = 


--r-^- r-( 23(0 +be v ^)—d(a + be v ^') 

{ad - be) (e^ - e-^J V v 7 v 7 


0 V\(l-s) 


a + be 


-vA 


o vA(s-i) 


y{s)ds 


2 7 (c + de^' /r )+&(c + (fe~' /r )^ e vA(i-«) 


y{s)ds 


— b (c + de / e 


VX(s —1) 


y(s)ds 


C 2 = 


-7-^-r- (^25 ("a + be'^'\ — d (a + be^*\ 

{ad - be) (e-'K - e^J V v 7 v 7 


0 vA(i- s ) 


a + be 


vA 


c vA(s-l) 


y{s)ds 


27 (c + de'^j + b + de'^ S ) / e v ^ < - 1 S ^y{s)ds 


— b (c + de'^j / e 


VX(s —1) 


y(s)ds 


C 3 = 


{EH - FG) (e^ - e~^j 


E + Fe~^ - h(e + Fe ~^ 


o vA(i-s) 


E + Fe 


-vA 


o vA(s-i) 


z{s)ds 


— \/a 


- 2/1 (g + He~^ +F (g + He 


- F(G + He ~^ / e VX(3 ~ 1) z{s)di 


a vA(l —s) 


z(s)ds 


(3.1.8) 
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and 


C 4 = 


Proof. Choose 


- y - V (^2C ( E + Fe- H (e + Fe^) 

(EH - FG ) (e-^ - e' 7 *) V v 7 v 7 

J e' /r * < ' 1 ~ s ^ z(s)ds + F + Fe^^j j e v W s -i) z (s)ds 
2/x (g + He^ +f(g + J e^ {1 ~ s) z(s)ds 

F (G + He^ e^ s ~ l) z(s)d^j. 


y(t) = Ul (t) - u^O) - A / Ui(s) ds, 
Jo 


z(t) = u 2 (t) — tt 2 (0) — A / M 2 (s)ds, 

J 0 


7 = aui(O) + 6 ui(l), 


6 = cu 2 ( 0 ) + du 2 ( 1 ), 


/i = £/ui( 0 ) + Fin(l), 


C = Gu 2 (0) + i?u 2 (l). 

With the help of (3.1.9)-(3.1.14), (3.1.7) can be written as 

[L 1 (u 1 ,u 2 )](t) = ( y(t ), z(t), ( 7 ,5), (/z, 0) ■ 

Differentiating (3.1.9) with respect to t. we have 

y'(t) = u'[(t) — \u\(t), A > 0. 

The general solution of (3.1.16) can be easily determined using variation of parameters 
along with integration by parts and taking limits of integration from 0 to t, we have 

Ul (t) = C 1 e VXt + C 2 e~' /Xt + l f* e' /X ^y(s)ds + \ f* e' /X ^y(s)ds. 

* Jo * Jo 

Ci and C 2 can be easily determined with the help of (3.1.11) and (3.1.12) as 

7 = (a + 6 e^) C x + (a + be~C 2 + h - (^J e' /X( 1 - s) 2 /(s)ds + e VX ^y(s)ds^j , 

5 = (c + de ^) Gi + (c + de _v/X ) C2 + i(f Q e' A{1 ~ s) y(s)ds + e' /X(s " 1) 2 /(s)ds^ . 


(3.1.9) 

(3.1.10) 

(3.1.11) 

(3.1.12) 

(3.1.13) 

(3.1.14) 

(3.1.15) 

(3.1.16) 
technique 

(3.1.17) 

(3.1.18) 
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Solving the system of equations (3.1.18), we have 

c, = 


(ad — 


g\/A g—\/A 


25 


(a + fee"' 71 ) - d(a + be 


- b 

and 

C 2 = 


e^ l -^y(s)ds + d(a + be- VX ) e^ s ~ 1) y(s)ds 
2 7 (c + de + b (c + de^' 71 ) J e^ 1 - fl) y(s)ds 
(c + de _>A ) J e' /X{s - 1) y(s)ds^j, 


(ad — 


& — g\/^ 


^25 + fee' 77 ^ — d + fee v 

J e' /X{1 - s) y(s)ds + d (a + fee' 71 ) J e VX{s - 1) y(s)ds 

- fefc + de' 71 )^ e v/X(s - 1) 2 /(s)ds^. 

Similarly on the same line, it can be easily shown that 

u 2 (t) = Cse^ + C A e-^ + \ f e^ {t ~ s> > z(s)As + \ f e' 7I(s - t) ^(s)ds, 

2 Jo 2 J 0 


with 


G, = 


E + Fe - ' 71 ) -h(e + Fe"' 71 ) 


- F 


(EH - FG) (e' 7 * - e-' 77 ') 

J e^ 1- ®)z(s)ds + F (f + Fe - ' 7 *) j e^ s ~ 1) z(s)ds 

^(0 + H^) + F(0 + H^)j\^-M^ 

(G + He-^j e^ s ~ 1) z(s)ds^, 


(3.1.19) 


(3.1.20) 


(3.1.21) 


(3.1.22) 
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and 


C 4 = 


- -4 - 2((E + Fe' /X ) -H (E + Fe 

(EH - FG) (e~^ - e^j V V J V 

J e' /X{1 - s) z(s)ds + F^E + Fe^ j e' A(s ' 1) 0(s)ds 

2 /j, (G + +F(G + He^ J e' /X{1 ~ s) z(s)ds 

\ 

z(s)ds 


F (g + He^4 J e^ 8-1 ) 


(3.1.15) can also be written as 

(ui(t), u 2 (t)) = [L _1 ( : y(t ), ^r(t), ( 7 ,5), Ob C))]- 
Hence, (3.1.17)-(3.1.24) prove the required result. 


(3.1.23) 

(3.1.24) 

□ 


3.1.2 Coupled lower and upper solutions 


In this section, we develop the dehnition of CLUSs that assists us to unify the treatment of many 
linear and non linear BVPs by imposing some monotonicity assumptions on the arguments of the 
boundary functions <p and ip. Moreover, monotonicity assumptions on <p and ip deals the classical 
existence criterion ((0.0.11)-(0.0.12)) for Dirichlet and Neumann BVPs as a particular case. Now, 
let us define the concept of CLUSs: 

Definition 3.1.2. We say that a couple of functions ( 01 , 0 : 2 ), (Pi,p 2 ) G C 2 [0,1] x C 2 [0,1] are 
coupled lower and upper solutions for the coupled BVP (3.1.1)—(3.1.2), if ( 01 , 0 : 2 ) is a lower solution 
and (/3i,/3 2 ) is an upper solution for the system (3.1.1) with ( 01 , 02 ) 7 ? (Pi, @ 2 ), such that 

^’(/3 1 (0),/? 2 ( 0 ),/3i( 1),(1),/3i( 0 ),/? 2 ( 0 )) d ( 0 , 0 ) 7 ? </»(o 1 ( 0 ),O 2 ( 0 ),oi(l),o 2 (l),o , i( 0 ),o;( 0 )), 
< p(p 1 (0),p 2 (0),a 1 (l),a 2 (l),p' 1 (0),p' 2 (0)) 7 ? (0,0) 7 ? 0(oi(O), o 2 (0), ft(l), p 2 (l), o'i(0), o^O)), 
^(/3i(0),^ 2 ( 0 ),/?i(l),/? 2 ( 1 ),/?i(l),/3 2 (1)) ^ ( 0 , 0 ) ^ ' 0 (oi(O),O 2 (O),ai(l),o 2 (l),ai(l),a 2 (l)), 
^(ai(0),o 2 (0),/3i(1),/? 2 (1),/?i(1),/? 2 C 1 )) ^ (0,0) ^ t^(/3i(0),/? 2 (0), oi(l), a 2 (l),Oi(l), o 2 (l)). 

(3.1.25) 


3.1.3 Existence of at least one solution 

Theorem 3.1.2. Assume that a couple of functions ( 01 , 02 ), (Pi,p 2 ) nre coupled lower and upper 
solutions for the coupled BVP (3.1.1)-(3.1.2). Suppose that the functions cp and ip are monotone 
nondecreasing and nonincreasing in the fifth and sixth arguments respectively. In addition, suppose 
that the functions 


4>a(x,y) ■= (p ^oi(0), o 2 (0), x, y, Oi(0), o 2 (0)^ 
<Pp(x,y) ■= (/M0),/3 2 (0),x,y,/3i(0),£ 2 (0)) 



62 


are monotone on [oi(l), /?i(l)] x [ 0 : 2 ( 1 ), ^ 2 ( 1 )] and that the functions 

f)<x(x,y) ■= i> (x,y,a i(l),o 2 (l),a , i(l), 02 ( 1 )) * 

i>p{x,y) ■= (x,y,/3 i(1),^ 2 (1),/3i(1),/3 2 ( 1 )) , 

are monotone on [oi(0),/?i(0)] x [ 02 ( 0 ), 0 2 (O )}. 

Then there exists at least one solution ( ui,u 2 ) G [ai,/?i] x [ 02 ,/3 2 ] of the coupled BVP (3.1.1)- 
(3-1.2). 

Proof. Let A > 0 and consider the following modified coupled BVP 

%(*) - Aui(t) = Ff(t,ui(t),u 2 (t)), te [0,1], 
u 2 (t) - Xu 2 (t) = F 2 (t,ui(t),u 2 (t)), te [0,1], 

<t>* (wi(0),u 2 (0),ui(l),u2(l),ui(0),u 2 (0)) = (ui(0), u 2 (0)), 

ip* (0), U 2 (0), tii(l), u 2 (1), fo (1), a 2 (!)) = (wi(l),w 2 (l)), (3.1.26) 


with 


Ff(t,ui(t),u 2 (t)) = < 


fi (t, 02(f)) - X0i(t) 
fi (t,u 2 (t)) - X0i(t) 
fi (t,a 2 (t)) - XPi(t) 
fi (t, 02(f)) - Xm(t) 
fi (t,u 2 (t)) - Xui(t) 


fi (t,a 2 (t)) - Xm(t) 
fi (f 02(f)) - Aoi(t) 
fi (t,u 2 (t)) - A 01 (t) 
fi (t,a 2 (t)) - Aoi (t) 


if (m(t),u 2 (t)) ^ (0i(t),0 2 (t)), 
if o 2 (t) < u 2 (t) < 0 2 (f),ui(t) > 0i (t), 
if u 2 (t) < a 2 (t),ui(t) > 0i(t), 
if u 2 (t) > 0 2 (t),ai(t) < ui(t) < 0i(t), 
if a 2 (t) < u 2 (t) < fa (t), 

011 (f) < ui(t) < 0i(t), 
if u 2 (t) < a 2 (t),ai(t) < ui(t) < 0i(t), 
if u 2 (t) > 0 2 (f),ui(t) < oi(f), 
if a 2 (t) < u 2 (t) < 0 2 (t),ui(t) < Oi (t), 
if u 2 (t) < a 2 (t),ui(t) < oi(t), 


and 


F 2 (t,ui(t),u 2 (t)) 


h(t,0i(f))-X0 2 (t) if (ui(t),u 2 (t)) £ (0i(t),0 2 (t)), 
h (t,ui(t)) - X0 2 (t) if oi(t) < ui(t) < 0i(f),u 2 (t) > 0 2 (t), 

h (t,oii(t)) - X0 2 (t) if ui(t) < ai(t),u 2 (t) > 0 2 (t), 

h (t,0i(t)) - Xu 2 (t) if ui(t) > 0i(t),a 2 (t) < u 2 (t) < 0 2 (t), 

f 2 (t,ui(t)) - Xu 2 (t) if oi (t) < ui(t) < 0i(t), 

o 2 (t) < u 2 (t) < 0 2 (t), 

h (t, 011 (f)) - Xu 2 (t) if m(t) < oi (t),a 2 (t) < u 2 (t) < 0 2 (t), 

f 2 (t, 0i(t)) - Xa 2 (t) if ui(t) >0i(t),u 2 (t) <a 2 (t), 

f 2 (t,m(t)) - Xa 2 (t) if oi (t) < ui(t) < 0i(t),u 2 (t) < 012 (f), 

h (t,aii(t)) - Xa 2 (t) if Ui(t) < ai(t),u 2 (t) < a 2 (t), 


<P* (j,k,l,m,n, o) 
ip* (j, k, l, m, n, o) 


P (0, (j, k) + cp (j , k, l , m, n, o)), 
P (1, (h m)+ip (j, k , l, m, n , o)) 
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and 


P(t, (x,y)) 


( (fii(t),fi 2 (t)) if (x,y) ^ (fi 1 (t),fi 2 (t)), 

< (x,y) if (ai (t),a 2 (t)) A (x,y) A (fii(t), fi 2 (t)), 

[(ai(f),a 2 (f)) if (®,y) ^ (ai(t),a 2 (t)), 


Note that if (ui,u 2 ) G [ai,/?i] x [a 2 ,/3 2 ] is a solution of (3.1.26), then (mi,m 2 ) is a solution of 
(3.1.1)-(3.1.2). 

For the sake of simplicity we divide the proof in three steps: 

Step 1: We define the mappings 


L, N : C 2 [ 0,1] x C 2 [ 0,1] -> Cg[0,1] x Cg[0,1] x ! 2 x K 2 , 


by 


[L(u 1 ,u 2 )](t) 



u 1 ( 0) — A / ui(s) ds, u 2 (t) 
Jo 


u 2 ( 0) — A / u 2 (s) ds, 
Jo 


(mi(0),u 2 (0)),(ui(1),m 2 (1)) 


and 


[iV(ui,u 2 )](t) = 


[ Ff(s,ui(s),u 2 (s))ds, f F 2 (s, ui(s),u 2 (s))ds, 
to Jo 

</>*(ui(0), u 2 ( 0), «i(l), « 2 (1), u[(0), u' 2 (0)), 
i>*{u 1 (0),u 2 (0),u 1 (l),u 2 (l),u' 1 (l),u 2 {l))\ 


Since FT(s, ui(s), u 2 (s)) and F 2 (s, ui(s),u 2 (s)) are bounded on [0, l]xR 2 and integral is a continuous 
function on [0,1]. Further <p* and ip* being constant functions are continuous. Therefore [ N(u\,u 2 )] 
is continuous on [0,1]. Further, the class { N(u\,u 2 ) : u\,u 2 G C 2 [0,1]} is uniformly bounded 
and equicontinuous. Therefore in view of Arzela-Ascoli theorem {N(ui,u 2 ) : ui,u 2 G C 2 [0,1]} is 
relatively compact. Consequently N is a compact map. Also from Lemma 3.1.1 with a = 1,6 = 0, 
c = 1, d = 0 and E = 0, F = 1, G = 0, H = 1, F _1 , exists and is continuous. 

On the other hand, solving (3.1.26) is equivalent to find a fixed point of 

L~ 1 N : C 2 [0,1] x C 2 [ 0,1] -► C 2 [0,1] x C 2 [0,1], 

Now, Schauder’s fixed point theorem guarantees the existence of at least a fixed point since L~ l N 
is continuous and compact. 

Step 2: If [u\,u 2 ) is a solution of (3.1.26), then (ui,u 2 ) G [ai,/?i] x [a 2 ,/3 2 ]. 

We claim (ui,u 2 ) A (fi\,fi 2 ). If (u\,u 2 ) ^ (/3i,f} 2 ), then either u\ ^ fii and/or u 2 ^ fi 2 . If 
ui id fii, then Mi — fii attains a positive maximum at some to G [0,1]. Clearly, ui(to) — fii(t o) > 0. 
Thus (ui — fii)' (t 0 ) = 0 and (tq — fii)" (to) < 0. But, 

(ui - fii)"(t 0 ) > i 7 )* (t 0 ,Ui(to),U 2 (t 0 )) + AUi(to) - fi (t 0 ,fi 2 (t 0 )) 

= /l (to, fil(to)) ~ A/?i(t 0 ) + Aui(f 0 ) — fl (to, fi 2 (t 0 )) 

= a (ui(t 0 ) - fii(to)) > o, 

a contradiction. Similarly, one can show that (cq, a 2 ) A (u\,u 2 ). Hence (mi,m 2 ) G |cq, /3i] x [a 2 , fi 2 \. 
Step 3: If (mi,m 2 ) is a solution of (3.1.26), then (u\,u 2 ) satisfies (3.1.2). 
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We claim 

(ai(0),a 2 (0)) ^ (mi(0),u 2 (0)) + </> (ui(0),« 2 (0),ui(l),u 2 (l),%(0), « 2 (0)) ^ (ft(0),ft(0)). 

(3.1.27) 

If (ui(0),tt 2 (0)) + </> ^mi(O), w 2 (0), ui(l), it 2 (l), itftO), w 2 (0)) ^ (ft(0), ft(0)), then 

(ui(0),u 2 (0)) = 4>* ^ui(0),u 2 (0),ui(l),u 2 (l),w 1 (0),u 2 (0)j 

= p(0, (ui(0),u 2 (0)) + ftui(O), it 2 (0), ui(l),u 2 (l), Ui(0), u 2 (0))) 

= (ft(0),ft(0)). 

From Step 2, we know that (u\,u 2 ) -< (ft, ft), and this together with 

(«i ~ Pi,u 2 - /3 2 ) G C' 2 [0,1] x C 2 [0,1] and (ui(0),u 2 (0)) = (ft(0),ft(0)) yields %(0) < @[(0) and 
w 2 (0) < ft(0). If 4>p{x,y) is monotone nonincreasing, then we have 

(iti(0),u 2 (0)) + <f> ^ui(0),u 2 (0),ui(1),u 2 (1),u 1 (0),u 2 (0)^ 

= (/5i(0), ft(0)) + <j> ^ft(O), ft(0), «i(1), w 2 (1), Wi(0), u 2 (0)) 

^ (ft (0), ft(0)) + <j> (/3i(0), ft (0), «i(l), u 2 (l), ft (0),^(0)) 

= (/5i(0), ft(0)) + fti(l), u 2 (1)) 

^ (/5i(0), ft(0)) + (ai(l),a 2 (l)) 

= (ft (0), ft(0)) + <j> (ft (0), ft (0), ai(l), a 2 (l), ft (0), ft (0)) 

^ (ft(0), ft(0)), (3.1.28) 

a contradiction. Similarly, if <j>p(x,y) is monotone nondecreasing, then we get same contradiction. 
Consequently, (3.1.27) holds. Similar reasoning shows the other boundary condition. Consequently, 
(ui,u 2 ) satisfies (3.1.2). Hence the system of the coupled BVPs (3.1.1)-(3.1.2) has at-least one 
solution (ui,u 2 ) G [aq,ft] x [a 2 ,ft]. 

3.1.4 Application of the theoretical results 

In this section we study some examples to demonstrate the validity of our theoretical results. 
Example 3.1.1. Consider the following nonlinear second order coupled system 

u 1 (t)=u 2 (t) + sin(f), t G [0,1], 

w 2 (f) = u\ (f) + cos(t), t G [0,1], (3.1.29) 

with the following nonlinear CBCs 

(^1)- U2 (0),uf^i:i)- Wl (0)) =(0,0), n G N, 

(u[(l)-ul n+ \l),u 2 (l)-ul n+ \l))=(0,0), «GN. (3.1.30) 

Choose oq (t) = sin(t) — 5, a 2 (t) = sin(t) — 6 and ft(t) = sin(t) + 5, ft(f) = sin(t) + 6. In view 
of the Definition 3.1.1, it is evident that a couple of functions (aq ,a 2 ) and (ft, ft) are lower and 
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upper solutions of the coupled system (3.1.29), respectively. Further, ( 011 , 0 : 2 ) and (/? i,/3 2 ) satisfies 
the system of inequalities (3.1.25). Hence by Theorem 3.1.2, the coupled system of BVPs (3.1.29)- 
(3.1.30) has at least one solution ( 111 , 1 x 2 ) G [oi,/3i] x [a 2 ,/3 2 ]. 

Example 3.1.2. Consider the following nonlinear second order coupled system 

u ± (t) = u\{t) — 10\/sin(i), t € [0,1], 

u 2 {t) = uf(t) — 8\/cos (t — 1), iG[0,1], (3.1.31) 

with the following nonlinear CBCs 

(ixi(O)ixi(l) - M 2 (0)u 2 (l),tii(0)u 1 (0) - 1 x 2 ( 0 ) 112 ( 0 )) = (0,0), 

(111 (O)^(1) - ix 2 (0)u 2 (1),iii(1)ix'i( 1) - ix 2 ( 1)4(1)) = (0,0). (3.1.32) 

Choose «i (t) = oi 2 (t) = — 3 - and fi\(t) = t + 2, fi 2 {t) = t + 3. In view of the Definition 
3.1.1. it is evident that a couple of functions (oi,o 2 ) and (/?i,/3 2 ) are lower and upper solutions 
of the coupled system (3.1.31), respectively. Further, ( 01 , 02 ) and (/?i,/3 2 ) satisfies the system of 
inequalities (3.1.25). Hence by Theorem 3.1.2, the coupled system of BVPs (3.1.31)-(3.1.32) has at 
least one solution ( 1 x 1 ,u 2 ) G ]oi,/?i] x [o 2 ,/3 2 ]. 

□ 

3.2 Second-order coupled system with nonlinear coupled bound¬ 
ary conditions (II) 


In this section, we study the existence of solutions of the following second-order ordinary nonlinear 
coupled system of BVPs of the type 

—“l (*) = fl{t,U2{t)), t G [0, 1], 

-u 2 {t) = f 2 {t,u 1 (t)), t G [0,1], (3.2.1) 

subject to nonlinear CBCs 

p (iii(O), it 2 (0), iii(O),ix 2 (0),ix 1 (l),ix 2 (l)) = (0,0), 
v (iXi(0), ix 2 (0)) + (iti(l), ix 2 (l)) = (0, 0), (3.2.2) 

where /i,/ 2 :[0,l]xl4l,/x:l 6 -) R 2 and v : M 2 — \ R 2 are continuous functions. We establish 
the existence of at-least one solution of (3.2.1)-(3.2.2) with the help of LUSs method in the presence 
of the existence of well ordered LUSs such that (ai,a 2 ) di (/3i, /3 2 ) on [0,1] x [0,1]. 

The monotonicity assumptions on (f> and 1 f as stated in Theorem 3.1.2 does not address to the 
existence criterion for a very important BVPs like periodic and anti-periodic. But in this study, 
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we introduce a new generalized BCs 3.2.2 that satisfy the existence criterion for a periodic and 
anti-periodic BVPs as a particular case under some nronotonicity assumptions on the arguments of 
the boundary functions p and v as discussed in detail in the introductory chapter of this thesis. 

3.2.1 Lower and upper solutions method 

In this section we recall the idea of LUSs for the nonlinear coupled system (3.2.1). 

Definition 3.2.1. We say that a couple of function ( 01 , 02 ) G C 2 [0,1] x C 2 \ 0,1] is a lower solution 
of (3.2.1) if 

-a[(t) < fi(t,a 2 (t)), iG [0,1], 

-02 (i)</ 2 (i,oi(t)), t e [0,1]. (3.2.3) 

In the same way, an upper solution is a function (Pi,p 2 ) G C 2 [0, 1] x C 2 [0, 1], that satisfies the 
same inequalities in reverse order. For ui,u 2 G C 2 [ 0, 1], we define the set 

= {w 2 G C 2 [0,1] : m(t) < w 2 (t) < u 2 (t),t G [0,1]} . 

3.2.2 Coupled lower and upper solutions 

This section deals to the concept of CLUSs that unifies the treatment of many linear and nonlinear 

BVPs under some nronotonicity assumptions on the arguments of the boundary functions y, and v. 

Definition 3.2.2. We say that a couple of functions ( 01 , 012 ), (Pi,P 2 ) G C 2 [0,1] x C 2 [0,1] are 
coupled lower and upper solutions for the problem (3.2.1) and (3.2.2), if ( 01 , 02 ) is a lower solution 
and (Pi,P 2 ) is an upper solution for the system (3.2.1),(oi, 02 ) ^ (Pi,P 2 ), if «i < pi and 02 < P 2 , 
and 

M/3i(0)^2(0),^(0),/3;(0),^(1),/3;(1)) ^ (0,0) ^ M (oi(0),o 2 (0),o'^O), o^O),o'^l),o^l)), 
MA(0),A(0),^(0),^(0),o;(l),o;(l)) d (0,0) ^ /i(oi(0), o 2 (0), oi(0), a' 2 (0),p’i(l), P 2 (l)), 

(oi(l), o 2 (l)) + 1/ (pi(0),p 2 (0)) = (0,0), 

(/3i(1),/3 2 (1)) + 1/(01 (0),o 2 (0)) = (0,0), 

(oi(l), o 2 (l)) + v (oi(0), o 2 (0)) = (0,0), 

(/5i(l), /3 2 (l)) + v (Pi(0), P 2 (0)) = (0,0). (3.2.4) 

3.2.3 Existence of at least one solution 

Theorem 3.2.1. Assume that a couple of functions ( 01 , 02 ) and (Pi, p 2 ) are coupled lower and 
upper solutions for the coupled BVPs (3.2.1)-(3.2.2). Suppose that y is nondecreasing in the third 
and fourth arguments. In addition suppose that the function v in [oi(0),/?i(0)] x [ 02 ( 0 ),/^(O)] is 
monotone and the functions 

/^(ai, a.2 ){x,y) ■■= y ^oi(O), a 2 (0), a 1 (0), o 2 (0), x, y^j , 

M(/3i ,fa)(x,y) ■= P ^Pi(0),p 2 (0),p'i(0),p' 2 (0),x,yj , 
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have got the same kind of monotonicity as v, then there exists at least one solution ( ui,u 2 ) G 
[an, A] x [a 2 ,A] of the coupled BVP (3.2.1)-(3.2.2). 

Proof. Let A > 0, and consider the following modified coupled BVP 

~u 1 (t) + Aui(t) = F* (t,u 1 (t),u 2 (t)), t G [0,1], 

-u 2 (t) + \u 2 {t) = F 2 (t,ui(t),u 2 (t )), t G [0,1], 

A (ui(0),'U2(0),Mi(l),M2(l),w , i(0),M2( 0 )) = («i(0), w 2 (0)), 

(ui(l),« 2 (l)) + («i(0),«2(0)) = (0,0), (3.2.5) 


with 


F*{t,ui{t),u 2 (t)) = < 


fi (LAW) + AA W 

fi (t,u 2 (t)) + AA (t) 

fi (t,a 2 (t)) + AA W 
fi (LAW) + Awi(t) 

fi (t,u 2 (t)) + Aui (t) 


fi {t,a 2 (t)) + Xui(t) 

fl (L AW) + A«i(t) 
fi (t,u 2 (t)) + Aai(t) 
fi (t,a 2 (t)) + A«i(t) 


if (m(t),u 2 (t)) ^ (AW, AW), 
if a 2 (t) < u 2 (t ) < fi 2 (t),ui(t) > AW, 
if u 2 (t) < a 2 (t),ui(t) > AW, 
if u 2 {t) > AW, a i(t) < ux{t) < AW, 
if a 2 (t) < u 2 (t) < AW, 

Oil W < Ui(t) < /3i(t), 
if u 2 (t) < a 2 (t),ai(t) < ui(t) < AW, 
if u 2 (t) > AW,«iW < ai(t), 

if a 2 (t) < u 2 (t) < /3 2 (t),ui(t) < ai (A), 
if u 2 (t) < a 2 (t),ui(t) < ai(t), 


and 


F 2 (t,ui(t),u 2 (t)) = < 


h (LAW) + AAW 
h (L«iW) + AAW 
f 2 (t,ai(t)) + AAW 
/ 2 (LAW) + A u 2 (t) 
f 2 (t,ui(t)) + A u 2 (t) 


/ 2 (t,ai(t)) + A u 2 (t) 

/ 2 (LAW) + A a 2 (t) 
/ 2 (L u i(t)) + Aa 2 W 
h (t,ai(t)) + Aa 2 W 


if («iW,u 2 W) ^ (AW,AW), 
if aiW < tiiW < AW,w 2 W > AW, 
if uiW < aiW,u 2 W > AW, 
if uiW > AW,“ 2 W < « 2 W < AW, 

if aiW < tiiW < AW, 

a 2 {t) < u 2 {t) < AW, 
if uiW < an W, a 2 W < m 2 W < AW, 
if uiW > AW,m 2 W < a 2 (t), 
if aiW < ui (t) < (3i(t),u 2 (t) < a 2 (t), 
if uiW < aiW,u 2 W < a 2 W, 


H* U,k,l,m,n,o) = <j> (0, (j, k) + fi (j, k, l, m, n, o )), 

A(j,fc) = u(<j>(0,(j,k))), 

( (AW, AW) if (x,y) (A, A) 

(x,y)) = l (x,y) if (ai,a 2 ) ^ (x,y) ^ (A, A) 

[(aiW,a 2 W) if (x,y) ^ (ai,a 2 ), 

Note that if (ui,u 2 ) G [ai,A] x [a 2 , A] is a solution of (3.2.5), then («i,u 2 ) is a solution of 
(3.2.1)-(3.2.2). 
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For the sake of simplicity we divide the proof in three steps: 
Step 1: We define the mappings 


by 


L,N : C 2 [ 0,1] x C 2 [ 0,1] -► C&[ 0,1] x C$[0, 1] x f x M , 


[L(ui, U 2 )](t) = u 1 (t) — u 1 ( 0) — A / ui(s) ds, u 2 (t) — u 2 (0) — A / u 2 (s)ds, 


(«i(0),u 2 (0)),(ui(l),u 2 (l)) ), 


and 


[N(ui, u 2 )](i) = (/ F? { s,u l{ s),u 2 (s))ds,J F2(s,u l{ s),u 2 (s))ds, 

H* («i(0), u 2 (0), ui(l), u 2 (l), «i(0), u 2 (0)), 

-i/*(ui(0),u 2 (0))V 


Since (s, Ui(s), u 2 (s)) and F£{s, Wi(s), w 2 (s)) are bounded on [0, l]xK 2 and integral is a continuous 
function on [0,1]. Further /i* and z/* being constant functions are continuous. Therefore [IV(ui,u 2 )] 
is continuous on [0,1]. Further, the class {N{ui,u 2 ) : ui,u 2 G C 2 [0,1]} is uniformly bounded 
and equicontinuous. Therefore in view of Arzela-Ascoli theorem {N{ui,u 2 ) : ui,u 2 G C 2 [0,1]} is 
relatively compact. Consequently A” is a compact map. Also from Lemma 3.1.1 with a = 1,6 = 0, 
c = 1, d = 0 and E = 0, F = 1, G = 0, H = 1, A -1 , exists and is continuous. 

On the other hand, solving (3.2.5) is equivalent to find a fixed point of 

L~ 1 N : C 2 [0,1] x C 2 [ 0,1] C 2 [ 0,1] x C 2 [ 0,1], 

Now, Schauder’s fixed point theorem guarantees the existence of at least a fixed point since L~ l N 
is continuous and compact. 

Step 2: If («i,u 2 ) is a solution of (3.2.5), then {u\,u 2 ) G [ai,/3i] x [a 2 ,/3 2 ]. By definition of //, 
we see that (wi(0), w 2 (0)) G [ai(0),/?i(0)] x [a 2 (0),/3 2 (0)]. Thus, if v is nondecreasing, we have by 
condition (3.2.4) 

(ai(l),a 2 (l)) = —v (/?i(0), /? 2 (0)) ■< -v (ui(0), it 2 (0)) = (iti(l),u 2 (l)) ■< -v{a i(0),a 2 (0)) 
(ai(l),a 2 (l)) A (ui(l),u 2 (l)) A (A(l),/3 2 (1)), (3.2.6) 

similarly, if u is nonincreasing, then (3.2.6) holds. Hence («i(l),u 2 (l)) G [ai(l), /3i(l)] x [a 2 (l), /3 2 (1)]. 
Now, it remains to show that (wi,w 2 ) G [ai,/3i] x [a 2 ,/3 2 ] for t G (0,1). 

We claim {ui,u 2 ) r< {(3i,f3 2 ). If {ui,u 2 ) ^ {fii,/3 2 ), then u\ ^ /3i and/or u 2 ^ f3 2 . If Ui ^ (3 1 , 
then there exists some t 0 G [0,1] such that Ui{t 0 ) — f3i{t 0 ) > 0. So, U\ — fix attains a positive 
maximum at t 0 G [0,1]. Thus (u\ — fi\) (to) = 0 and (ui — (3 1 ) (t 0 ) < 0. But, 

(id - Pi)" {t 0 ) > —F* (t 0 ,ui(to),u 2 {t 0 )) + Aui(to) + fi {to,fa{t 0 )) 

= —fi {to, p 2 {to)) — A/3i(to) + A«i(to) + fi {to, ^2 (to)) 

= A (ui(to) - Pi{to)) > 0, 
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a contradiction. Similarly one can show that (<* 1 , 0 : 2 ) ^ (ui,u 2 ). Hence (ui,it 2 ) G [oi,/0i] x [a 2 ,/3 2 ]. 
Step 3: If (u\,U 2 ) is a solution of (3.2.5) then (ui,U 2 ) satisfies (3.2.2). 

We claim 

(O!(0),o 2 (0)) ^ ( Wl (0), w2 (0))+ / r(u 1 (0),« 2 (0),u , 1 (0), u ;(0), u ' 1 (l), w ;(l)) ^ (ft(0),/3 2 (0)). (3.2.7) 
If (oi(0),o 2 (0)) (ui(0),u 2 (0)) + p (u 1 (0),u 2 (0),u 1 (0), u 2 (0), <4(1).^ 2 ( 1 )); then 

(ui(0),w 2 (0)) = p* (ui(0),u 2 (0),u 1 (0),u 2 (0),u' 1 (1),u 2 (1)) 

= 0(0, («i(0), u 2 (0)) + p{u\ (0), m 2 (0),u 1 (0),u 2 (0), Ui(1), « 2 (1))) 

= (oi(0),o 2 (0)). (3.2.8) 

If v is nondecreasing, we have 

(ui(l), w 2 (l)) = v («!(0), U 2 (0)) = -v (oi(0),a 2 (0)) = (/3i(l),/3 2 (1)). (3.2.9) 

Using (3.2.8),(3.2.9) and Step 2, we have ^u 1 (0),u 2 (0)^ y ^oq(O), a 2 (0)^ and ^(1), w 2 (l)^ y 
(/3i( 1),^ 2 (1)). But 

(ui(0),u 2 (0)) + p (u 1 (0),u 2 (0),m , 1 (0),u 2 (0),u' 1 (1),u 2 (1)) 

= (oi(0),o 2 (0)) + p ^o 1 (0),o 2 (0),u 1 (0),u 2 (0),u 1 (1),u 2 (1)^ 

h (oi(0),o 2 (0)) + p ^o 1 (0),o 2 (0),o 1 (0),o 2 (0),u 1 (1),u 2 (1)^ 

= (Oi(0), O 2 (0)) + P( ai ,a 2 ) (% (1), U 2 (l)) 

h («i(0),o 2 (0)) + p (o 1 (0),a 2 (0),a 1 (0),o 2 (0),^ 1 (l),^ 2 (l)) 
h (ai(0),o 2 (0)), (3.2.10) 

a contradiction. Similarly if v is nonincreasing we get same contradiction. Consequently, (3.2.7) 
holds. By definition of v* and Step 2, the second boundary condition is obvious. Consequently 
(wi,w 2 ) satisfies (3.2.2). Hence the coupled system of BVPs (3.2.1)-(3.2.2) has at least one solution 
(wi,w 2 ) G [ai,/3i] x [a 2 ,/3 2 ]. □ 


3.2.4 Application of the theoretical results 

Example 3.2.1. Consider the following nonlinear coupled system of BVP 

— u 1 (t) = ul(t) — 20-\/sin(t), t G [0,1], 

— u 2 (t) = u\{t) — 16-\/cos(t — 1), t G [0,1], 

with the following nonlinear CBCs 


(u 1 (0)u , 1 (0) - w 2 (0)u 2 (0), u 1 (0)m / 1 (1) - w 2 (0)u 2 (l)) = (0,0), 

^V / tan(«i(0))'Ui(l) + tan(u 2 (0))it 2 (l), \/tan(wi(0)) + tan(w 2 (0))w 2 (l)^ = (0, 0). 


(3.2.11) 


(3.2.12) 
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Choose ai (t) = —^ 7 . a 2 (i) = — y and /3i(0 = p 2 {t) = y. In i/ie light of the Definition 

3.2.1, It is evident that a couple of functions ( 01 , 02 ) and (/ 3 i,/? 2 ) are lower and upper solutions of 
the coupled system (3.2.11), respectively. Further, ( 01 , 02 ) and (/?i,/? 2 ) satisfy (3.2.4). Hence by 
Theorem 3.2.1, the coupled system of the BVPs (3.2.11)-(3.2.12) has at least one solution ( 161 , 1 ( 2 ) € 
[oi,/3i] x [a 2 ,/3 2 ]. 


Example 3.2.2. Consider the following nonlinear coupled system of BVP 

- u 1 (■ t ) = t 2 + sin (u 2 (t)), t € [0,1], 

— u 2 (t) = 5 1 2 + cos(iii(0 + 4) + sin(ui(f) — 4), t € [0,1], 

with the following nonlinear CBCs 

(ni(O)ui(l) - n , 2 ( 0 )n , 2 (l),n / 1 ( 0 )'u , i(l) - n 2 (l)) = ( 0 , 0 ), 
(sin(«i( 0 )) + 111 ( 0 ) 112 ( 1 ),sin(u 2 ( 0 )) + u 2 ( 0 )ui(l)) = ( 0 , 0 ). 


(3.2.13) 


(3.2.14) 


Choose Oi (t) = t 2 , 02 (t) = t 3 and = t, fait) = t. In the light of the Definition 3.2.1, It is 
evident that a couple of functions ( 01 , 02 ) and are lower and upper solutions of the system 

(3.2.13), respectively. Further, ( 01 , 02 ) and satisfy the coupled system (3.2.4). Hence by 

Theorem 3.2.1, the system of BVPs (3.2.13)-(3.2.14) has at least one solution ( 1 ( 1 , 112 ) €E [ai,/3i] x 

[a 2 ,/3 2 ]. 



Chapter 4 

Nonlinear Coupled Boundary 
Value Problems With Derivative 
Dependent Nonlinearity 


In this chapter, we extend the results of Chapter 3 for a nonlinear second-order coupled BVPs 
subject to nonlinear CBCs with dependence of the nonlinearity on derivative. The existence results 
are established assuming the existence of well ordered LUSs, such that (on, 012) ^ (/3i, P2), if ol\ < /?i 
and «2 < P 2 on [0,1] under Nagurno conditions. The classical LUSs approach is adopted to ensure 
the existence of at least one solution in the closed interval defined by a well ordered LUSs. 

The study of BVPs with nonlinearity dependent on the derivative with certain type of BCs was 
initially studied by Schmitt in [167] and developed the classical results to ensure the existence of 
solutions of these types of BVPs using LUSs approach under Nagumo condition. Currently, extensive 
study has been done on these types of BVPs including, Dirichlet, periodic, and multi-point BCs, 
see for example [55, 80, 91, 110, 111, 166, 175, 189, 194] and Chapter III of [102]. In [13], Chapter 
III, IV, the new classical results are established with multi-point linear coupled BCs for the case 
when the derivative dependent nonlinear terms are allowed to be a singular at the end points of 
the closed interval [0,1]. In contrast of these studies we develop the existence results for derivative 
dependent nonlinear coupled BVPs corresponding to a generalized nonlinear CBCs in the presence 
of the Nagumo condition. 

This chapter is organized as: In Section 4.1.1, the idea of LUSs is recalled for the second-order 
nonlinear coupled BVPs having dependence of the nonlinearity on the first-order derivative. In 
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Section 4.1.2, the Nagumo conditions are presented that ensure the existence of the solutions of the 
BVPs having dependence of the nonlinearity on the derivative. In Section 4.1.3, the idea of CLUSs 
is discussed. In Section 4.1.4, the existence of at-least one solution is ensured using LUSs method. In 
Section 4.1.5, an example is considered to check the validity of the developed theoretical results. In 
Section 4.2.1, the concept of CLUSs is discussed for the more generalized BCs. In Section 4.2.2, the 
LUSs methodology is used to investigate the existence of at-least one solution for the second-order 
BVPs with dependence of the nonlinear terms on the first-order derivative. In Section 4.2.3, an 
illustrative example is discussed to check the validity of the developed theoretical result. 

4.1 Second-order coupled system with derivative dependent 
nonlinearity (I) 

In this section, we discuss the existence of solutions of the following second-order nonlinear coupled 
system of BVPs having nonlinearity in the first-order derivative of the type 

~ul(t) = fi(t,u 1 (t),u 2 (t),u' 1 (t),u 2 (t)), t€ [0,1], 

-u'^it) = f 2 (t,U 1 (t),U 2 (t),U 1 (t),U 2 (t)), t £ [0,1], (4.1.1) 

< ('(' u i(0), u 2 (0), ui(l), u 2 (l), 1 / 1 ( 0 ), i4(0)) = (0,0), 

^(mi(0),U2(0),Ui(1),u 2 (1),u , 1 (1),U2( 1 )) = (0,0), (4.1.2) 

where / 1 , f 2 : [0,1] x R. x R. x R x R — > R, and <fi, ip : R 6 — > R 2 are continuous functions. 

The BCs (4.1.2) are very much generalized in the sense that they generalize most of the linear and 
nonlinear BCs under some monotonicity assumptions on the arguments of the boundary functions 
< f> and if) as discussed in detail in the introductory chapter of this thesis. Moreover (4.1.2) handle 
the classical existence criterion for the Dirichlet and Neumann BCs as a particular case. 

4.1.1 Lower and upper solutions method 

In this section, we recall the concepts of LUSs for the nonlinear coupled BVPs of the type (4.1.1) 
and Nagumo condition. 
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Definition 4.1.1. We say that a couple of function ( 0 : 1 , 02 ) G C 2 [0,1] x C 2 { 0,1] is a lower solution 
of (4.1.1) if 

—«i (*) < /1 (t,a 1 (t),a 2 {t),a 1 (t),a 2 (t^J , i G [0,1], 

—“2 (*) < h (t,on(t), a 2 (t),a' 1 (t), a 2 (t)) , [0,1]. (4.1.3) 

In the same way, an upper solution is a function (fii, f3 2 ) G C 2 { 0,1] x C 2 [0,1], that satisfies the 
same inequalities in reverse order. For u\,u 2 £ C 2 [ 0,1], we define the set 

lui,u 2 j = {w 3 G C 2 [0,1] : ui(t) < w 3 (t) < u 2 (t),t G [0,1]} . 

4.1.2 Nagumo condition 

In 1954, Nagumo [137] observed that existence of LUSs is not sufficient to ensure the existence of 
solutions of Dirichlet BVPs when nonlinear terms involving the first-order derivative of the solution. 
This is due to the fact that nonlinear term grows very rapidly with respect to the gradient of the 
solution. Later on Habets and Pouse [87] observed the same problem for periodic and separated 
BVPs and found that method of LUSs is not valid without considering the Nagumo condition. 
Therefore to ensure the existence of solutions of the system (4.1.1), we need to find the prior bounds 
on u'i and u 2 . 

Now we define the concept of Nagumo condition for the nonlinear second order coupled system 
(4.1.1). 

Definition 4.1.2. We say that /1 and f 2 satisfies a Nagumo condition relative to the intervals 
[ai(t),/3i(f)] and [a 2 (t), fi 2 (t)\ respectively, if for 

Po = max {/?i(0) - «i(l) +/3 2 (0) - a 2 (l), /3i(l) - «i(0) + /3 2 (1) - «2(0)} , 

there exists a constant N t such that 


Ni > max < po, sup |a-(t)|, sup |/3-(t)| > , i = 1,2, 
( te [ 0 , 1 ] te[o,i] J 

and a continuous function g\ : [0, 00 ) —> (0, 00 ) such that 

\fi{t,u 1 {t),u 2 {t),d,e\ < gi{\d + e|), t G [0,1], 
\f 2 (t,u 1 (t),u 2 (t),d,e\ < g 2 (\d + e\), t G [0,1], 


whenever (m (t),a 2 (t)) ^ (ui(t),u 2 (t)) < /3 2 (t ,)), d, e G R, and 



1 

9i(y) 


dy > 2. 


(4.1.4) 
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4.1.3 Coupled lower and upper solutions 

In this section we define the concept of CLUSs that unify the treatment of many linear and nonlinear 
BVPs by imposing some monotonicity assumptions on the arguments of the boundary functions (p 
and ip. 

Definition 4.1.3. We say that a couple of functions ((on, a 2 ) and {pi,p 2 )) G C 2 [0,1] x C 2 [0,1] are 
coupled lower and upper solutions for the coupled BVP (4.1.1)—(4.1.2), if (on, < 32 ) is a lower solution 
and (/?i,/ 3 2 ) is an upper solution for the system (4.1.1) with (ai,a 2 ) p {Pi,p2), if «i < Pi and 
a 2 < P 21 such that 

ct>(Pi(0),P2(0),Pi(l),p2(l),P' 1 (0),P' 2 m p (0,0) ^ P( ai ( 0),a 2 (0),a 1 (l),a 2 (l),a' 1 (0),a;(0)), 
^(/3 1 (0),/3 2 (0),a 1 (l),a 2 (l),/3;(0),/3;(0)) ^ (0,0) ^ 0(ai(O),a 2 (0),ft(l), p 2 (l), a^O), a 2 (0)), 
^(/3i(0),^ 2 (0),/3i(1),/3 2 (1),/3 1 (1),/3 2 (1)) ^ (0,0) ^ 0(01(0),o 2 (0),ai(l),a 2 (l),Oi(l),a 2 (l)), 
■0(01(0),a 2 (0),/3i(l),/3 2 (l),^ 1 ( 1 ),/3 2 (l)) ^ (0,0) ^ 0(A(O),/3 2 (0), ai(l), o 2 (l), Oi(l), a 2 (l)). 

(4.1.5) 


Now we state and prove theorem that establishes the existence of at least one solution for BVPs 
(4.1.1)-(4.1.2). 

4.1.4 Existence of at least one solution 

Theorem 4.1.1. Suppose the nonlinear coupled BVP (4.1.1)-(4.1.2) having (cu, o 2 ) and (P\,p 2 ) 
coupled lower and upper solutions respectively and the nonlinear functions f\ and f 2 relative to the 
intervals [ 01 (f), Pi(t)\ and [o 2 (f), p 2 (f)] satisfy the Nagumo assumptions respectively. Moreover, as¬ 
sume that ip and (p are non-increasing and non-decreasing in the fifth and sixth arguments respectively 
and 

fpa{x,y) := ip [x,y,a i(l), o 2 (l), o^l), o 2 (l)) , 
tpp{x,y) ■= *P (x, y, /5i( 1 ), /? 2 ( 1 ), /3i(f), /3 2 (i)), 

are monotone on ]oi(0),/?i(0)] x [a 2 (0),/? 2 (0)] and 

ipa(x,y) ■■= ip (ai( 0 ),a 2 ( 0 ),a:, 2 /,a 1 ( 0 ),a 2 ( 0 )^ , 

*Pp{x,y) ■■= Ip (^pi(0), p 2 {0),x,y, P^O), P 2 (0) S j , 
are monotone on [oi(l), /?i(l)] x [a 2 (l),/? 2 (1)]. 

Then the nonlinear coupled system of BVPs (4.1.1)-(4.1.2) have at least one solution (ui,u 2 ) € 
[ai,/?i] x [a 2 ,/J 2 ] and (-Vi,-7V 2 ) p (u[{t),u 2 (t)) p {Nx,N 2 ),t G [0,1]. 
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Proof. Let A > 0, and let the modified coupled BVP 

-w'i'(t) + Aui(t) =A (t,u 1 (t),u 2 {t),u 1 (t),u 2 (t) S j , t S [0,1], 

-w 2 (<) + A** 2 (t) = F 2 ^t : u 1 (t),u 2 {t),u 1 (t),u 2 (t)^ , t e [0,1], 

cj) (ui(0),u 2 (0),ui(l),u 2 (l), %(()), *t 2 (0)) = (**i(0), ** 2 (0)), 

xj) (ui(0),u 2 (0),wi(1),u 2 (1),u' 1 (1),u 2 (1)) = («i(l), u 2 (l)), (4.1.6) 

with 

Fi (t,ui{t),u 2 (t),a(t),b{t)) = /i (t,T 1 (t,u 1 (t)),T 2 (t,u 2 (t)),a 1 (a(t)),a 1 (b(t))) + Ari(<,**i(f)), 
i ?2 (f,**i(f),?* 2 (t),a(f),&(f)) = / 2 (f, Ti(t, Ui(t)),T 2 (t, u 2 (t)) : tii(a(t)), a 2 (b{t))) + Ar 2 (f, *t 2 (t)), 


(f>(j,k,l,m,n,o) = (n(0,j),r 2 (0,/c)) + <M«i(0 ),u 2 (0 ),ui(1),m 2 (1 ),w , 1 (0 ),u 2 (0)) , 
xp(j,k,l,m,n,o ) = (ri(l,Z),T 2 (l,m)) + V’(wi(0),u 2 (0),ui(1),u 2 (1),u , 1 (1),u 2 (1)) , 

fft(i) ifx^/3j(t), * = 1,2 

7 y(f, z) = < x if a.i(t) < x < /3j(t), *=1,2 
[a*(t) if Xjfaii(t), * = 1,2 

and 

ifs^iVi, * = 1,2 

(T,:(a;) = < x ii — Ni < x < N. t , * = 1,2 

ifz^-TVi, * = 1,2. 

Note that if (**i,?* 2 ) € [cn,/3i] x [a 2 ,/3 2 ] is a solution of (4.1.6), then (**i,*t 2 ) is a solution of 
(4.1.1)-(4.1.2). 

For the sake of simplicity we divide the proof in four steps: 

Step 1 : We define the mappings 

L, N : C 2 [ 0,1] x C 2 [ 0,1] Cg[0,1] x C o 2 [0,1] x R 2 x R 2 , 


by 


[L(ui,u 2 )]{t) 



** 1 (0) — A / ?*i(s) ds, ** 2 (t) 
Jo 


** 2 (0) — A / *t 2 (s) ds, 
Jo 


(**i(0),** 2 (0)) ,(**i(l),*t 2 (l))^, 


and 


[A/(**i,** 2 )](t) = yj F 1 (s,ui(s),u 2 (s),u 1 (s),u 2 (s))ds, J F 2 (s,u 1 (s),u 2 (s),u 1 (s),u 2 {s))ds, 

</>(**! (0), *t 2 (0), **i (1), u 2 (1), u[ (0), u 2 (0)), 

V’(**i(0),u 2 (0),ui(l),u 2 (l), u^l), *t 2 (l))V 



76 


Since Fi(s,ui(s),u 2 (s),u 1 (s),u 2 (s)) and -F 2 (s, Mi(s), m 2 (s), m^s), m 2 (s)) are bounded on [0,1] x M 4 
and integral is a continuous function on [0,1]. Further </> and i/j being constant functions are contin¬ 
uous. Therefore [N(ui,u 2 )] is continuous on [0,1]. Further, the class {N(ui,u 2 ) : m x ,m 2 € C 2 [0,1]} 
is uniformly bounded and equicontinuous. Therefore in view of Arzela-Ascoli theorem {N(ui, u 2 ) : 
m x ,m 2 € C 2 [0,1]} is relatively compact. Consequently N is a compact map. Also from Lemma 3.1.1 
with a = 1, b = 0, c = 1, d = 0 and E = 0, F = 1, G = 0, H = 1, L _1 , exists and is continuous. 

On the other hand, solving (4.1.6) is equivalent to find a fixed point of 

L- X N : C 2 [ 0,1] x C 1 [0,1] C^O.l] x C^O.l]. 

Now, Schauder’s fixed point theorem guarantees the existence of at least a fixed point since L~ l N 
is continuous and compact. 

Step 2: If (ui,u 2 ) is a solution of (4.1.6), then (i(i,it 2 ) G [aq,/3i] x [a 2 ,/3 2 ]. By definition of 
( j ) and ip , it is obvious that, (ui(0), m 2 (0)) G [«i(0),/?i(0)[ x [a 2 (0),/? 2 (0)] and (ui(l), m 2 (1)) G 
[ai(l),^i(l)[ x [a 2 (l),/3 2 (l)] respectively. We claim (ui,m 2 ) ^ (Pi,p 2 ). If (ui,u 2 ) ^ (/3i,/3 2 ), then 
Mi ^ /3i and/or u 2 ^ /? 2 . If Mi ^ /3 X , then Mi — f3\ has positive maximum at some r 0 G (0,1) such 
that (mi — /3i) (r 0 ) = 0 and (mi — pi)" (ro) < 0. But, 

(mi-/3i) (r 0 ) > -F (r 0 ,M 1 (ro),M 2 (ro),M 1 (ro),M 2 (r 0 )^ +AMi(r 0 ) + /i (r 0 ,/?i(r 0 ),&(r 0 ), ^(ro),/3 2 (r 0 )) 
= ~fi (ro,Pi{ro),p2{ro),u 1 (r 0 ),u' 2 {r 0 ) S j - A/3i(r 0 ) + AMi(r 0 ) 

+ fi (r 0 , Pi (ro ), P 2 (r 0 ) , Pi (r 0 ), P 2 (r 0 )) 

= ~fi (ro, Pi(ro), p 2 (r 0 ), Pi(ro), P 2 ( r o)^ - A/3i(r 0 ) + Aui(r 0 ) 

+ h (ro , Pi {r 0 ), P 2 {ro), P[ (r 0 ), P 2 (r 0 )) 

= A (Mi(r 0 ) - pi (r 0 » > 0, 

a contradiction. Similarly we can show that (a i,a 2 ) ;< (u,v). Hence (mi,m 2 ) G [ai,/3i] x [a 2 ,/I 2 ]. 

Step 3: If (mi,m 2 ) is a solution of (4.1.6), then (ui,u 2 ) satisfies (4.1.2). 

For this, we claim that 

(«i(0),a 2 (0)) ■< (mi(0),m 2 (0))+</> (^mi(0),m 2 (0),Mi(1),m 2 (1), iq(0),M 2 (0)^ ■< (/3i(0), #>(0)). (4.1.7) 

If (mi(0),m 2 (0)) + <j> ^mi(0),m 2 (0),mi(1),m 2 (1),m' 1 (0),M2(0)^ ■£ (/3 X (0),/3 2 (0)), then 

(mi(0),m 2 (0)) = ( j ) ^mi(0),m 2 (0),Mi(1),m 2 (1),m 1 (0),m 2 (0)^ 

= (ti(0, mi(0)), t 2 (0,m 2 (0))) + </>(mi(0), m 2 (0), tq(l), m 2 (1), m^O), m 2 (0))) 

= (/3 1 (0),/3 2 (0)). 


From Step 2, we know that (ui,m 2 ) ^ (Pi,p2), and this together with 

(mi ~ P 1 ,u 2 - p 2 ) G C 2 [0,1] x C 2 [0,1] and (mi(0),m 2 (0)) = (/3 X (0),/3 2 (0)) yields u'^O) < p[(0) and 
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it 2 (0) < ^ 2 ( 0 )• K <!>p(x,y) is monotone nonincreasing, then we have 

(mi(0),m 2 (0)) + < t > (■Ui(0),m 2 (0),Wi(1),m 2 (1),m , 1 (0),u 2 (0)) 

= (0i (0), # 2 ( 0 )) + 0(A(O), ^ 2 ( 0 ),ui( 1 ),u 2 ( 1 ),u 1 ( 0 ),u 2 ( 0 )^ 

^ (0,(0), fo(0)) + </> (/3i(0), 02(0), «i(l), « 2 (1), 0[(O), ^(0)) 

= (0l(O),02(O)) + 4>p(U!(l),U2(l)) 

^ (0i(O),02(O)) + (f>p(a 1 (l),a 2 (l)) 

= (0i(O), 02(0)) + 0 (0i(O), 02(0 ), ai(l), a 2 (l), ^(0), ^(0)) 

^ (0,(0), 02(0)), (4.1.8) 

a contradiction. Similarly, if cj)p(x,y) is monotone nondecreasing, then we get same contradiction. 
Consequently, (4.1.7) holds. Similar reasoning shows the other boundary condition. Consequently, 
(wi,« 2 ) satisfies (4.1.2). Hence the system of BVPs (4.1.1)-(4.1.2) has a solution (u,,u 2 ) £ [aq,/?i] x 
[a 2 ,/3 2 ]- 

Step 4 : If («i,« 2 ) £ C 2 [0,1] x C 2 [0,1] is a solution of (4.1.6) then (—N,,—N 2 ) -< (u',(t) ,w 2 W) "S 

(N,,N 2 ). We claim -< (N,,N 2 ). If (u'-^u^ -ft (N,,N 2 ), then u, N, and/or u' 2 N 2 . If 

u, N,, then there exists Sq € [0,1] such that u,(so) > N,. Moreover using the Lagrange theorem, 
there exists t 0 £ [0,1] with iq(to) = «i(l) — «i(0). So, 

-Ni < -p 0 < ai(l)-0 1 (O)+a 2 (l)-02(O) < u 1 (to)+n 2 ( i o) < /3i(0)-ai(0)+/3 2 (0)-a 2 (0) < p 0 < N,. 

Now consider an interval [si,s 2 ] or [s 2 ,Si] such that tti(si) + w 2 (si) = po and u 1 (s 2 ) + u 2 (s 2 ) = N,, 
either 

Po = Wi(si) +u 2 (si) < %(*) + u 2 (t) < u,(s 2 ) +u 2 (s 2 ) = N,, te(si,s 2 ), or, 

/ / 111 / (4.1.9) 

Po = %(si) +u 2 (si) < u,(t) + u 2 (t) < u,(s 2 ) +u 2 (s 2 ) = N,, t £ (s 2 ,s 1). 


In the first situation we obtain from Definition 4.1.2 that 

,i*;(8 2 )+t4(* 2 ) dy _ rN 1 dy 

j (si)+u2(si) diy') J pq 

Using (4.1.4), Step 2 and N, > u 1 (t) + u 2 (t) > po > 0, for all t £ (s,,s 2 ), we get the contradiction 


rU,(s 2 )+U 2 (s 2 ) 

I u, (si)+ u((si) 9(y) 


r 2 u,(t)+u 2 (t) dt 
L, g(u’i(t) +u' 2 (t)) 

r __ dt+ r _^_ dt 

Is, sKW+^W) Js, g(u\(t) + u 2 (t)) 

rs 2 —F, (t, U\(t ), U 2 (t ), U, (t ), U 2 (£)) + Xu,(t) 

Is, 9 (u'i(t) + u' 2 (t)) 

r s 2 -F 2 {t,ui(t),u 2 (t),u l (t),u 2 (t)^ + A u 2 (t) 

Is, g(u\(t) + u 2 (t)) 

r s 2 -/l (t,Ul(t), u 2 (t), u,(t),u 2 (t)} rs 2 -f 2 {t,U\(t),U 2 (t),U,(t),U 2 

Is, g(u'i(t) + u 2 (t)) Js, g(u'i(t) + u 2 (t)) 
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< 


< 


'Si 

rS2 


'si 

rS2 


l/l (t,u 1 (t),u 2 (t),u 1 (t),u' 2 (t)^ I rs 2 |/ 2 (t,U 1 (t),U 2 {t),U 1 {t),U 2 {t)^j 

g( u i(t) + u' 2 {t)) J 81 

giu'iit) + u 2 {t)) dt + j~ S2 gju^jt) + u 2 (t)) dt 


g(u'i(t) + u' 2 (t)) 


dt 


g(u'i(t) + u 2 (t)) 


Sl sKW + <4(t)) 


/ dt+ dt 
J S l J Si 

2(s 2 — Si) < 2. 


Similarly in the second situation we get a contradiction. Hence ^w 1 (t), u 2 {tyj -< (Ni,N 2 ). The proof 
of other inequality is similar. □ 


4.1.5 Application of the theoretical results 

In this section, we discuss the practical applicability of our theoretical results with the help of a 
example. 

Example 4.1.1. Consider the nonlinear coupled system of BVPs 

-mi (t) = - (f 2 + 3) (e ui(t) + e“ 2(t) ) - 2 (Vi) 2 + (m 2 ) 2 ) + 2 (t 2 + l) , t £ [0,1], 

-u '2 (t) = - (f 4 + 9) (e 2ui(t) + e 2u2(t) ) - 4 ((m)) 4 + (w 2 ) 4 ) + 4 (t 4 + 2) , te [0,1], (4.1.10) 

subject to nonlinear CBCs 

(mi(0) sin (u 2 (0)) + sin ^m(O)) u 2 (0),ui(0) sin (m 2 (0)) + Mi(0)u 2 (0)) = (0, 0), 

(«i(0) tan (m 2 (0)) - sin (mi(1)) u 2 (1), -Mi(1) tan (n 2 (l)) + m 2 (0)mi(1)^ = (0, 0). (4.1.11) 

The nonlinearities fi(t, xo, xi, x 2 , xf) = — (t 2 + 3) (e x ° + e xi ) — 2 ((® 2 ) 2 + ( *3 ) 2 ) + 2 (t 2 + l) and 
f 2 (t, xo,xi,x 2 ,X 3 ) = — (t 4 + 9) (e 2x ° + e 2xi ) — 4 (fx 2 ) 4 + (X 3 ) 4 ) + 4 (t 4 + 2) are continuous on [0,1] x R 4 . 
Then the functions (ai,a 2 ) , (/3i,/ 3 2 ) : [0,1] x [0,1] -> I x 1 defined by 
ai(t) = —2, a 2 (t) — —4, /3i(t) = t 2 , and f) 2 (i) = t 4 
are, respectively, lower and upper solutions of system (4.1.10), according to Definition 4-1.1. 

Moreover, it can be easily check that the functions a\{t) = —2, a 2 (t) = —4, fii(t) = t 2 , and fd 2 (t) = t 4 
satisfy system of inequalities (4.1.5). Consequently, these functions are, respectively, coupled lower and upper 
solutions of (4.1.10)-(4.1.11), according to Definition 4-1-3. 

Also notice that the functions fi(t, xo, xi, x 2 , X3) = — (t 2 + 3 ) ( e x ° + e xi ) — 2 ((x 2 ) 2 + (*3) 2 ) +2 (t 2 + l) 
and f 2 (t, xo, xi, x 2 , *3) = — (t 4 + 9) (e 2 * 0 + e 2xi ) —4 ((ir 2 ) 4 + (*3) 4 ) +4 (t 4 + 2) satisfies the Nagumo condi¬ 
tion (4.1.4) with g[x 2+X3) = 4 (1 — 2e)—2 ((x 2 ) 2 + (X3) 2 ) , x 2 , X3 £ R and g(x 2 +X3 ) = 4 (3 — 5e 2 — (*2) 4 — (*3) 4 ) 
with x 2 ,X3 £ R, respectively. 

Hence, by Theorem 4-1-1, there is at least one solution (ui(t),u 2 (t)) of the BVPs (4.1.10)-(4.1.11) such 
that, for every t £ [0,1], 

(-2,-4) X ( Ml (t), M2 (t)) X (t 2 ,t 4 ). 
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4.2 Second-order coupled system with derivative dependent 
nonlinearity (II) 

In this section, we investigate the existence of solutions of nonlinear coupled system of BVPs of the 
type 

= f 1 (t,u 1 (t),u 2 (t),u 1 (t),u' 2 (t)), t<= [0,1], 

-U 2 (t) = f 2 (t,u 1 (t),u 2 (t),u 1 (t),u 2 (t)), t € [0,1], (4.2.1) 

/x(mi(0),u 2 (0),u 1 (0),u 2 (0),w'i(1),U2(1)) = (0,0), 

y(«i(0),«2(0)) + (ui(l),«2(l)) = (0,0), (4.2.2) 

where / 1: f 2 : [0,1] x R. x R. x R x R — > R, /.t : R 6 — > I 2 and v : R 2 —> R 2 are continuous functions. 

We establish the existence of at-least one solution of (4.2.1)-(4.2.2) with the help of LUSs method 
in the presence of the existence of well ordered LUSs such that (ai,a 2 ) ^ (/3i, f3 2 ) on [0,1] x [0,1]. 

The monotonicity assumptions on <p and ip as stated in Theorem 4.1.1 does not address to the 
existence criterion for a very important BVPs like periodic and anti-periodic. But in this study, 
we introduce a new generalized BCs 4.2.2 that satisfy the existence criterion for a periodic and 
anti-periodic BVPs as a particular case under some nronotonicity assumptions on the arguments of 
the boundary functions /r and v as discussed in detail in the introductory chapter of this thesis. 

This section deals to the concept of CLUSs that unifies the treatment of many linear and nonlinear 
BVPs under some monotonicity assumptions on the arguments of the boundary functions /i and v. 

4.2.1 Coupled lower and upper solutions 

Definition 4.2.1. We say that a couple of functions (ai,a 2 ), (/3 i,/? 2 ) G C 2 [0,1] x C 2 [0,1] are 
coupled lower and upper solutions for the coupled BVP (4.2.1) and (4.2.2), if (ai,a 2 ) is a lower 
solution and (/?i,/ 3 2 ) is an upper solution for the system (4.2.1),(on, a 2 ) -< (/3i,/? 2 ), if op < /3i and 
a 2 < /3 2 , and 

M/?i(o),/3 2 (o),/3;(o),/3;(o),/3;(i),/?;(i)) r (o,o) r mm o), « 2 (o), ^(o), «;(o), «;(i), «;(i)), 

MA(0),/3 2 (0),/?;(0),/3;(0),a' 1 (l), a ;(l)) ^ (0,0) ^ / z(« 1 (0), CK2 (0),« , 1 (0),«;(0),/3;(l),/3;(l)), 

(ar(l), o 2 (l)) + i/ (M0), & (0)) = (0,0), 

(/3 1 (1),/3 2 (1)) + i/(ai(0),a 2 (0)) = (0,0), 

(a 1 (l),a 2 (l)) + ^(a 1 (0),a 2 (0)) = (0,0), 

(/?!(!),/3 2 (1)) (ft (0),/3 2 (0)) = (0,0). (4.2.3) 
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4.2.2 Existence of at least one solution 

Theorem 4.2.1. Suppose the nonlinear coupled BVP (4.2.1)-(4.2.2) having ( 011 , 0 : 2 ) and (/3i,/3 2 ) 
coupled lower and upper solutions respectively and the nonlinear functions fi and f 2 relative to 
the intervals [a\(t), f3i{t)\ and [ 02 (f), (3 2 {t)\ satisfy the Nagumo assumptions respectively. Moreover, 
assume that p is nondecreasing in the third and fourth arguments. In addition suppose that the 
function v in [oi(0),/3i(0)] x [ 02 ( 0 ),/A (0)] is monotone and the functions 

h(otl ,a 2 ){x,y ) := p (o^O),o 2 (0),o 1 (0), 02 ( 0 ),a ;,yj , 

M(/ 3 i ,02)( x ,y) '■= h (/?i(0),^ 2 ( 0 ),/3i(0),/? 2 (0), x,y'j , 

have got the same kind of monotonicity as u, Then the nonlinear coupled system of BVPs (4.2.1)- 
(4.2.2) have at least one solution (ui,u 2 ) G [oi,/3i] x [ 02 ,^ 2 ] an d (— Ni,—N 2 ) S (u 1 (t),u 2 {t)) S 
(N u N 2 ),te [0,1]. 

Proof. Let A > 0, and consider the modified system of BVP 

-«i(t) + Xui(t) = Ff (t,u 1 (t),u 2 (t),u 1 (t),u 2 (t)^ , t G [0,1], 

-u 2 {t) + A u 2 {t) = F 2 (t,u 1 (t),u 2 (t),u 1 (t),u 2 (tfj , t G [0,1], 

p* ^Wi(O), w 2 (0), Wi(l), w 2 (l),«i(0),*4(0)) = (ui(0),u 2 (0)), 

( Ul (l),u 2 (l)) + v* (u 1 (0), u 2 (0)) = (0,0), (4.2.4) 


with 

F* (t,ui(t),U2(t),r(t),s(t)) = A (t 1 A 1 {t,u 1 {t)),A 2 (t,u 2 (t)),T 1 (r{t)),T 2 (s(t))) + AA 1 (i,rt 1 (t)), 
F 2 {t,ui{t),u 2 {t),r{t),s{t)) = f 2 (t, Ai(t,tti(t)), A 2 (t,u 2 (t)),T 1 (r(t)),T 2 (s(t))) + AA 2 (t,u 2 (t)), 


and 


p* ( j,k,l,m,n,o ) 
v* (j, k) 

A i(t,z) = 

Uz) = 


= (Ai(0, j),A 2 (0, k)) + p(j,k,l,m,n,o ), 
= v (Ai(0, j), A 2 (0, k )), 

Pi{t) if z£Pi(t), i — 1,2, 

-z if oa(t) < z < Pi(t), i = 1, 2, 
oti(t) if z ^ Oj(t), 2 = 1,2, 

N., if z£N it 2 = 1,2, 
z if — IVj < z < Ni , 2 = 1,2, 

-Ni if z^-Ni, 2 = 1,2. 


Note that if (ui,u 2 ) G [oi,/3i] x [ 02 , fi 2 \ is a solution of (4.2.4), then (u\,u 2 ) is a solution of 
(4.2.1)-(4.2.2). 

For the sake of simplicity we divide the proof in four steps: 

Step 1: We define the mappings 


L, N : C 2 [ 0,1] x C 2 [ 0,1] Cg[ 0,1] x C o 2 [0, l]xl 2 x R 2 , 
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by 


[L(u 1 ,u 2 )](t) 



'Ui(O) —A / ui(s) ds, u 2 (t) 
Jo 


w 2 (0) — A / 1 ^ 2 ( 5 ) ds, 
Jo 


(ui(0),u 2 (0)),(ui(l),u 2 (l)) 


and 

[N(u 1 ,u 2 )](t) 


F*(s, tii(s), u 2 (s), %(«), u 2 {s))ds, / F 2 * (s, ui(s), u 2 (s), ti^s), u 2 (s))ds, 
3 Jo 

/z*(wi(0), u 2 (0), ui(l), w 2 (l), «i(0), w 2 (0)), 

-i/*(ui(0),u 2 (0))Y 


Since i 7, ^(s,'iii(s),U 2 (s),u 1 (s),u 2 (s)) and F 2 (s,ui(s),u 2 (s),u 1 (s), w 2 (s)) are bounded on [0,1] x ]R 4 
and integral is a continuous function on [0,1]. Further p* and v* being constant functions are 
continuous. Therefore [iV(wi,u 2 )] is continuous on [0,1]. Further, the class {N(ui,u 2 ) : ui,u 2 G 
C 2 [0,1]} is uniformly bounded and equicontinuous. Therefore in view of Arzela-Ascoli theorem 
{N(m,u 2 ) : ui,u 2 G C 2 [ 0,1]} is relatively compact. Consequently TV is a compact map. Also from 
Lemma 3.1.1 with a = 1,6 = 0, c = l,d = 0 and E = 0, F = 1, G = 0, H = 1, L _1 , exists and is 
continuous. 

On the other hand, solving (4.2.4) is equivalent to find a fixed point of 
L~ 1 N : C 2 [0,1] x C 2 [0,1] C 2 [0,1] x C 2 [0,1], 

Now, Schauder’s fixed point theorem guarantees the existence of at least a fixed point since L~ 1 N 
is continuous and compact. 

Step 2 : If ( u\,u 2 ) is a solution of (4.2.4), then (ui,u 2 ) € [ai,/3i] x [a 2 ,/3 2 ]. By definition of /x* , 
we see that (ui(0),u 2 (0)) G [ai(0),/?i(0)] x [a 2 (0),/3 2 (0)]. Thus, if v is nondecreasing, we have by 
condition (4.2.3) 

(ai(l),a 2 (l)) = —v (/?i(0), /? 2 (0)) ■< —v (ui(0), u 2 (0)) = (ui(l), u 2 (l)) ^ -v («i(0), a 2 (0)) 

(o 1 (l),a 2 (l)) =< («i(l),U2(l)) d (/3 i(0),/? 2 (0)), (4.2.5) 

similarly, if v is non-increasing, then (4.2.5) holds. Hence (wi(l), u 2 (l)) G [oq(l), /?i(l)] x [a 2 (l), /3 2 (1)]. 
Now, it remains to show that (ui,u 2 ) G [ai,/3i] x [a 2 ,/3 2 ] for t G (0,1). 

We claim (ui,u 2 ) A (Pi,p 2 ). If (ui,u 2 ) ^ (Pi,p2), then either u\ ^ Pi and/or u 2 ^ p 2 . If 
Mi ^ Pi, then there exists some to G [0,1] such that ui(t-o) — /?i(to) > 0. So, u\ — pi attains a 
positive maximum at to G [0,1]. Thus (m — Pi) (to) = 0 and (ui — 0 i) (to) < 0. But, 

(«1 - Pi)" (to) > -Ff Yo,Ml(fo),M 2 (to),M' 1 (to),M 2 (fo)) + Xui(t 0 ) + fl (to, Pi (to), &2 (t 0 ), Pi (fo), P 2 (fo)) 

= -fl (to,Pl(to),p2(to),U 1 (to),U 2 (to)^J - XPi (to) + Xui(to) 

+ fl (to, Pl(to), P2O0), Pl(to), P 2 (to)^ 

= — fl (to, Pl(to), p2(to), Pl(to), P 2 (to)^j — XPi(to) + Xui(to) 

+ fl (to, Pl(to), P2O0), Pl(to), P 2 (to)^ 

= X (til (to) — Pi(to)) > 0, 
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a contradiction. Similarly one can show that (<* 1 , 0 : 2 ) ^ (ui,u 2 ). Hence (ui,it 2 ) G [oi,/0i] x [a 2 ,/3 2 ]. 
Step 3: If (u\,U 2 ) is a solution of (4.2.4) then (ui,U 2 ) satisfies (4.2.2). 

We claim 

(o 1 (0),o 2 (0)) ^ (u 1 (0),u 2 {0))+p (u 1 (0),u 2 (0),u 1 (0),u 2 (0),ii 1 (l),u' 2 (l)^ ^ (ft(0),/3 2 (0)). (4.2.6) 
If (oi(0),o 2 (0)) ^ («i(0),u 2 (0)) + p (« 1 (0),u 2 (0),ui(0),U2(0),w , i(l))%(!)) 1 Bren 
(ui(0),u 2 (0)) = p* (u 1 { 0 ),u 2 ( 0 ),u 1 ( 0 ),u' 2 ( 0 ),ii 1 (l),u 2 (l)^ 

= (Ai(0, «i(0)), A 2 (0, w 2 (0))) + p(ui(0),u 2 (0), Wi(0), w 2 (0), Wi(l), u^l))) 

= (o 1 (0),o 2 (0)). (4.2.7) 

If v is non-decreasing then we have 

(ui(l), w 2 (f)) = -^(mi(0),u 2 (0)) = —v (oi(0), o 2 (0)) = (/3i(l),/3 2 (l)). (4.2.8) 

Using (4.2.7),(4.2.8) and Step 2, we have ^u 1 (0),w 2 (0)^ y ^aq(O), a 2 (0)^ and ^(1), w 2 (l)^ y 
(/^C 1 ),^ 1 ))- But 

(u 1 (0),u 2 (0)) + p (u 1 (0),u 2 (0),m 1 (0),u 2 (0),u 1 (1),u 2 (1)^ 

= (oi(0),o 2 (0)) + p ^oi(0),o 2 (0),'u 1 (0),m 2 (0),m 1 (1),m 2 (1)^ 

h (ai(0),o 2 (0)) + p (o 1 (0),a 2 (0),o' 1 (0),o 2 (0),u' 1 (l),u 2 ( 1 )) 

= (ai(0),a 2 (0)) + p { Qii,ac2) (wi(l),u 2 (l)) 

h (oi( 0 ),a 2 ( 0 )) + p ^o 1 ( 0 ),a 2 ( 0 ),a 1 ( 0 ),a 2 ( 0 ),^ 1 (l),^ 2 (l)) 

h (ai(0),a 2 (0)), (4.2.9) 

a contradiction. Similarly if v is non-increasing we get same contradiction. Consequently, (4.2.6) 
holds. By definition of v* and Step 2, the second boundary condition is obvious. Consequently 
(wi,w 2 ) satisfies (4.2.2). 

Now the rest of the proof is analogous to the proof of the Theorem 4.1.1. □ 

4.2.3 Application of the theoretical results 

In this section we check the validity of the Theorem 4.2.1 by considering an example. 

Example 4.2.1. Consider the following nonlinear coupled system of BVP 

-«i (t) = -5 (u^t) + u 2 (t )) - (ui(t) + l) 2 - (u 2 (t) + l) 2 + sin 2 (nt), t € [ 0 , 1 ], 

— u 2 (t) = —10 (uxit) + 2 u 2 (t)^ — (m(£) + 2) 4 — (u 2 (f) + l) 4 + sin 2 (7rt), t € [ 0 , 1 ], ( 4 . 2 . 10 ) 

corresponding to nonlinear CBCs 

^Mi(0)ui(0) - m 2 (0)m 2 (0),mi(0)mi(1) - m 2 (0)w 2 (1)) = (0,0), 

(«i(0)tt 2 (0) + mi(0)m 2 (1),mi(1)m 2 (0) + ui(0 )m 2 (1)) = (0,0). 


(4.2.11) 
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The nonlinearities fi(t,xo,xi,X 2 ,X 3 ) = —5 (*2 + * 3 ) — (xo + l) 2 — (xi + l) 2 + sin 2 (7rf) and 
f 2 (t, * 0 , * 1 , X 2 , * 3 ) = —10 (x 2 + 2a; 3 ) — (x 0 + 2) 4 — (xi + l ) 4 + sin 2 (7rt) are continuous on [0,1] xl 4 . Then 
the functions ( 01 , 0 : 2 ), (/3i,/?2) : [0,1] x [0,1] ->1 xR defined by 
oi (t) = —t 2 — t, 02 (t) = —t and fii(t) = t 2 + t,p 2 {t) = t 
are, respectively, lower and upper solutions of the coupled nonlinear system (4.2.10), according to the Defi¬ 
nition f.1.1. 

Moreover, it can be easily check that the functions oi (t) = —t 2 —t, 02 (t) = —t and fii(t) = t 2 +t,fi2(t) = t 
satisfy (4.2.3). Consequently, these functions are, respectively, coupled lower and upper solutions of the 
coupled BVP (4.2.10)-(4.2.11), according to the Definition 4 - 2 . 1 . 

Also notice that the functions fi(t, xo, xi,X 2 , * 3 ) = —5 (*2 + * 3 ) — (xo + l) 2 — (*1 + l) 2 + sin 2 (7rt) and 
f 2 (t, xo, xi, * 2 , * 3 ) = —10 (X 2 + 2 x 3 ) — (x 0 + 2) 4 — (x! + l ) 4 + sin 2 (7rt) satisfies the Nagumo condition (4.1.4) 
with g(x 2 + X3) = —5(*2 + * 3 ) — 1, 12,13 and g(x 2 + * 3 ) = — 10 ( 2:2 + * 3 ), * 2 , *3 € R, respectively. 

Hence, by Theorem 4 -%-T there is at least one solution (ui(t),U2(t)) of the coupled BVP (4.2.10)-(4.2.11) 
such that, for all t £ [0,1], 

(— t 2 - t, — t) ■< (ui(t),U2(t)) ■< (t 2 + t,t) . 



Chapter 5 

Coupled Systems Of Multi-terms 
Fractional Order Partial 
Differential Equations 


Partial differential equations of fractional order have rich applications in almost every field of science 
based including, but not limited to, fluid dynamics, physics, chemistry, aerodynamics, signal and 
image processing, chemical engineering, economics, and even in psychology, [9, 81, 84, 101, 142, 161, 
176, 177], 

The applications of coupled systems of fractional order partial differential equations (FOPDEs) 
have been widely discussed in the existing literature. For example in Biomechanics, coupled systems 
are used to model the electrical activity of the heart [168, 169], in chemical and material engineering 
they are used to model the systems consisting of a plug flow reactor and a continuous stirred tank 
reactor [10, 135], in solid mechanics, the physical phenomenon used to discuss the dynamics of 
multi-deformable bodies are modeled using coupled systems of FOPDEs [112, 152]. 

Based on a preliminary mathematical modeling, determination of exact and/or approximate 
solutions of fractional order differential equations has proved to be a major emerging research area 
that has captured the interest of the researchers round the world [35, 36, 58, 92, 93, 94, 100, 127, 136, 
155, 187]. For instance, many scientists have devoted themselves to establishing efficient and reliable 
numerical schemes for finding the approximate solutions of FOPDEs using orthogonal polynomials 
which reduce them into an easily tractable system of algebraic equations. 
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Certain types of FOPDES where the non-integer order is left arbitrary, have been solved numer¬ 
ically applying such ideas and types of orthogonal polynomials, [23, 51, 113, 114, 115]. Recently, 
particular types of FOPDES have been solved by constructing new operational matrices of frac¬ 
tional derivatives and integrals based on such schemes using both orthogonal and non-orthogonal 
polynomials, [57, 115, 155, 187]. In [113, 114], the numerical schemes are developed with the help 
of operational matrices based on orthogonal Legendre polynomials. 

Based on above works, we opted to develop an efficient numerical scheme able to treat a gener¬ 
alized class of coupled systems of FOPDEs, by developing the operational matrices for a particular 
kind of polynomials, namely, the Jacobi polynomials, (JPs). Unlike in the case of the Legendre poly¬ 
nomials, the JPs have the ability to approximate the solutions with the aid of two parameters, which 
is indeed a more generalized way of approximation, since in comparison to the works, [113, 114], 
herein, the operational matrices developed can deal with the FOPDEs having mixed partial deriva¬ 
tives of fractional order. Moreover, in [115], the operational matrices has been developed using one 
dimensional JPs. In our case, operational matrices are exhibited for two dimensional JPs. 

This chapter is organized as: In Section 5.1.1, the two dimensional orthogonal shifted JPs are 
presented along-with their orthogonality expressions. In Section 5.1.2, we develop the analytical 
relation to determine the error bound for a sufficiently smooth function. In Section 5.1.3, the 
operational matrices of fractional order integrals and derivatives are established. In Section 5.2, 
the applications of the developed matrices of fractional order integrals and derivatives are discussed 
by developing the new schemes for the multi-terms FOPDEs and its coupled systems. In Section 
5.2.3, the applicability of the developed method is checked by considering some test examples. In 
Section 5.3.1, the notion of three-dimensional orthogonal shifted JPs is discussed along-with their 
orthogonality relationship. In Section 5.3.3, the analytical relation is developed to determine the 
error bound for a sufficiently smooth function. In Section 5.3.4, the applications of fractional order 
integrals and derivatives is discussed by developing the numerical scheme for finding the approximate 
solutions of high dimensional multi-terms FOPDEs. In Section 5.3.5, some illustrative examples are 
tested to check the reliability and efficiency of the proposed method. 
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5.1 Multi-term fractional partial differential equations mixed 
type Coupled Systems 

In this section, we are interested in developing the numerical scheme based on operational matrices 
of fractional order integrals and derivatives of the following multi-term FOPDEs and its coupled 
systems of the type 


d ai U(x,y) d 11 U{x,y) d ei U(x,y) 
d . =ai dyn + ° 2 dx<*/ 2 dyei/ 2 

corresponding to the following initial conditions with j = 0,1, ...,n 


Gi(x,y), 


U^(0,y) = Hj(y), 


and a coupled system of a generalized class of FOPDEs 


d ai U(x,y) d 11 U{x,y) d l2 V(x,y) d l3 V{x,y) 


dx ai 


=a i 


a 2 - 


dx T' 2 


a 3 - 


dy^ 


d 3i U{x,y) d e3 V(x,y) 

4 dx ei / 2 dye ^/ 2 5 dx e 3 / 2 dy ^ 3 / 2 1 ’ ’ 

d° 2 V(x,y) d pl V{x,y) d P2 U(x,y) d p3 U(x,y) 
=0i -r o 2 —-h o 3 - 


dx a2 


dxP 2 


dy p 1 


dypi 

d 33 U(x,y) d 32 V(x,y) 

4 dx ei / 2 dy ei / 2 5 dx e2 -l 2 Qyeil 2 2 ’ ’ 


corresponding to the following initial conditions with j = 0,1, ...,n 


UW(0,y) = Hj(y), y«(0 ,y) = G^y), 


(5.1.1) 


(5.1.2) 


(5.1.3) 


(5.1.4) 


where n < a\, er 2 < n+1, j = 1,..., 5 are all real constants and F\{x, y), E 2 (a;, y), U(x, y), V (x, y), € 

C([0, A] x [0, A]). 

We recall the concept of two-dimensional shifted JPs in two parameters. 

5.1.1 Two parametric shifted Jacobi polynomials in two dimensions 

This section elaborates the idea of orthogonal shifted JPs in two dimensions along-with its or¬ 
thogonality relationship with the aid of a weight function. The following relation defines the one¬ 
dimensional shifted orthogonal JPs on the interval [0, A] as, (see [115] ) 

0 (Tl,T2) (x) = V + T 2 + 1 ) r (i +1 + n + r 2 + 1) , 

A ’ j h ra + r 2 + l)r(j+T 1 +r 2 + l)(j-0!/!A' ’ 3 ’’ 


(5.1.5) 
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The orthogonality expression is 

[ Q { L)f 2 \ X )Q { L]r\ X ) R A UT2 \ X ) dx = W A4 T2)S 3,U 
Jo 

where weight function R^ 1,T2 ^ is defined as follows 


R 


(ri ) / \ _ 


(a:) = (A - x) 7 


and 


A^+ir(i + r 1 + i)r(i + T 2 + i) 

A ’ 4 (2i + n + 72 + l)T(i + l)T(i + n + 72 + l) ' 


(5.1.6) 


(5.1.7) 


(5.1.8) 


It means that any integrable function r{x ) defined on [0, A] can be approximated in terms of the 
series of orthogonal shifted JPs as follows 


r{x) « J2 D bQ { A,b 0*0 


(5.1.9) 


6=0 


The series approximation approaches to the exact function as n —► oo. With the help of (5.1.6), 

(5.1.7) and (5.1.8), the coefficient D b can be easily computed. The vector expression of (5.1.9) is as 
following 

r(x) « N (x), (5.1.10) 

where < hjv(2 ; ) represents the N terms vector function with N = n+1, and G indicates the coefficient 
vector. On the same fashion the two-dimension JPs of order N defined on [0, A] x [0, A] can be 
described using the following relation 

Q ( L]Z 2 \ x ^y) = ( ( ?A, 1 i r2) ( a; ))( ( 5A, 1 i’ r2) (y)) 5 m = Nj+i+ 1, j = 0,1,2,...,n, i = 0,1,2, ...,n. (5.1.11) 

The orthogonality expression of Q^m\ x ^v) I s as following 

[ A / A (0i T r 2) W)(0iT 2) (2/))(OiT 2) (*))(Qiy 2) (y)) R( A’ T2 \ x ) R( A 1,T2) (y)dxdy = fPi7/ 2) ^, d wi i ; i c ’ r2) ( 5 b , c 
Jo Jo 

(5.1.12) 

So, the function r(x, y) defined and square integrable in the region [0, A] x [0, A] can be expressed 
in terms of the series of two-dimensional shifted JPs as follows 

n n 

r{x, y) ~^^ D ba{Q { l]i T2 \x)){Q { £a 2 \y))’ 

b —0 a—0 


(5.1.13) 



Db a is computed using the following relation 


pA p A 


Dba = 




r(x,y)(Q i £f 2 \x))(Q i ^ T2 \y))R% uT2 \x,y)dxdy. 


(5.1.14) 


The weight function R < ^ 1,T2 \x,y) is given by 

R^ T2 \x,y) = R^\x)R^\y). 


(5.1.15) 


(5.1.13) can also be expressed in the form of vector notation as 

N 2 

r(x, y) « Y D mQt£\ x ,v) = G^ 2 ^ N 2 {x,y), with m = Nb + a + 1, and D m = D ba , (5.1.16) 

m—1 

where Gn 2 refers to coefficient column vector of N 2 x 1 dimensions and 4 >jv2 (x,y) indicates the 
column vector of dimensions N 2 x 1, defined as 


&N*(x,y) = ( 4>u(x,y) ••• 4>iN(x,y) (fai (x,y) ■■■ 4>2N(x,y) ••• ^nn{x,v) ) 


(5.1.17) 


where (j) j+ + i+1 (x,y) = (Q^ r2) (x))(Q^l’ T2j (y)), j,i = 0,1,2, ...,n. 


(ti,T 2 ) I 


5.1.2 Error Analysis 


In this section we develop an analytical relationship for the error of approximation for sufficiently 
smooth function r(x, y) € C([0, A] x [0, A]). For this assume that N (x, y ) the space of N terms 
JPs and g(N,N)( x :V) is its best approximation in YI(n n)( x ^v)- Th en according to the definition of 
the best approximation for any polynomial f(N,N)( x , y) of degree < N in variable x and y, we can 
write 

\\r(x,y) - r {N}N) (x,y)\\ 2 < \\r{x, y) - f( N>N )(x, y)\\ 2 . (5.1.18) 


Then by following the same procedure as adopted in [82], we can find the bounds for the errors given 
below 


\r(x,y) - f(N,N)( x ,y )I < ^ 


1 / (A — 0) 


N +1 


1 ^ (A — 0) \ 
N J 


N+l 


max (x,y)\ 


N 

d N+1 r(x, y) 


max 


d N+1 r(x,y) 
dx N+1 


16 


1 / (A — 0) \ iV+1 / (A — 0) \ 


d y N + l 

N+l 


Q2N+2 


N ) V N ) ™ a j)'dx N + 1 dy N + 1 


r(x,y) |. 


(5.1.19) 
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T i \ I d N + 1 r(x,y) | \ 

Let Ai = max {XtV) \ dx it + i y, \, A 2 = 
Then, we have the following relation 


fTlCLX^x ,y) 


d N+1 r(x,y ) 
dy N + l 


and A 3 


max (x,y) 1 dx N+iQyN+i r {x, y) 1 


kOuy) - f(N,N)( x J y)\ < \ ( ^ A n Ai 


Finally using (5.1.18), we have the following relation 

\\r(x,y) - r {N>N) (x,y )|| 2 <$/(A - 0 )(A - 0)(^ A, 

(5.1.21) 


For the detail study of the above analytical findings, we refer the reader to see [141]. 

The results of the next section are very important to develop the numerical scheme. Our main 
purpose in this section is the development of the new operational matrices of integrals and derivatives 
using two parametric orthogonal JPs. Using these operational matrices the fractional order problems 
(5.1.1)-(5.1.4) are reduced into a problems of solving a system of algebraic equations of Sylvester 
type which are simple in handling and can be solved by any computational software. 


5.1.3 Operational matrices of fractional order integrals and derivatives 

Lemma 5 . 1 . 1 . Let 4>jv 2 (x,y) be the function vector as defined in (5.1.17), then the fractional 
integration of order p of 4*^2 ( x, y) with respect to x is given by 

I£{* N >(x,y)) ~ Q p ^ xN Am(x,y), (5.1.22) 

where Q% 2 xN 2 indicates the operational matrix of integration of order p and is defined as 




^1,2 ,1 

^1 ,q,l 

^1 ,N 2 ,l \ 


^ 2,1 ,1 

^2,2 ,1 

^2 ,q,l 

^ 2 , at 2 ,; 

^N 2 xN 2 ~ 


^ r,2,l 

^ r,q,l 

^r,N 2 ,l 


\ ^N 2 ,l,l 

^N 2 ,2,l ’ ’ ’ 

^ N 2 ,q,l 

^N 2 ,N 2 ,l ) 

Nj + i + 1, v — Nb +1 + 1, ^ r , 

Q 1 1 — 

: for i,j, a,b = 

0 ,1, 2,..., m 


b 

~\j,i,b,a,l = ^2 1+J dj tit a, 

1=0 b,l,p 
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, _ <■ (-i) J fc ( 2 i + n + t 2 + i)r (j + i)tq + fc + n + t 2 + i)r(z + p + k + r 2 + i)r(ri + i)A p 

jia l ' a r(j + n + i)r(fc + r 2 + i)(j — fc) !/c!r(Z + p + fc + r 2 + n + 2) 


(-l) b -'r(6 + r 2 + 1)T(6 + l + n + r 2 + l)r(l + Z) 
T(Z + r 2 + 1)T(Z> + t 1 +t 2 + 1 )(b - Z)!Z!r(l + Z + p)A*' 


(5.1.23) 


Proof. To prove the above result take the fractional integration of (5.1.11) of order p with respect 
to x, we have 


ip Qa% 2 \ x ^v) = 


(~l) fc l T(b + r 2 + 1)T(5 + l + Ti + r 2 + 1) ;„(ti,t 2 )/ x 

T(Z + r 2 + 1)T(6 + n + r 2 + 1)(6 — l)\l\A l x 


Applying fractional integral, above relation becomes 

rp 0 (Ti,r 2 )/ x = (^l) b - f r(b + r 2 + l)T(fe + Z + Ti + r 2 + 1)T(1 + Z) j +(1 ( n ,r 2 ) , x 
x^A,m 1 . y> ^r(Z + r 2 + l)r(fo + r 1 +T 2 + l)(6-Z)!Z!r(l + Z + p)A i A ’° 

(5.1.24) 

The expression x k+v Q^^ 2 \y) in view of two-dimensional JPs becomes 


xl+p Q [ A,a 2 \y )~ ^J2 d iP a Q ( Af 2) ^QaT 2) 

3 =o *=0 


(5.1.25) 


where cZ i>ija = (T1 . T2) 1 , T1 . T3 , f 0 f 0 x l +PQ^ T2 \y)Q^ 1 f 2> (x)Q^ 1 J T2> (y)R^ 1 ’ T2> {x,y)dxdy. And in 

VV A ,j i vv A ,i 

the light of orthogonality condition, it may be written as 


di i a — 


W^f 2 > Jo 


x 1+p QaT ) ,T2 \ x )dx. 


(5.1.26) 


With the help of (5.1.5) and (5.1.7), (5.1.26) can also be written as 


dpi a — r. 

UA 


p— V (^ 1 ) J kT (j + T 2 + 1 ) r (j + fc + Ti + 7-2 + 1) r A /+p+fc+T2 , _ xriw x 

W 2 ) “ r(fc + t 2 + l)T(j + n +7-2 + 1)0 - fc)!fc!A fe Jo 


"Aj fe=0 v ~ 7 w ‘ 7 

(5.1.27) 

Applying convolution theorem of Laplace transformation to the integrand of the above relation 
and then taking inverse Laplace of the obtained relation, we can write the generalized form of the 
coefficient involved in (5.1.25) as 

6p a (-1 y~ k T(j + r 2 + i)r(j + k + n + t 2 + i)r(z + P + fc + r 2 + i)r(Tj + i)A d+P+^+^+i) 
Jia u/( ri > T2 ) T(k + To + HIT? + ti + T9 + 111? — fcllfc!A fe rfZ + o + fc + T9 + ti +21 


wt-' n) kf~o r(fc + r 2 + l)T(j + n + r 2 + l)(j — fc)!fc! A fe T(Z + p + k + r 2 + n + 2) 

(5.1.28) 

Using(5.1.8) in (5.1.28) the most generalized form of the coefficient is as following 


djia — ^ ' 


(—1 y k (2j + n + r 2 + l)T(j + l)T(j + Zc + n + r 2 + 1)T(Z + p + k + r 2 + 1 )T(ti + 1)A P 
T(j + Ti + l)T(fc + r 2 + l)(j - fc)!fc!T(Z + p + fc + r 2 + ly + 2) 


(5.1.29) 
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Now, for convenience in handling the notations, let 

II {-i) b ~ l r(b + T 2 + i)r(& + l + n+T 2 + i)r(i + z) 

U - r(Z + r 2 + l)r(6 + t 1 +t 2 + 1)(6 - Z)!Z!r(l + l + p) A*' 

In the light of (5.1.25), (5.1.29) and (5.1.57), (5.1.24) can be written as 


IP xQ { A,m\xiy) = |+) 




( r i ,t 2 ) / 
A,i V 


x )QAj T2 \y)- 


l —0 6,/,p j=0 i—0 


(5.1.30) 


(5.1.31) 


Let 

b 



(5.1.32) 

/—0 b,l,p 


we get 

71 71 


W) = EE^w^ T2) w^ T!) (^ 

j— 0i —0 

(5.1.33) 

Using the notations q = Nj + i + 1, r = AZ> + a + 1, = Ij.i.a.&.z for i, j, 

a, b = 0,1,2,3,..., n, 

we get the desired result. 

□ 


Lemma 5.1.2. Let dE> ^2 (a:,?/) be the function vector as defined in (5.1.17), then the derivative of 
order r of & N 2 ( x , y) with respect to x is as follows 


Dx{&n 2 (x,y)) A! D T f} 2 xN 2 ^N 2 {x,y), 


(5.1.34) 


where C^f xAr2 


indicates the operational matrix of derivative of order r and is defined as 



/ ri 1;1)l 

^ 1 , 2 ,; 

^1 ,q,l 

^l,N 2 ,l \ 



^2,2,Z 

^2 ,q,l 

^2,N 2 ,l 

f) T > x _ 

^N 2 xN 2 ~ 

^r, 1 ,1 

^r,2,Z 

^ r,q,l 

^r,N 2 ,l 


\ ^V 2 ,1,Z 

^N 2 ,2,1 ■■■ 

^N 2 ,q,l ’ ' ' 

' ^N 2 ,N 2 ,l J 

Nj + i + 1, c[ — A/"6 + c + 1, 

\q,l — ~~\j,i,b,a,l 

for i,j, a,b = 

0,1,2,..., n, 


*~\j,i,b,a,l — ^ ^ (+J dj,i, a , 

Mr | b,l,r 


(5.1.35) 


(5.1.36) 


djia 


And 


Si,a ^2 


(— 1 y k (2j + n + t -2 + i)T{j + i)r(j + k + ti + t - 2 + i)r(z — t + k + r 2 + i)r(n + i)A T 


k =0 


r (j + ti + i)r(fc + t 2 + i)(j — fc)!fc!r(z — r 


(5.1.37) 


w 

b,l,r 


(~i) b ~ l r(b + t 2 + i)r (6 +1 + n + r 2 + i)r(i + z) 
r(Z + r 2 + l)r(6 + Tx + r 2 + 1)(6 - Z)!Z!L(1 + Z - r)A ; ' 


(5.1.38) 
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Proof. To prove the above result take the fractional derivative of (5.1.11) of order r with respect to 
x, we have 

D tq(ti,t 2 ), x == y- (-l) b ~ l T(b + r 2 + 1 )T(b + l + n + r 2 + 1) r , (ri , T2 ), , 

X ^A,m l*.!/! A r(« + T 2 + l)r(6 + T 1 +r 2 + l)(6-0!^A i x ^ A ’° l2/j ’ 

In the light of the definition of the fractional derivative, above relation takes the following form 

nr 0 (ri,T 2 ) ( x = y-' {-l) b ~ l T(b + t 2 + l)T(b + l + n + r 2 + 1)T(1 + l) 1 -tq(ti,t 2 ), x 

iVa ’™ 1 ,2/j ^ r(; + T2 + i)r(6 + T 1 + T 2 + i)(6-0^!r(i + ^-r)A ( ^ ^ A ’ a 


The expression x l r Q^' a T2 * (y) in view of two-dimensional JPs becomes 


x'~ T Q { L]a 2) (y) « EEdM.OkT ) (*)0& T,,) (»)> 

j—02=0 

r A r A j-rn( T i’ T 2)/,„\n( Tl ' T! )/x.\n( ri ' r! )^ l 'i p(ri,T 2 ) 


(5.1.39) 


(5.1.40) 


where dj-, <ia = (T1 , T2 > f 0 f 0 x l T Q { ^ 2 \y)Q^)f 2 \x)Q^: T2) {y)R^' T2 > {x,y)dxdy. 

VV A ,j i vv A ,i 

And in the light of orthogonality condition, it may be written as 


dj,i, a = } f x 1 T Q^ 1 j T2) (x)R i l UT2 \x)dx. 

’ Jo 

With the help of (5.1.5) and (5.1.7), (5.1.41) can be written as 


(5.1.41) 




by_ (-l)J k T(j + r 2 + l)T(j + k + Tl + r 2 + 1) /~ A , T+fc+T2 t . 

^ttr(k + T 2 + mj + T 1 + T 2 + i)(j-k)\k\A^J 0 x [ ^ x) ■ 


A ,j k— 0 


Applying convolution theorem of Laplace transformation to the integrand of the above relation 
and then taking inverse Laplace of the obtained relation, we can write the generalized form of the 
coefficient involved in (5.1.40) as 


djia — ^ ' 


(—1 y k {2j + n + t -2 + l)T(j + l)T(j + k + ti + t -2 + 1)T(Z — r + k + r 2 + l)L(n + 1)A T 


T (j + n + 1 )T(fc + r 2 + l)(j — k)\k\T(l — r + k + r 2 + n + 2) 


Now, for convenience in handling the notations, let 

■ | (—l) h_i T(& + r 2 + 1)T(6 + 1 + t 1 +t 2 + i)r(i + l) 

“ T(l + t 2 + 1 )T(6 + n + r 2 + 1)(6 - Z)!Z!T(1 + l - t)A 1 ' 
In the light of (5.1.40), (5.1.43),(5.1.44), (5.1.39) takes the form 

b n n 

D lQ ( L% 2 \x,y) = Y U YY d ^ i ,aQ { ^l T 2 \x)Q { ^i T 2 \y)- 

/=[>] b,l,T j= 0 i— 0 


(5.1.43) 


(5.1.44) 


(5.1.45) 


— y ~ |dj 

i=frl b ,1 ,t 


(5.1.46) 
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we get 


^iE T2 W) = EE^> 6 

j= 0 i—O 


,iQaT 2) ( X )Q { L]i T 2 \y)- 


Using the notations q = Nj + i + 1, r = Nb + a + 1, ^ r ,q,i = for L j , a, b 

we get the desired result. 


(5.1.47) 

0,1,2,3,..., n, 
□ 


Lemma 5.1.3. ZeZ $^2 (a:, y) Z>e ZZie function vector as defined in (5.1.17), ZZien f/ie derivative of 
order r of $^2 (x, y ) with respect to y is as follows 


D T y {$ N 2 {x,y)) ~ D^ 2 xN 2 &n 2 (x, y), 


(5.1.48) 


where D T jf% N2 indicates the operational matrix of derivative of order r and is defined as 



^1.2 ,1 

^1 ,q,l 

^l,N 2 ,l \ 


^2,1 ,1 

^2,2,1 

^2 ,q,l 

^2,N 2 ,l 


^r,"l . ! 

~^r,2,l 

^ f,q,L 

^ r,N 2 ,l 

, (5.1.49) 

\ ^iV 2 ,l,Z 

^N 2 ,2,1 ’ 1 

"^N 2 ,q,l ’ 

’ ’ ^ N 2 ,N 2 ,l ) 



and r = Nj + i + 1, g = ZVZ> + c + 1, D r , 9) j = Wj ti ,b, c ,i for j, i,b,c = 0,1, 2,..., n, 

b 

^ j ,i,b,c,l — ^ ^ ^b,l : r\j,i,ct (5.1.50) 

i=M 


Xj'ic — < 5 . 7,6 E 


(_l)j-m( 2 j + n + t 2 + l)T(j + l)T(j + TO + n + 72 + 1)T(Z - r + m + r 2 + i)r(n + 1)A T 


m —0 


r(j + n + l)r(m + r 2 + l)(j - to)!to!T(Z - r + l + r 2 + n + 2) 


and 


A b,lr — 


(-I) b - Z r(6 + r 2 + i)r(6 + z + n + T 2 + I)r(i + Z) 
r(Z + r 2 + l)r(6 + n + r 2 + 1)(6 - Z)!Z!r(l + Z - r)A k ‘ 
Proof. The proof of this lemma is similar as Lemma 5.1.2. 


(5.1.51) 

(5.1.52) 

□ 


Lemma 5.1.4. Let $n 2 {x, y) be the function vector as defined in (5.1.17), then the mixed partial 
derivative of order t of ^^(x, y) is as follows 

dx (T/2)dyir/V (^ 2 fo V)) ~ ^N 2 xV &N 2 ( X ’ V)> ( 5 ' L53 ) 

where P> T fjf xN 2 is the operational matrix of mixed partial derivatives of order r, and is defined as 



( ^1,1 ,/,m 


^1 ,q,l,m 

^ l,N 2 ,l,m 


^2,1 ,Z,m 

^2,2,Z,m 

^2 ,q,l,m 

^ 2,N 2 ,l,m 

jy x ’V - 

^N 2 xN 2 ~ 


^ r,2,L,m 

^ r,q,l,m 

^r,N 2 ,l,m 


\ ^N 2 ,l,l,m 

^ N 2 ,2,l,m ’ ‘ 

^ N 2 ,q,l,m * ‘ 

^ N 2 ,N 2 ,l,m 


(5.1.54) 
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and q = Nj + i + 1, r = Nb + c+1, H r , q ,i,m = to j,i,b,c,l,m,T for j, i,b,c= 0,1, 2, n, 


A — ® AW ® A .X(j,t,l,m), 

, rx , rT , (6,i,Ti,T 2 ,A) (c,m,Ti,T 2 ,A) 

;=fjl m=r51 


(5.1.55) 


1 J * 

X(j,i,l,m) = T/f/ (ri,r 2 ) M/ (ri,T 2 ) ^ Q(i,m',n ,r 2 , A) Qfi.i'.n ,-r 2 .A) tym' ,t,Ti ,t 2 .A) ^Igi , m .Ti ,t 2 .AP 

(5.1.56) 


^7 ^17 m , =0 , =0 


0, 


(_i)j-m r (j + T2 + i)r(j + m > 


+ i) 


and 


(j,m ,ti,t 2 ,A) p( TO / + t- 2 + l)r(j + T\ + r 2 + 1)(6 — m'y.m'l A m ' ’ 

o _ r(m' + 72 + 1 - t/2 + i)r(n +1) . m ' +T2+/ _ T / 2 +ri+ i 

“(m'.J.n.Ti.A) - rK + r^l-T^ + n + l) 

(—l) &_i r(& + r 2 + l)r(6 + l + Ti + r 2 + 1) 


(5.1.57) 

(5.1.58) 

(b,i,rur 2 ,A) ~ r(/ + r 2 + 1)T(6 + n + r 2 + 1)(6 - 0!A'r(l + fc - r/2)' 1;1 1 

Proof. To prove the result take the fractional derivative of order r of Q^f£ 2 ^{ x ,y) as defined in 
(5.1.11), the following relation is obtained 

ax >JL(rm Q, ™P^ = ™,T) t, 

y i=o ™-n 


0 


m— 0 


where 


0 


(6,;,ti,t 2 ,T) 


(—l) b T(& + r 2 + 1)T(6 + l + n + r 2 + 1) 
T(Z + r 2 + l)r(6 + n + r 2 + 1)(6 - V)WA l ’ 




d x {T/2)dy( T /2) 


x (l-T/2)y(m-T/2) ^ 


Q A.n ( X 1 y) 'y ' 0(6,i,Ti,T 2 ,A) ^ 1 0(c,m,Ti,T 2 ,A) X 
z=m m=rii 

r(i + or(i + m) 

T(1 + l — r/2)T(l + m — r/2) 

The series approximation of a ;( (_T / 2 )?/ m_T / 2 ) with two dimensional JPs, yields 

n n 

= xo,MQZr > WQZr > (y), 

j=0 i=0 


(5.1.60) 


(5.1.61) 


where 


1 


X(j,i,l,m) i Ti ,t 2 )ti ^(n,r 2 ) 

W Aj ^A.i 


J j 0,0—(ar)^--) (y)^ 1 - 2) (or, 

(5.1.62) 

Above relation further implies that 


X(j,i,l,m) (t 1iT2 ) ( ri ,T 2 ) ®(j> m '> T i-' r 2,A)0(i,J' > ri,r 2 ,A) f / / ) ( A j/)idj/ 

"A,j ” A,i m'=(H'=0 J0 

/■A 

/ j;! 1 " +T 2 +m-r/ 2 )(A _ a .)^ c J a ,_ 

Jo 


(5.1.63) 
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Applying convolution theorem of Laplace transformation to the integrands of the above relation and 
then taking inverse Laplace of the obtained relations, we have 


[ x (m ' +T2+l - T/2) {A -x) 7 

JO 


dx 


r(m' + t 2 + i-t/ 2 + i)r(n + i)A( m '+ r2 + i - r / 2+Ti + 1 ) 
r(m' + r 2 + l — t /2 + n + 1) 


and 


2/ 


(i'+T^+m—r/2)( A _ yY^dy = 


r(r + r 2 + m - r /2 + l)r(n + l)Ad'+r 2 +m-r/2+r 1+ l) 
r(Z' + r 2 + m — r/2 + n + 1) 


With the help of above relations, the coefficient X(j,i,i,m ) of (5.1.61) can be written in a more compact 
form as 


1 3 

X(j,i,l,m) = w (T 1 ,r 2 ) TT/ (r 1 ,r 2 ) ® U,m',n ,t 2 ,A)©(»,/', Tl ,r 2 ,A)^(m',J,n,r 2 ,A)^(V ,m,Tx, t 2 ,A)> 


where 


and 


MA 1 ’ 2 'T/W 

"AJ "A,i m'=0i'=0 


o _ r (ra' + 72 +i - r/2 + i)r(n +1) m ' +T2+/ _ r/2+ri+1 

(m ' , ’ T 1 ’ T2 ’ A ) r(m' + r 2 + Z - r/2 + n + 1) 

ry _ r(/' + r 2 + m — r/2 + l)r(r! + 1) ;' +T2+m _ T / 2+Tl+1 

(i'.m.n.Ts.A) r(Z' + T2 + TO — r/2 + n + 1) 


(5.1.64) 

(5.1.65) 

(5.1.66) 


Now using (5.1.64) along with (5.1.61),(5.1.65) and (5.1.66) in equation (5.1.60), we get 


( g )Q£c (y) 

cte( T / 2 )<9y( T / 2 ) 


= EE E E e 


(6,Z,ri ,ro,A) (c.m.ri ,ro. A) 


X(i,i,J,m)QA,i T2) ( a; ) ( 3A n / T2) (2/)- 
(5.1.67) 


where 


and 


Let 


we get 


e 


(—l) ft T(6 + r 2 + l)r(6 + l + Ti + r 2 + 1) 
(b,/,Ti,r 2 ,A) r(z + t 2 + i)r(6 + t\ + t 2 + 1)(6 — 0!A*r(i +1 — t/ 2 ) 

(—l) a_m r(a + t 2 + l)T(a + m + n + r 2 + 1) 

(o,m,ri,r 2 ,A) r(m + T 2 + l)F(a + Ti + T 2 + l)(a — m)!A m r(l 4- m — r/2) 


w 3 ,hb,c,l,m,A ^2 ® AW ® 

, rT1 (b,i,ri,r 2 ,A) (a,m,ri,T 2 ,A) 

*=r§i m=r§i 


0 ' 


(5.1.68) 

(5.1.69) 

(5.1.70) 




cU’( T / 2 )<9y(r/2) 


Q^b ( X )Q A,civ) ™j,iAa,l,m,AQ ( AJ T2) ( X )QA )i (j/)' (5.1.71) 

j—0 i—0 


Using the notations q = Nj+i+1, r = Nb+c+1, D r ,q,i,m = j,i,c,b,i,m,T for i, j, b, c = 0,1, 2, 3,..., n, 
we get the desired result. □ 
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5.2 Applications of the operational matrices of integration 
and derivative of FOPDEs 

In this section we are interested to generalize a new numerical technique to approximate the solutions 
of a generalized class of multi terms FOPDEs and its coupled systems. 

5.2.1 Solution of multi-terms of FOPDEs by JPs 

Firstly, we develop the numerical technique for the multi-terms FOPDEs with mixed partial deriva¬ 
tives of fractional order of the form 


d ai U{x,y) _ d' 11 U(x,y) d ei U(x,y) 

dx° 2 ~ ai dtfn + ° 2 dx^dytn/* + G ' {x ' y >' 


(5.2.1) 


corresponding to the following initial conditions 


U^(0,y) = H j (y),j = 0,l...,m. 


(5.2.2) 


We are interested in finding the solution of (5.2.1) in the form of orthogonal two-parametric shifted 
JPs, such that the following exists 


d <7l U(x,y) 


= K JV 2$ Ar 2(a;,2/). 


(5.2.3) 


By application of fractional integral of order a and making use of operational matrix of integration, 
above expression can be written as 


U(x,y) = K N2 Q^ N A N *(x,y). 


(5.2.4) 


Using (5.2.2), we have 


U(x,y) = K N2 Q^ N A^(x,y) + F 1 N A N2 (x,y). 


(5.2.5) 


By application of operational matrices of derivatives, we have 


dll ^ V) = ^-N 2 Q G ’2x N 2 ^N^'xN 2 ^N 2 {x,y) + F 1 N2 i)^l’^ N 2&N2(x,y), 


d ei U(x, y) 
dx ei / 2 dy ei / 2 


— K - -.rW’U r»(ei>*,i/) 

— x\.n 2 *<4 N2 x N2 u N 2 x N2 


®n 2 (x,y) + ’*’*1 4 > jv 2 (x,y). 


(5.2.6) 


(5.2.7) 
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Now using the estimates (5.2.7),(5.2.6) and (5.2.4) in (5.2.1) we get the following system of algebraic 
equations 


Kjv 2 = + a 2 D^’ x x $) + + a 2 D^$) + F% 2 . (5.2.8) 

The unknown matrix Kjv 2 can be easily computed by solving the above system of algebraic equations. 

5.2.2 Solution of Coupled systems of multi-terms FOPDEs by JPs 

Now we develop the numerical technique for the following coupled system of multi-terms FOPDEs 
with mixed partial derivative terms of fractional order as 

d ai U(x,y) d 11 U{x,y) d l2 V(x,y) d l3 V(x,y) 


dx a1 01 dy ^ 


o 2 - 


d e *U{x,y) 

°4 ^ „./o + a 5 


dx T' 2 
d B3 V(x,y) 


a3 dy'r s ' + 


dx g 4 / 2 dy Q4 / 2 dx^/^dy 63 / 2 


+ Gi_{x,y), 


d° 2 V{x,y) d Pl V(x,y) d p2 U(x,y) d p3 U(x,y ) 

=0 1 7 —- h 0 2 - —— - h 0 3 - 


dx a2 


dx p2 


dy p 3 


dy pl 

d Pl U(x,y) d S2 V(x,y) 

4 dx ei / 2 dy ei / 2 5 dx e 2 / 2 dy P2 / 2 2 ’ ’ 


corresponding to the following initial conditions 


(5.2.9) 


C (j) (0, y) = hj(y ), F 0) (0 ,y) = gj(y), j = 0,1...,to, 


(5.2.10) 


where to < ay, cr 2 < m + 1, aj,bj , j = 1, 2,..., 5 are all real constants and U = U(x,y), V = V(x, y) 
are the unknown solutions of the system to be determined. Moreover U(x,y),V(x,y) G C([0, A] x 
[0, A]), and Fi(x, y), F 2 (x, y) € C([0, A] x [0, A]). We are interested in finding the solution of the 
problem (5.2.9) in the form of orthogonal shifted JPs, such that the following exists 

d^U(x,y) ^ d a2 V(x,y) A , , roilA 

— =K N 2 <P N 2 {x,y) ,— 7 ^ = L N 2$ N 2 (x,y). (5.2.11) 

Now applying integral of order a i and cr 2 on the corresponding equations, we get 

m 

U(x,y) -£c,V = K N 2 Q { x 2 'xN 2 $N 2 {x,y), 
j=o 

m 

V(x,y) -^2 djX 3 = Ljv 2 Q^ 2 xn 2 ®n 2 {x, y). 
j=o 


(5.2.12) 
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The generalized relations for the solution of the problem (5.2.9) and constants of integration can be 
obtained by applying the initial conditions (5.2.10) as 

U{x,y) = K N 2Qfy’* ] N 2 $ N 2(x,y)+ F 1 N 2$ N 2{x,y), 

(5.2.1.3) 

V(x,y) = L N 2Q^axlfi& N 2 (x, y) + F 2 N 2$ N 2 (x,y), 

where ( x )2/) = Ej=o h j(v) x °-, f 'n2 <&N 2 (x,y) = J27=o 9j(y) xj - For simplicity of notations 

we can write 

K n* Q&xlr* + F ^ 2 = f n 2 Qn 2 7n 2 + f1 n2 = O- (5.2.14) 

Further, it can also be written as 

U(x,y) = K N 2 <t N 2 (x,y), V(x,y) = L N 2 $ N 2 (x,y). (5.2.15) 


Using relation (5.2.15) and operational matrices of derivatives we can get the following estimates 

d^U{x,y) , d^V(x,y) ^ A(72 , x ) A 4 

-Kjv 2 D jv2xiv2 $jv2(2:,!/), — F N 2F > N 2 xN 2®N 2 {x,y), 

-&*&&**<** £££»> 

d ei U(x,y) _ .i. / \ d e 2 V(x,y) _ r \ 

Q xei / 2 dyei /2 - K N 2 F > N 2 xN 2^N 2 {x,y), dxe2/2gyg2 /2 - F N 2 F > N 2 xN 2^N 2 (x,y), 

Fi{x,y) = F %2® N 2(x,y), F 2 {x,y) = F%A N 2(x,y). 


(5.2.16) 


( K N 2<& N 2 (x, y) \ _ f fl i K^D^x^iv 2 ^,!/) \ / a 2 L N2 T)^’ 7 n 2 ^n 2 (x, y) \ 

\ F N 2 &N 2 {x, ) \ h Ljv 2 t)^l’^ N2 ^N 2 (x,y) J \ b 2 K N2 f)^V xN 2<&N 2 (x,y) ) 

+ I a 3 L N2 f)7 3 2' X N 2 ^N 2 (x,y) \ / a 4 K,v2 D^lv^.j/) \ + / F ^ 2 $ N 2 (x,y) \ 

V ) V b 5 L^t)^ 2 2 ’^ N 2(x,y) ) \ f n 2 ®n 2 (x,y) ) 

+ ( a 5 L N 2l)xl x ’N2$ N 2 (x,y) \ 

V h1&ii ( ° 1 a ™l* N *{x,y) ) 


(5.2.17) 
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K 


N 2 


J N 2 


where 


)Mk 


N 2 


J N 2 


K 


N 2 


K 


N 2 


K 


N 2 


JN 2 


J N 2 


J N 2 


K 


N 2 


J N 2 


A(7i,J/) 

ailJ N 2 xN 2 

On 2 xN 2 

( On 2 xN 2 
a2L, N 2 xN 2 


On 2 xN 2 
\(puv) 


N 2 x N 2 


o_ 

( „ rS (Bi,x,y) 
N 2 xN 2 

On 2 xN 2 

On 2 xN 2 

«5-L' Af 2 x JV 2 



'N 2 xN 2 

O n 2 xN 2 


■ N 2 


N 2 


A, 

(5.2.18) 


A = 


$N 2 (x,y) 0 N 2 \ 

V °n 2 ®N 2 (x,y) ) 

On 2 is a column vector of zeros, and Oj<[ 2 xn 2 is matrix having all entries equal to zero. 

Now canceling out the common term and after a short simplification we can write (5.2.18) as 


K 


N 2 


JN 2 


)-( K 


N 2 


2 N 2 


H -(F^ 2 F^ 2 ) =0, 


(5.2.19) 



„ fS(7i.y) 

UllJ N 2 xN 2 


„ fS(C4 ,x,y) 
U4lJ N 2 xN 2 


a2U N 2 xN 2 


n fS(73 ,V) 

a 3 U N 2 xN 2 


I „ f\(P3,x,y) 
+ U 5 1J N 2 XN 2 


u r\{P2,x) , f.{p 3 ,v) , , pv(ei ,x,v) 

b ? U N 2 xN 2 + b 3 U N 2 xN 2 + b4L) N 2 xN 2 

f\(B 2 ,x,y) , I j\(pi,y) 

°57V2 V V2 + 0iU N 2 xN 2 


N 2 x N 2 

Using the values of Kjv 2 and Ljv 2 we can write the above equations as 

( K N 2 L n 2 ) — h N 2 ^ Hi — (Fjy.F 2 ^) H — (F^ 2 F^ 2 ) 


= 0, 


(5.2.20) 


H, = 


where 

n(<ri.U r> 

'4n 2 xN 21J4 

,(<T1,X) 

l N 2 xN 2 

D 2 = 


Q \° 1 (C71 

N 2 xN 2 *~ J 3 ^2-^4 

,(P2,x) , fvO 

N 2 xN 2 + b 3 U N 

D 4 = ( 656 ^’*^ + b 1 t) ( f uv) 


i.U n \ 

^ N 2 XN 2U2 \ Di = (aif> (7i, B ) 


(<Tl,x) 

N 2 xN 2 

^ 3 ^N 2 xN 2 + *40$’*$), D 3 = (a 2 D) V 2 XJV 2 


N 2 xN 2 


a4lJ Af 2 xAf 2 l 


(ei ,*,!/)> 


»(72,x) 


„ fv(73,y) 

a 3- U Ar 2 xAT 2 


„ fdes.z.y) 
“S-L'AraxAf 2 


), and 


N 2 xN 2 ^ t N 2 xN 2 r 

The unknowns (Kjv 2 Ljv 2 ) can be easily determined using (5.2.20) which is a system of easily 
solvable algebraic equations. The approximate solution of the problem (5.2.9)-(5.2.10) can be easily 
computed by inserting (K N 2 L^ 2 ) into (5.2.15). 
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5.2.3 Numerical test of the proposed scheme 


By solving some test problems, the validity of the proposed method is examined in this section. 
Note that all the simulations are carried out using 5 Ghz processor. All the results are displayed 
using plots and tables. The algorithm is designed in such a way that it can be easily simulated using 
any computational software. We use MAT LAB for calculation and simulations. 

Example 5.2.1. As a first test case consider the following FOPDEs with mixed partial derivative 

tCVTflS 

d ai U(x,y ) d 7l U(x,y ) d ei U(x,y) 


= a i- 


a 2 


dx 2 dy 


t + G\{x, y), 


(5.2.21) 


dx ai dy 71 

where source term G\ is given as under 
Gi{x,y) = 3923 - z 2 cos (y) - 2;r cos (y)- 

Example 5.2.2. As a second test case assume the integer order coupled problem of PDEs of the 
form 


d 2 U dU dV dV r d 2 U „ d 2 V „ , , 

wy —2——b 2——b 3——b 5 — tt— b 6— —— \-Gi(x,y), 
dx z dy dx dy dxdy dxdy 

d 2 V dV dU dU n d 2 U d 2 V „ , , 

~A ~2 ^A —^ ^A —^ 2 A —^ ^A~A A~A —b G 2 (x,y), 

dx z dy dx dy dxdy dxdy 


(5.2.22) 


where the source terms G 1 and G 2 are defined as 

G± = — 8 x 4 y 3 — 80a; 3 y 3 + 15a; 3 y 2 + 12 x 2 y 4 + 6 x 2 y 3 +99 x 2 y 2 — 6 x 2 y — 6xy 3 —4xy 2 — 24a :y, 

(5.2.23) 

G 2 = — 8x 4 y 3 — 16a: 3 y 4 — 48 x 3 y 3 +15 x 3 y 2 +42 x 2 y 3 +18 x 2 y 2 — 6x 2 y—6xy 3 +4xy+2y 2 . (5.2.24) 
corresponding to the following initial conditions 

U (0, y) = U'(0, y) = 0, E(0, y) = V'( 0 , y) = 0. (5.2.25) 

Example 5.2.3. As a third test case consider the following coupled system of FOPDEs 


d 18 U _d°- 8 U d 2 U d 2 V 

dx 1 - 8 dy 0 - 8 dxdy dxdy 1 ’ 

d 18 V d 08 V d 2 U d 2 V 


(5.2.26) 


9a; 1 - 8 dy 0 - 8 dxdy dxdy 

where source terms G 1 and G 2 are given below 


+ G 2 {x,y), 


Gi (x, y) =3 x 2 y 3 {2y — 2) (x — 1) + 3 a; 3 y 3 (2 y - 2) (a: - 1) + 4 a: 4 y 3 {2x - 2) (y - 1) + 
3 a : 4 y 4 (2 x — 2) (y — l) 2 + 9 x 2 y 2 (x — l) 3 (y — l) 2 + 9 a: 3 y 2 (x — l) 2 (y — 1) 2 + 

16 a; 3 y 3 (x - l) 2 (y - l) 3 + 12 a; 3 y 4 (x - l) 2 (y - l) 2 - 
89181460669789625 a;T y 4 {y — l) 3 (125 a; 2 - 175 a; + 56) 

504403158265495552 + 

\ 2 


445907303348948125 x 4 y^ (x - l) 2 (4375 y 3 - 11625 y 2 + 10075 y - 2821) 


406548945561989414912 


(5.2.27) 
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G 2 (x,y) =3x 2 y 3 (2 y — 2) (x - l) 3 + 3 x 3 y 3 (2y - 2) (x — l) 2 + 12 x 4 y 3 (2x - 2) (y - 1) 3 + 
9 x 4 y 4 (2 x — 2) (y - l) 2 + 9x 2 y 2 (x - l) 3 (y - l) 2 + 9x 3 y 2 (x - l) 2 (y - 1) 2 + 

\3 

„ .. am i i £4-m innM / AMn/.n n 5 i r — i 

36x 3 y 4 (x — l) 2 (y — l) 2 + 


3 4 2 2 89181460669789625 x 3 yT (x-l ) 3 (125 y 2 - 210 y +84) 


3026418949592973312 


17836292133957925 x? y 3 (y - l ) 2 (1250 x 3 - 2625 x 2 + 1680 x - 308) 


1008806316530991104 
+ 48x 3 y 3 (x — l ) 2 (y — l) 3 . 

Example 5.2.4. As a fourth test case consider the following coupled system of FOPDEs 

d 18 U _dU 
dx 1 - 8 dy 

d 2 V d 08 U dU d 2 U 

dx 2 dx 0 - 8 dy dxdy 2 

where the source terms G\ and G 2 are given as under 

39239842694707435 xs e y 


dV d 2 V 

Tr x + 33Ty+ G ^ 


G\{x,y) =- 


18014398509481984 
G 2 {x,y) =y 2 e x — x 2 e y — 8xe v . 


-y 2 e x ~x 2 e y -2ye x 


(5.2.28) 


(5.2.29) 


(5.2.30) 

(5.2.31) 


5.2.4 Results and Discussions 


The exact solution for the first test Example 5.2.1 for a± = 1.8, 71 = 1, and Qi = 2 with a\ = 1, and 
a 2 = 2 is U(x,y) = x 2 sin{y). At first glance we compare the approximate solution of (5.2.1) with 
the exact solution of (5.2.1) by fixing a 1 = 1.8, 71 = 1, and Q\=2 with 77 = 2 = t 2 , the parameters 
of JPs at scale level M = 5. We observe that at scale level M = 5, our approximate solution are in 
a good agreement with the exact solution as displayed in Figure 5.1. We also investigate the error 
analysis at different scale levels and observe that the absolute errors between exact and approximate 
solution is minimizing as the scale levels are increasing. The results are displayed in Figure 5.2. We 
fix 71 = 1, Q\ = 2 and observe the behaviour of the approximate solutions of (5.2.1) for different 
values of 07 , it is found that the approximate solutions for different values of 07 approaching to the 
exact solution as 07 —> 3. The results are displayed in Figure 5.3. 

The exact solutions of Example 5.2.2 are given as under 

U (x, y) = (xy) 4 - (xy) 3 , V (x, y) = (xy) 2 - (xy) 3 . 

We approximate the solution of the second test Example 5.2.2 of integer order at different scale 
levels. We compare exact and approximate solution of the problem depicting in Figures (5.4) and 



102 


(5.5) for parameters ri = 2, 72 = 2 and scale level M = 5, observing that the exact solution matches 
well with the approximate solution. We also approximate the absolute error at different scale level 
M = 5,6, and 7 and at different points of the plane showing in Tables (5.1) and (5.2). Figure 

(5.6) showing that with the increase of the scale level the error of approximation (absolute error) 
decreases significantly. We also observe the behaviour of the numerical scheme for different values of 
<Ti and find that as <j\ —> 2, the approximate solutions approach to the exact solution. The results 
are displayed in Figures 5.7 and 5.8. 

The exact solution of Example 5.2.3 is given below that can be easily verified by direct substitu¬ 
tion 

U(x,y) = -x 4 y 4 (x- l ) 2 (y - l) 3 , V(x,y) = -x 3 y 3 (x - l ) 3 (y - l) 2 . 

We analyze the test Example 5.2.3 using our proposed numerical scheme and observe that the 
numerical results are in a good agreement with the exact solution. At different scale levels we check 
the validity of our numerical results and examine that our proposed numerical scheme exhibits great 
resemblance with the exact solution. We carry out the error analysis and note that with the increase 
in scale level the absolute error decreases in a great extent. We display the exact and approximate 
solution of (5.2.3) with parameters t\ = 1, T 2 = 1 and scale level M = 6. The results are displayed 
at Figures 5.9 and 5.10. The results regarding to error analysis with parameters t\ = 1, 72 = 1 and 
scale level M = 10 are displayed in Figures 5.11 and 5.12. 

The exact solution of fractional order test Example 5.2.4 is given below that can be easily verified 
by direct substitution 

U(x,y) = x 2 e y , V(x,y) = y 2 e x . 

We investigate the numerical accuracy of our proposed scheme with the help of Example 5.2.4 at 
different scale levels and examine that the proposed numerical scheme provides highly efficient and 
accurate results. We note that results generated by proposed scheme are highly accurate even for low 
scale level. That shows the accuracy of our developed scheme. We also do the error analysis and note 
that the absolute error decreases remarkably with the increase in scale level. The results obtained 
using the proposed scheme are compared with the exact solution of (5.2.4) choosing parameters 
ri = 1, and t -2 = 2 at scale level M = 7 and note that results are in a good agreement with the 
exact solutions as displayed in Figures 5.13 and 5.14. The graphical results regarding to the error 
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Figure 5.1: Comparing approximate solution with the exact solution of Example 5.2.1. 


Table 5.1: At various scale level M and for different values of x and y 1 the absolute error for Example 
5.2.2 is computed for the unknown function U(x,y). 


(x, y) \ M 

M = 5 

M = 6 

M = 7 

(x,y) 

= (0.i, 0.1) 

1.707 E - 13 

1A25E - 13 

1.567 E - 13 

(x,y) 

= (0.1,0.5) 

1.956E - 13 

1.596.E — 13 

1.235.E — 13 

{x,y) 

= (0.1,0.9) 

2.236E - 13 

1.773E — 13 

1.756^-12 

{x,y) 

= (0.5,0.1) 

1.873E- 12 

6.380.E - 13 

1.857 E — 13 

(x,y) 

= (0.5,0.5) 

1.880E - 12 

7.858 E - 13 

1.648^-12 

(x,y) 

= (0.5,0.9) 

2.092E - 13 

9.619F1 — 13 

1.829E- 12 

(x, y) 

= (0.9,0.1) 

1.270.E — 11 

9.363E - 12 

1.472E- 11 

{x, y) 

= (0.9,0.5) 

1.307 E - 11 

1.029.E 1 — 11 

1.497E - 13 

(x, y) 

= (0.9,0.9) 

1.306E- 11 

1.418E- 12 

1.652E - 13 


analysis with parameters t\ = 1, and 72 = 2 at scale level M = 10 are displayed in Figures 5.15 and 
5.16. 


5.3 Multi-terms partial differential equations of fractional 
order in high dimensions 

The main aim of this section is to extend the idea for high dimensional multi-terms FOPDEs. The 
new operational matrices are developed for a function vector involving three variables and with 
the help of these operational matrices, we solve the FOPDEs in high dimensions. A multi-term 
three dimensional heat conduction problem is chosen to check the efficiency and applicability of the 
proposed numerical scheme. 
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Figure 5.2: For various values of M (scale level), absolute error is computed for Example 5.2.1. 



Figure 5.3: For multiple values of <j\, the comparison is presented between approximate and exact 
solutions of Example 5.2.1. 






V(x,y) 
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Figure 5.4: Comparing approximate solution with the exact solution of Example 5.2.2. 



Figure 5.5: The comparison between approximate and exact solutions of Example 5.2.2 is presented. 
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Figure 5.6: Error of approximation in U(x,y) and V(x,y) of Example 5.2.2 


Table 5.2: At various scale level M and for different values of x and y, the absolute error for Example 
5.2.2 is computed for the unknown function V(x,y). 


(x, y) \ M 

M = 5 

CO 

II 

II 

(x, y) 

= (0.1,0.1) 

1.826.E - 13 

1.678E - 13 

1.231E- 13 

(x, y) 

= (0.1,0.5) 

2.638.E - 13 

1.839E - 13 

1.578E - 13 

(x, y) 

= (0.1,0.9) 

2.126E- 13 

2.019E - 13 

1.378E - 13 

(x, y) 

= (0.5,0.1) 

4.45LE - 12 

6.368E - 13 

1.582E - 12 

{x,y) 

= (0.5,0.5) 

3.781E- 12 

7.520E - 13 

1.673E- 12 

(x,y) 

= (0.5,0.9) 

3.933E - 12 

9.051E - 13 

1.567E - 11 

(x,y) 

= (0.9,0.1) 

1.650E- 11 

8.394E - 12 

1.482E- 12 

(x,y) 

= (0.9,0.5) 

1.820E- 11 

9.996E - 12 

1.785E - 13 

(x,y) 

= (0.9,0.9) 

1.746E- 11 

1.209E- 11 

1.652E- 13 














U(x,y) 
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Figure 5.7: For multiple values of a i, the comparison between approximate and exact solutions of 
Example 5.2.2 is done. 




V(x,y) 
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Figure 5.8: For multiple values of (J-\ , the comparison between approximate and exact solutions of 
Example 5.2.2 is done. 



Figure 5.9: Comparison is given between approximate solution and exact solution of Example 5.2.3 
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j jb'xact V ( x , 5 


Approximate V ( x , y) 



Figure 5.10: Comparison is given between approximate solution and exact solution of Example 5.2.3 



Figure 5.11: At scale level M = 10, the error of approximation is computed for U(x,y) of the 
Example 5.2.3. 
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Figure 5.12: At scale level M = 10, the error of approximation is computed for V(x,y) of the 
Example 5.2.3. 

High dimensional FOPDEs have numerous applications in the existing literature. Many en¬ 
gineering phenomenon are being modeled applying high dimensional FOPDEs, for example, the 
phenomena of cooling and heating of building structures, the phenomena of water reservoirs and 
heating lakes by radiations, the phenomena of heating of solid surfaces in materials processing, the 
phenomena of cyclic heating of laminated steel during pickling, and many more, see [19, 70, 83, 170]. 

The study of high dimensional FOPDEs is relatively very much complicated as compared to 
2-dimensional FOPDEs. Many researchers have devoted their intellect in the development of the 
numerical schemes of high dimensional FOPDEs. Kulish [116] developed a proficient numerical 
method based on laplace transform for the approximate solution of high dimensional integer order 
PDEs. M. Akbarzade [3] applied Homotopy analysis method to develop an approximate scheme for 
high dimensional transient state heat conduction equation. But the developed method [3] based on 
very much complicated and complex algorithms and is not easy to apply on multi-terms FOPDEs. 
Povstenko [153] applied the integral transform method to construct the axi-symmetric solutions for 
a time fractional heat conduction equation corresponding to Robin BCs. Currently, Wu [182, 183] 
studied fractional diffusion equation and developed the numerical scheme for it with the aid of 
variational iteration technique. The reader can also consult [4, 144, 154] to see more results on the 
development of the high dimensional FOPDEs. 
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For application purpose, we consider the following multi-terms three dimensional heat conduction 
equation in a generalized form as 


d ai U(t,x,y) d 01 U(t,x,y) d 02 U(t,x,y) d 03 U{t,x,y ) ^ ~ 

— m — = Cl —d^i — + C 2 ~d^ — + Cs 9^% + Gl(M ’ y) ’ (5 - 3 - 1} 

where C-\. C '2 and C 3 are constants, 0 < <j\ < 1, and t £ [0, A], x £ [0, A], y £ [0, A], subject to the 
initial condition 


U(0,x,y) = f(x, y). 


(5.3.2) 


The key idea is to transform the given problems into a system of easily solvable algebraic equations 
that can be easily handle by any computational software. We use the software Matlab for simulation 
purpose. We start the development of the scheme with the brief introduction of three dimensional 
shifted JPs. The notion of two dimensional shifted JPs has been discussed in Section 5.1.1. 


5.3.1 Shifted JPS in three dimensions 


The shifted three dimensional JPs of order N on the domain [0, A] x [0,A] x [0, A] can be expressed 
in terms of the product of function of three JPs as following 


QAMc(t, x ,y) = (Q£r\mQ£r\mQ { ir\y)u= o, 1,2,&=0,1,2,...,«, c =0,1,2, 

(5.3.3) 

its orthogonality condition can be determined as 

r A /*A rA 


III 

^0 Jo Jo 

R^’ T 2 \y)dtdxdy = W^/ 2 ) 5 a , d W^ T 2 ) 6 b , e W^f 2 ) 6 cJ . 


(5.3.4) 


5.3.2 Function approximation with three dimensional shifted JPs 

In this section, we assert on the idea that any function r(t, x, y) defined on [0, A] x [0, A] x [0, A] 
can be approximated in terms of the series of three-dimensional JPs as follows 

n n n 

r(t,X,y) -^^^ D abc{QA]a 2 \0){Q { L]b 2 \ X ))iQ { A]c 2 \y))- 

a—0 6—0 c—0 


(5.3.5) 
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The coefficients D a t c can be deduced using the following relation 

-j I'A p/\ 

D a bc 


W 


(i-l ,r 2 ) ^(n ,T 2 ) w (n ,r 2 


w\ 


r(t,x,y)(Q^\mQ ( Ar\^(QAT\y)) 


’A,a vr A,6 ^A.c 

R^ 1,T2 \t, x, y)dtdxdy. 


< 0 ^0 JO 


Where 


R ( £' T2) (t,x,y) = R { £’ T2 \t)R { £’ T2 \x)R { £' T2 \y). 
In terms of vector notation (5.3.5) can also be written as 

r{t,x,y) « ${t) T G N xN 2 &N 2 x 1) 


(5.3.6) 

(5.3.7) 

(5.3.8) 


where G represents the coefficient matrix and $ indicates the function vector related to the variables 
x,y , and <f>(t) shows the one dimensional Jacobi function vector related to the variable t. 


5.3.3 Error Analysis 


In this section we develop an analytical relationship for the error of approximation for sufficiently 
smooth function r(t, x, y) € C([0, A] x [0, A] x [0, A]). For this assume that Hat N N (t, x, y) the space 
of N terms JPs and g(N,N,N)(t,x,y) is its best approximation in II(Ar n n)(^i V)- Then according 

to the definition of the best approximation for any polynomial /( n,n,n ) (i> x , u) of degree < N in 
variable t, x and y, we can write 


II r{t,x,y) - r( N<NjN ){t,x,y )|| 2 < \\r(t,x,y) - f( N ,N,N){t,x,y)\\ 2 - (5.3.9) 

Then by following the same procedure as adopted in [76], we can find the bounds for the errors given 
below 

in ^ r n u ^ 1 / ( A “ °) ^+i t d N+ 1 r(t,x,y) 

\r(t, x, y) - f( N ,N,N)(t,x,y)\ < ^—) rnax^ iXtV) |- q x n+i -1 

i 1 1 (A — 0) s^jv+i 

+ T (—^—) max 


4 v N 


d N+1 r(t, x, y ). 1 (A-0) N+1 d N+1 r(t,x,y) 

d,x,y)\ dy N +1 l + N > ^ x ’ y '>'dt N + 1 (N+ 1)V 

_1 , (A - 0) ,,y -i. (A - 0) JV+1 . (A - 0) . jv+i | d 3N+3 r(t,x,y) 

64 N [ N ’ N (t’X’V) Q x N+iQ y N+ig t N+i 


1 

— max 


Q2N+2 


Mt,x,y) |. 


16 {t,x, v ydx N + 1 dy N + 1 

T i \ \ d N+1 r(t,x,y) \ -v \ d N+1 r(t,x,y) \ x \ d N+1 r(t,x,y) \ \ 

Let Ai = max (t ,x,y)\ — 9 x n+i U a 2 = max( tiXi y) \— ^+1 A 3 = max( t ,x,y)\ — dtN \\ A 4 = 

maX{t,x, y ) 1 9 x N+%lyN+igtN+% I, and A 5 = max (MiS) x -> V)\- Then, we have the following 
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relation 


\r(t,x,y) - f( N ,N,N)(t,x,y)\ < N °' > ) jV+1 Ai 

l,(A-0) yv +1 1 (A-0) JV +1 
+ 4 IV Aa 4 TV As 

1 / (A — 0 ) ^jv+i /(A — 0 ) , JV+1 , (A — 0 ) xAr+ix 1 / (A — 0 ) , N+ i ,(A — 0) ,jv +1 

- at - ) \ - AT - > \ - rr - ) A 4 + Tul- rr - > 1- at - > A 5' 


64 ' N ' K N ' N 

Finally using (5.3.9), we have the following relation 


16 v N 


N 


II r(t, x , y) - r {N , N>N) (t, x, y) || 2 < y/ (A - 0)(A - 0)(A - 0 )(|() N ^ 


N 

_1 / (A ~ 0) ,Af+i/ (A — 0) ,at + 1 / (A — 0) .jy +1 . 

CIA v AT / \ AT ) \ AT ) 4 


64 v N 
1 ,(A-0) N+1 .(A-0) N+1 , 
+ I6 ( N As ^ 


N 


N 


N 


(5.3.11) 


(5.3.12) 


5.3.4 Application of the operational matrices for generalized FOPDEs 

d ai U(t, x, y) d 01 U(t,x,y) d^U{t,x,y) d 03 U(t,x,y) 

— = Cl dxt>. + ° 2 V 3 aZ-tsfi l( • ’” ) ’ ' 1 

where C \, C% and C 3 are constants, 0 < o\ < 1. The time and space variables are defined as 
t € [0, A], x € [0. A], y £ [0, A]. Our main objective is the development of the numerical scheme 
subject to the following initial condition 


U(0,x,y) = f{x, y). 


(5.3.14) 


We seek the solution of the problem in terms of shifted JPs such that the following holds 

d^U(t,x,y) 


dt a i 


&{t)J N xN*®{x,y). 


(5.3.15) 


Applying fractional integration of order <ti with respect to t on (5.3.15), we have 

i ai dU( dP y) «$( *»$(*, y), 


which implies that 

U(t,x,y) - C!^(t)Ql^ M K MxM 2 ^{x,y). (5.3.16) 

The initial condition U(0,x,y) = f(x,y) yields c\ = f(x,y) and the above equation can be written 
as 

U(t,x,y)tt$(t)Qff xN J NxN 2 $(x,y) + f(x,y), ( 5 . 3 . 17 ) 
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which implies that 



U{t,x,y)^{t)(Q a ^ xN J NxN 2 + F NxN 2 )$(x,y). 

(5.3.18) 

Using the above estimate of U(t,x,y), we get 



q 0 ^^(^{QnxN^NxN 2 + F N x N 2 )G )( 1 >(x,y), 

(5.3.19) 

and 

dp2U ^ x,y) ^mmxNJ NxN 2+F NxN 2)G^ y) H x ,y)- 

dy 2 

(5.3.20) 

and 

dhU J? ,X £ ] ~Ht)(QNx N J N xN 2 + F NxN 2 )G&’ x ’ri*(x,y). 
dx 2 dy 2 

(5.3.21) 


Gi{t,x,y) = ${t)G NxN 2 $(x,y) 

(5.3.22) 


Using (5.3.15), (5.3.19), (5.3.20) in (5.3.13) we get 


&(t)JNxN 2 &( x > y) = Gi^(t)(Q^ xN J NxN 2 + F NxN 2)G < ' l3l ' a/> &(x, y) 

+ C 2 &(t)(Q^ xN J NxN 2 + F NxN 2)G^ 2 ’ y ^(x, y) 

+ C 3 ^(t)(Q^ xN J NxN 2 + F NxN 2)G^ 3 ' x,v ^(x, y) + <&{t)G NxN 2<b(x, y), 

( 5 . 3 . 23 ) 

which can be rewritten as 

^(t)(Cl(QNxN^NxN 2 + FnxN^G^ 31 '^ + C2{Q^ xN JnxN 2 + x N 2 )G^ 2 '^ + 

C 3 {Q^ xN Jnxn 2 + Fnxn 2 )G^ 3 ' x ' v ^ — Jnxn 2 + G NxN 2)&(x, y) = 0 . 

Hence it follows that 

Ci{Q^ n Jnxn 2 + F NxN 2)G^ + C 2 (Q^ n J N xN 2 + F NxN 2)G&’V) + 

( 5 . 3 . 24 ) 

C 3 {Q a ^ xN JNxN 2 + FnxN 2 )G^ 3,V ^ — JnxN 2 + G NxN 2 = 0, 

Which can be written as 

A-Jnxn 2 ^ ~~ Jnxn 2 = C ( 5 . 3 . 25 ) 

Which is lyapanove type matrix equation, where 


A — O aiA 

~ ^NxN 
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B = CiG (/3l>a:) + C 2 G i02 ’ v) + C 3 G i03 ’ x ’ y) , 


and 


C — F NxN 2B + G nxn 2 , 

which is an algebraic equation of Sylvester type and can be easily solved for the unknown matrix 
Jnxn 2 using the value of J in (5.3.18) we can get the approximate solution of the problem. 

5.3.5 Numerical test of the proposed scheme 


In this section, some test problems are considered to examine the validity of the proposed method. 
Note that all the simulations are carried out using 5 Ghz processor. All the results are displayed 
using plots. The algorithm is designed in such a way that it can be easily simulated using any 
computational software. We use MatLab for calculation and simulations. 


Example 5.3.1. As a first test case consider the following multi-terms FOPDEs with mixed partial 
derivative terms 


d ai U(t,x,y) d^U{t,x,y) d 02 U(t,x,y) d 03 U{t,x,y) 


df 7 ' 


= C 1 


<9aA 


C 2 


dyP 2 


+ C 3 


0 ^3 o ^3 

OX 2 ay 2 


+ Gi(t, x, y), (5.3.26) 


corresponding to the following initial condition 


U(0,x,y) = x 3 y 3 , 


(5.3.27) 


where 0 < cy < 1, t £ [0,1], x £ [0,1], and y £ [0,1], The source term G\ is defined as under 

G\ ( t,x,y) = 3 1 2 — 6 x 3 y — 9 x 2 y 2 — 6 x y 3 (5.3.28) 


Example 5.3.2. As a second test case consider the following multi-terms FOPDEs with mixed 
partial derivative terms 


d ai U(t,x,y) = d pl U(t,x,y) d 02 U(t,x,y) 
dt a i 1 Ox?' 2 By? 2 


T C‘- 


d 03 U(t,x,y) 

O ^ 3 . o ^3 

OX 2 oy 2 


Gift, x 7 y), 


(5.3.29) 


corresponding to the following initial condition 


U{0,x,y) 


— e 2x+3y^ 


(5.3.30) 


where 0 < cy < 1, t £ [0,1], x £ [0,1], and y € [0,1]. The source term G\ is considered to be zero 
here. 
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5.3.6 Results and Discussions 

The exact solution of the first test Example 5.3.1 for a\ = 1 is given by t 3 + x 3 + y 3 . We compare 
the exact and approximate solution of the Example 5.3.1 at different values of t by taking scale level 
M = 5 applying our new method. We observe that our approximate solution is in a good agreement 
with the exact solution of Example 5.3.1 at t = 0.1, t = 0.4, t = 0.7, and t = 1. The results are 
shown in Figure 5.17. We also check the behaviour of the approximate solution of Example 5.3.1 
by fixing ( x,y ) = (0.1,0.1) and at some fractional values of o\. We observe that as a-\ —> 1, the 
approximate solution of Example 5.3.1 approaches to the exact solution. The results are displayed 
in Figure 5.18. We also do the error analysis of Example 5.3.1 at different values of t. x and y. We 
note that even at very low scale level M = 5, the absolute error between exact and approximate 
solution of Example 5.3.1 is very negligible. The results are plotted in Figures 5.19-5.21. 

The exact solution of Example 5.3.2 is e 2x+3y +At . Using our new numerical scheme, we compare 
the exact and approximate solution of Example 5.3.2 at different values of t. It is to be noted that 
the exact and approximate solution is in a good agreement at t = 0.1, t = 0.4, and t = 1. The 
outcomes are displayed in Figure 5.22. At some fractional values of (J \, the exact and approximate 
solution of Example 5.3.2 is compared by fixing ( x , y) = (0.5,0.5). It is evident that our approximate 
solution approaches to the exact solution as a\ —>• 1. The results are shown in Figure 5.23. The 
error analysis of the Example 5.3.2 is also done at different scale level and at different values of t. 
It is observed that as the scale level increases the absolute error between exact and approximate 
solution of Example 5.3.2 decreases in a great extent. The results are displayed in Figures 5.24-5.25. 
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Figure 5.13: Comparison is given between approximate solution and exact solution of Example 5.2.4 



Figure 5.14: Comparison is given between approximate solution and exact solution of Example 5.2.4 
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Figure 5.15: Error of approximation in U(x,y) of Example 5.2.4 



Figure 5.16: Error of approximation in V(x,y) of Example 5.2.4 











U{U,x,y) 
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'Exact U(t,x,y) att = 0.1 
Approximate U(t,x,y) att = 0.1 
lExact U ( t , x, y) att = 0.4 


Figure 5.17: At multiple values of t , approximate solution and exact solution of Example 5.3.1 are 
compared by choosing scale level M = 5, and JPs parameters t\ = 2 = r 2 . 
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Figure 5.18: For various values of <j\ and by fixing (x,y) = (0.3,0.3), comparison of approximate 
solution with exact solution is presented for Example 5.3.1. 





\U(t,x t ,y) - U{t,Xi,y)\ 
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Figure 5.19: Choosing scale level M = 5, and for various values of t, the absolute error for Example 
5.3.1 is computed. 



Figure 5.20: Choosing scale level M = 5, and for various values of x, the absolute error for Example 
5.3.1 is computed. 
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Figure 5.21: Choosing scale level M = 5, and for various values of y, the absolute error for Example 
5.3.1 is computed. 



Figure 5.22: At multiple values of t, approximate solution and exact solution of Example 5.3.2 are 
compared by choosing scale level M = 9. 
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Figure 5.23: For various values of a± and by fixing ( x,y ) = (0.5,0.5), comparison of approximate 
solution with exact solution is presented for Example 5.3.2. 



Figure 5.24: At various scala level M , and at t = 0.5, absolute error is computed for Example 5.3.2. 
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Figure 5.25: At various scala level M, and at t = 0.8, absolute error is computed for Example 5.3.2. 






Chapter 6 

Summary and Conclusion 


This thesis comprises of two parts. The first part addresses the existence results for the nonlinear 
ordinary coupled BVPs subject to nonlinear CBCs using LUSs approach. This part unifies the 
treatment of the existence results of many important linear and nonlinear BVPs such as periodic, 
anti-periodic, Dirichlet, and Neumann by applying some nronotonicity assumptions on the arguments 
of the certain type of BCs . The second part provides the efficient and reliable numerical methods 
to find the approximate solutions of the multi-terms FOPDEs and their coupled systems. In this 
part, the two-parametric orthogonal shifted JPs are used to establish the new operational matrices 
of fractional order integrals and derivatives. By means of the operational matrices, the considered 
FOPDEs are reduced to an algebraic ones. Being easily solvable, the associated algebraic systems 
lead to finding the solutions of the considered FOPDEs. This part of the thesis also deals with 
the development of the numerical schemes for the high dimensional multi-terms FOPDEs. The 
theoretical results included in the both part of the thesis are verified by taking many test examples. 

In Chapter 1, a comprehensive introduction on coupled LUSs approach is discussed along-with a 
short historical note on the subject fractional calculus. Some concepts form approximation theory 
and matrix theory are also studied in this chapter. Some classical results from Fixed Point Theory 
are included that are much useful to establish the existence of solutions of the considered nonlinear 
BVPs. Some properties of LUSs are also discussed. 

In Chapter 2, with the aid of Schauder’s fixed point theorem, we ensure the existence of at least 
one fixed point of a continuous and compact map, L~ l N : C' 1 [0,1] x C ,1 [0,1] —> C^O, 1] x C^O, 1]. 
Under the approach of LUSs, we ensure that the fixed point of the map L _1 N, which is a solution of 
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the first-order nonlinear coupled BVPs lies in a closed set defined by well ordered lower and upper 
solutions. We unify the treatment of the existence results of very important BVPs, for example 
periodic and anti-periodic by imposing some monotonicity assumptions on the arguments of the 
nonlinear boundary functions and using the idea of coupled LUSs. We observe that our developed 
results are very much generalized because they satisfy the classical existence criterion of many linear 
and nonlinear BVPs, like periodic and anti-periodic BVPs. The existence results are also studied 
when LUSs are in the reverse order. Some examples are studied to ensure the validity of our 
developed theoretical results. 

In Chapter 3, we extend the results of the Chapter 2 for the second-order nonlinear ordinary 
coupled BVPs with nonlinear CBCs. The existence of at least one fixed point of the continuous and 
compact map, L~ 1 N : C 2 [ 0,1] x C 2 [0,1] —> C 2 [ 0,1] x C 2 [ 0,1] is ensured by Schauder’s fixed point 
theorem, and the compactness of the map, L _1 N is claimed under the outcomes of Arzela-Ascoli 
theorem. We present the general approach of existence of solutions of many important linear and 
nonlinear BVPs, including periodic, anti-periodic, Diriclilet and Neumann in a unified way, which 
avoids treating these problems on a case-by-case basis. Several examples are treated to ensure the 
validity of the developed theoretical results. 

In Chapter 4, we extend the results of Chapter 3 for the nonlinear ordinary second-order coupled 
BVPs with dependence of the nonlinear functions on the first-order derivative. Under the approach 
of LUSs method, the existence of at least one solution is ensured in the closed set defined by well 
ordered LUSs. We also ensure that the nonlinear functions having derivative dependence satisfy 
the Nagumo condition relative to the closed intervals defined by ordered LUSs. Some examples are 
studied to ensure the validity of the developed theoretical results. 

It is very much evident that exact solutions of every fractional order problem are not always 
possible to determine. Therefore, there is a need to develop the reliable and efficient numerical 
methods that exhibit the very close agreement of the approximate solutions with the exact solutions 
of the considered fractional order problems. In this thesis we develop the efficient and reliable 
numerical methods based on operational matrices of fractional order integrals and derivatives to 
find the approximate solutions of multi-terms FOPDEs. The handsome amount of the test examples 
are taken to ensure the accuracy and efficiency of the proposed methods at different scale level. It is 
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to be noted that with the increase in scale level the absolute error between exact and approximate 
solution is decreased in a considerable extent. Also we check the behaviour of the proposed method 
at different fractional values of the order of the concerned FODPs and it is very much clear that 
as the order of the concerned problem approaches to the integer order, the approximate solution 
approaches to the exact solution with zero error. All the results have been shown in the form of 
Tables and Plots. The numerical analysis of the above discussion is elaborated in Chapter 5. 



Appendix 


Matlab Codes For Simulation of the Operational Matrices 
Algorithm to compute function vector: 

Function vector 4^'^ of orthogonal shifted JPs with parameters n and r 2 is generated by 
using the following algorithm: 

function K=funvec (t ,M, taul , tau2 , Delta ); 
syrns 1 taul tau2 t Delta; 
for j =0:M— 1; 

al =( —1-) (j — k ); 

a2=gamma( j+tau2 +1); 
a.3=gainma (j +k+t aul+tau2 + l); 
d l=gamma (k+t au 2 +1); 
d2=gamma (j +t au 1+t au2 + 1 ); 
d3=gamma (j —k +1); 
d4=gamma (k +1); 
d5=Delt a" k ; 

BB=(al* a2 *a3)/(dl*d2*d3*d4*d5); 

K(j+1)= symsum( BB* t' k , k , 0 , j ); 

end 

K=transpose (K) 

Where “M” is the size of the matrix and known as scale level of the numerical scheme. The 
parameters “taul, tau2” indicates the Jacobian parameters (ri and t 2 ), and “Delta” indicates the 
right end pint of the domain set [0, A]. 

Algorithm to generate weight function: 

The weight function is generated by the use of the following code: 
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function R=weight (t , taul , tau2 , Delta) 

R=t' tau2 * (Delta—t) ~ taul 

Where the parameters “taul, tau2” indicates the Jacobian parameters (t\ and r 2 ), and “Delta” 
indicates the right end pint of the domain set [0, A]. 

Algorithm to generate orthogonality constant: 

The orthogonality constant W^ 1 j’ T2 ' > is generated by the use of the following algorithm: 

function D=ort (i , taul , tau2 , Delta ); 

al=Delta ~ (taul+tau2 +1); 

a2=gamma(j +t au 1 +1); 

a3=gamma (j +t au 2 +1); 

dl = (2* j+taul+tau2 + l); 

d2=gamma(j +1); 

d3=gamma( j+taul+tau2 +1); 

C=(al*a2*a3)/(dl*d2*d3) 

Where the parameters “taul, tau2” indicates the Jacobian parameters (ly and r 2 ), and “Delta” 
indicates the right end pint of the domain set [0, A]. 

Algorithm to generate coefficient vector: 

The coefficient vector D ^ can be generated using the following algorithm: 

function F=coefvec (f,t ,M, taul ,tau2 ,Delta) 
syms t 

R= weight (t , taul , tau2 , Delta ); 

K=funvec (t ,M, taul , tau2 , Delta ); 

Kw=f *K*R; 
for j =0:M— 1; 

D=ort (j , taul , tau2 , Delta ); 
aa=KR( j +1); 

F( j+1) = (1/D)* int ( aa ,t,0,Delta); 

end 

Where the parameters “taul, tau2” indicates the Jacobian parameters ( t\ and r 2 ), and “Delta” 
indicates the right end pint of the domain set [0, A]. 
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Algorithm to generate operational matrix of fractional order integral: 

The operational matrix of fractional integral Q p ^ x M i of order p is generated with the aid of the 
following algorithm: 

function Q=j acobinta (rho ,M, taul , tau2 , Delta) 
syms p q 
for c=0:M—1; 

for e = 0 :M— 1; 

gg = ((— 1) * (c—q)) * gamma (c+t au 2 +1) * gamma (c+q+t au 1+t au2 +1); 
qq= gamma (q+t a u 2 +1) * gamma (c+t a u 1 +t au 2 +1) * gamma (c —q+1) * gamma (q+r h o +1); 
nlnl = (( — 1)' (e—p)) * gamma (e+p+taul+tau2+l)*gamma( taul +1); 
n2n2=gamma(p+rho+tau2 + 1)*(2*e+taul+tau2 + l)*gamma(e + l)*(Delta~rho ); 
nn=nlnl*n2n2 ; 

mlml=gamma( e+taul + I)*gamma(p+tau2 + l)*gamma( e—p + l)*gamma(p +1); 
m2m2= *garnma (q+p+t au 1 +t au 2+r ho + 2); 
mnr=mlml*m2m2; 

Q( c + l,e+l)=symsum ((( gg / qq) * symsum (( cc /mm) ,p,0 , e )) ,q,0 ,c); 

end 

end 

Where “M” is the size of the matrix and known as scale level of the numerical scheme. The 
parameters “taul, tau2” indicates the Jacobian parameters and 72 ), and “Delta” indicates the 
right end pint of the domain set [0, A], Moreover p represents the order of the integral. 

Algorithm to generate operational matrix of fractional order derivative: 

The operational matrix of fractional derivative °f or<; l er a is generated by using the 

following algorithm: 

function D=jacobder (sigma ,M, taul , tau2 , Delta ) 

syms p q 

bba=ce i 1 ( sigma ); 

for c=0:M—1; 

for e = 0 :M— 1; 
if c<bba; 

D( j +1, i +1)=0; 
else ; 

aa = (( — 1)" (j — 1 )) * gamma (j+tau2+l)*gamma( j+l+taul+tau2 + 1); 
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bb= gamma(l+tau2 + l)*gamma( j+taul+tau2 + l)*gamma( j — 1 + l)*gamma( 1— sigma + 1) 

clcl = (( — 1) * (i— k)) * gamma (i+k+taul+tau2+l)*gamma(a + l); 

c2c2 = (2* i+taul+tau2 + l)*gamma( i +1)*( Deltasigma )); 

c3c3=gamma(k+l—sigma+tau2 + 1); 

cc=clcl*c2c2*c3c3; 

dldl=gamma( i+taul + I)*gamma(k+tau2 + l)*gamma( i—k + 1); 
d2d2=gamma(k+l)*gamma( l+k+taul+tau2— sigma+ 2); 
dd=dldl*d2d2; 

D( j +1, i+l)=symsum ((( aa/bb) *symsum ((cc/dd),k,0,i)),l, bba , j ); 
end 

end 

end 

Where “M” is the size of the matrix and known as scale level of the numerical scheme. The 
parameters “taul, tau2” indicates the Jacobian parameters {ji and 72 ), and “Delta” indicates the 
right end pint of the domain set [0, A], Moreover, <7 indicates the order of the derivative. 

The codes for two and three dimensional function vectors along-with the rest of the other codes 
used in this thesis are left as an exercise for the reader. 
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